This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What Are Some Sanity Checks That Should Be Performed On Ngs Data?

I am compiling a wish list of analyses that should be run for every data set coming out of a sequencing facility - regardless of whether this sequence is for RNA-Seq, SNP calling, ChIP-Seq, or possibly de novo sequencing. The goal is to scan for potential red flags that would possibly indicate something has gone awry either in the lab or downstream. I want a list of "sanity checks" that will encompass both sequence quality analysis as well as what can be gleamed from alignments.

For example,

  • Sequence QA - basecalling bias, read quality, yield, throughput, GC bias, 5'/3' motifs?, restriction enzyme bias

  • Barcode distribution (if barcoded)

  • Alignment QA

    • chromosome bias,
    • annotational biases (whether experimentally induced or not)
      • genes, repeats, cpg islands, epigenetic markers, expression

I am sure this has already been implemented at a lot of the bigger sequencing cores - I just need a definitive list. Of course, many of these sanity checks will be triggered by the experiments themselves - the point is to develop a comprehensive checklist of analyses that will encompass both what we expect to see as well as what we don't.

next-gen sequencing pipeline quality

great question! also interested to see strand bias relative to annotations for RNA-Seq.

I suppose a kmer analysis for every dataset would not be unreasonable

@Bio_X2Y, that hexamer bias is visible in the FastQC output.

@brentp - thanks! I was aiming to highlight that platform-specific biases exist, I just used this as an example because I don't know of any others :)

great program, although I think raw tabular output would be welcome to developers

3 answers

We use FASTQC to perform a barrage of quality checks - you might get some useful ideas there.

We also quantify the amount of rRNA reads in our Illumina GA datasets - we hope to see around 4-6%.

+1 for FASTQC, it's the starting point for all our analyses.

People landing here, check out multiqc - it works with fastqc to make a nice combined report for all reads in a directory.

Yes I am seeing a few checks i hadn't listed as well as some great ideas for visualization -sequence length distribution -sequence duplication levels -overrepresented sequences

In addition to running FASTQC on every lane of sequencing, in my mapping pipeline I record the number of

  • Raw (purity filtered) reads
  • Unmappable reads
  • Multimapping reads
  • Uniquely mapping reads
  • Final number of reads after removing duplicates with Picard

These metrics tell us a few things eg certain types of experiments you expect to have more multimapping reads (DNA methylation pull downs), and the % of reads which are removed as duplicates really goes up when we're scraping the bottom of the tube when it comes to how much template we manage to get into library prep. Of course interpretation of these numbers really depends on the biological experiment going on.

I'm working with viral samples and I've running against this duplicate read problem. Do you know of a reference where this low concentration of DNA and high duplicate read count is discussed? Many thanks.

I definitely agree with the above post. Another very useful metric is the library complexity...this can be generated by running Picard's MarkDuplicates tool.

I agree that library complexity is an important metric. However, simply counting duplicates is perhaps an overly simplistic assessment of library complexity. At the very least the number of duplicates should be relative to the total number of reads in the library...

Log in to answer this question.