Answers.Body: <p><a href="http://multiqc.info/" rel="noreferrer">MultiQC</a> can merge all your different reports into a single one.
  Which could be useful once you manage to know which QC tools to use.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-10T15:18:41.337Z
Answers.Text: FWIW, Multiqc is also really pretty...
Answers.UserDisplayName: winni2k
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-01T09:56:56.927Z
Answers.LastActivityDate: 2017-06-01T09:56:56.927Z
Answers.UserDisplayName: MaxUlysse
Answers.Body: <p>The quality control of ngs reads is heavily dependent on type of the project. </p>
  
  <p><strong>For genome assembly</strong> projects based on short reads, beside already covered checking quality of sequencing, you would like to look at the kmer spectra to find out, if your reads are going to make sense when they will be translated to De Brujin graph. You will get also a clues about genome coverage, genome size, repetitive content and small clue about heterozygocity. A lot of useful info about interpretation you can find in the README of GenomeScope.</p>
  
  <p>A list tools I used:</p>
  
  <ul>
  <li><a href="http://www.genome.umd.edu/jellyfish.html" rel="noreferrer">jellyfish</a> - to count k-mer frequencies</li>
  <li><a href="https://github.com/schatzlab/genomescope" rel="noreferrer">GenomeScope</a> - a package for analysis and visualisation of k-mer frequencies (they recommend to use jellyfish for counting k-mer frequencies)</li>
  <li><a href="http://kmergenie.bx.psu.edu/" rel="noreferrer">kmergenie</a> - for prediction of optimal kmer for assembly</li>
  </ul>
  
  <p>Using these tools can save you a lot of frustration if you accidentally sequence a contaminated sample or if your colleagues / a sequencing company have sent you a wrong file!</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-04T17:58:23.480Z
Answers.LastActivityDate: 2017-06-04T17:58:23.480Z
Answers.UserDisplayName: Kamil S Jaron
Answers.Body: <p>We routinely run both <a href="http://www.bioinformatics.babraham.ac.uk/projects/fastqc/" rel="noreferrer">FastQC</a> and <a href="http://www.bioinformatics.babraham.ac.uk/projects/fastq_screen/" rel="noreferrer">FastQ Screen</a> on all of our raw sequencing reads. FastQ Screen is a tool for detecting cross-species contamination. <a href="https://github.com/crukci-bioinformatics/MGA" rel="noreferrer">MGA</a> is another similar tool.</p>
  
  <p>There are then lots of QC tools specific to different types of data, most of which run after alignment. For example <a href="http://rseqc.sourceforge.net/" rel="noreferrer">RSeQC</a> (RNA data), <a href="http://qualimap.bioinfo.cipf.es/" rel="noreferrer">Qualimap</a> and many many others. Without specifying what kind of data you have this is a bit difficult to make recommendations for though.</p>
  
  <p>Phil</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-03T09:06:32.320Z
Answers.LastActivityDate: 2017-06-03T09:06:32.320Z
Answers.UserDisplayName: ewels
Answers.Body: <p>The kind of QC you do routinely depends on what your lab's focus is. We do a lot of low-quality, multiplexed DNA and RNA. If you routinely do fresh frozen whole genomes, your QC will be different.</p>
  
  <p>Weighing in from the resequencing side of things (i.e. sequencing an organism that has a good reference), there are several things we are testing with quality control:</p>
  
  <ul>
  <li>Did the library preparation succeed? Was the capture successful? Is our input material of high enough quality? </li>
  <li>Did the sequencing work? For Illumina, did enough clusters pass filter? What's our Q30? Did we sequence into adapters?</li>
  <li>Can we proceed with analysis? Are the indices on multiplexed samples correct?
  Do we have the correct reference sequence? Do we have the right targets?</li>
  </ul>
  
  <p>You may use a combination of tools:</p>
  
  <ul>
  <li><a href="https://support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.html" rel="nofollow noreferrer">bcl2fastq</a> : if you multiplexed, check the undetermined indices file after converting to fastq to make sure you caught all of the indexed reads</li>
  <li><a href="https://www.bioinformatics.babraham.ac.uk/projects/fastqc/" rel="nofollow noreferrer">FastQC</a>: General sequencing quality, sequencing into adapters, handy flagging for obvious problems</li>
  <li><a href="https://github.com/Illumina/interop" rel="nofollow noreferrer">Illumina interop</a>/SAV: Cluster pass filter, Q30</li>
  <li>Alignment (with <a href="http://bio-bwa.sourceforge.net/" rel="nofollow noreferrer">BWA</a>, bowtie, etc), <a href="http://bedtools.readthedocs.io/en/latest/" rel="nofollow noreferrer">bedtools</a>, <a href="http://www.htslib.org/" rel="nofollow noreferrer">samtools</a>, <a href="https://broadinstitute.github.io/picard/" rel="nofollow noreferrer">Picard</a>: check for contamination, make sure we're using the right reference, overlap with target regions, depth of coverage, soft-clipping for low-quality sequence</li>
  </ul>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-06T15:04:52.107Z
Answers.LastActivityDate: 2017-06-06T15:04:52.107Z
Answers.UserDisplayName: morgantaschuk
Answers.Body: <p>If you are looking for NGS QC for your fastq, bam, bed and vcf files I would suggest a commercial tool called <a href="https://om.euformatics.com/apex/f?p=118:1:" rel="nofollow noreferrer">omnomicsQ</a>. </p>
  
  <p>It automates coverage analysis and sample pass, warn fail according according to your defined SOPs. Unlike fastqc it also comes with a database so that you can compare protocols, samples, and your performance to peer organisations. It also charts performance over time, exceptions,and  correlations between metrics</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-20T16:44:39.103Z
Answers.LastActivityDate: 2017-06-20T16:44:39.103Z
Answers.UserDisplayName: axthomas
Body: <p>What are good means for performing quality control (QC) or NGS reads?</p>
  
  <p>I'm aware of:</p>
  
  <ul>
  <li>FastQC</li>
  <li>NGS Screen</li>
  <li>Kraken (e.g., for screening against contaminants)</li>
  </ul>
  
  <p>What are other popular means for such QC?</p>
  
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-06-01T02:05:13.657Z
Comments.Text: This question is too broad. There are many available QC programs which all have their own advantages/disadvantages. Do you have a specific use case that you would like the QC to be applied to?
Comments.UserDisplayName: gringer
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-06-04T10:30:49.507Z
Comments.Text: agree with @gringer. Right now the answer to this question is, roughly, whatever you get by googling 'NGS QC'. If you can provide a detailed use case - e.g. maybe you have 200GB of paired-end illumina data from 12 individuals of a non-model species, and you have a distant reference available - it will be a much more useful question and one that google couldn't solve.
Comments.UserDisplayName: roblanf
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-06-08T04:32:16.280Z
Comments.Text: Agreed. I'd be happy to rescind my downvote if the question text were updated.
Comments.UserDisplayName: Daniel Standage
Tags: <ngs><quality-control>
Title: QC measures for NGS sequencing
ThreadText: # QC measures for NGS sequencing
  
  (Asked by: username_0 on 2017-05-31)
  
  What are good means for performing quality control (QC) or NGS reads?
  
  I'm aware of:
  
  * FastQC
  * NGS Screen
  * Kraken (e.g., for screening against contaminants)
  
  What are other popular means for such QC?
  
  **Tags:** ngs, quality-control
  
  > **Comments:**
  > **username_1**: This question is too broad. There are many available QC programs which all have their own advantages/disadvantages. Do you have a specific use case that you would like the QC to be applied to?  
  > **username_2**: agree with @username_1. Right now the answer to this question is, roughly, whatever you get by googling 'NGS QC'. If you can provide a detailed use case - e.g. maybe you have 200GB of paired-end illumina data from 12 individuals of a non-model species, and you have a distant reference available - it will be a much more useful question and one that google couldn't solve.  
  > **username_3**: Agreed. I'd be happy to rescind my downvote if the question text were updated.  
  
  ---
  
  ## 5 answers
  
  ### Answer by username_4 on 2017-06-01. Score: 12
  
  [MultiQC](http://multiqc.info/) can merge all your different reports into a single one.
  Which could be useful once you manage to know which QC tools to use.
  
  > **Comments:**
  > **username_5**: FWIW, Multiqc is also really pretty...  
  
  ---
  
  ### Answer by username_7 on 2017-06-03. Score: 7
  
  We routinely run both [FastQC](http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and [FastQ Screen](http://www.bioinformatics.babraham.ac.uk/projects/fastq_screen/) on all of our raw sequencing reads. FastQ Screen is a tool for detecting cross-species contamination. [MGA](https://github.com/crukci-bioinformatics/MGA) is another similar tool.
  
  There are then lots of QC tools specific to different types of data, most of which run after alignment. For example [RSeQC](http://rseqc.sourceforge.net/) (RNA data), [Qualimap](http://qualimap.bioinfo.cipf.es/) and many many others. Without specifying what kind of data you have this is a bit difficult to make recommendations for though.
  
  Phil
  
  
  ---
  
  ### Answer by username_6 on 2017-06-04. Score: 5
  
  The quality control of ngs reads is heavily dependent on type of the project.
  
  **For genome assembly** projects based on short reads, beside already covered checking quality of sequencing, you would like to look at the kmer spectra to find out, if your reads are going to make sense when they will be translated to De Brujin graph. You will get also a clues about genome coverage, genome size, repetitive content and small clue about heterozygocity. A lot of useful info about interpretation you can find in the README of GenomeScope.
  
  A list tools I used:
  
  * [jellyfish](http://www.genome.umd.edu/jellyfish.html) - to count k-mer frequencies
  * [GenomeScope](https://github.com/schatzlab/genomescope) - a package for analysis and visualisation of k-mer frequencies (they recommend to use jellyfish for counting k-mer frequencies)
  * [kmergenie](http://kmergenie.bx.psu.edu/) - for prediction of optimal kmer for assembly
  
  Using these tools can save you a lot of frustration if you accidentally sequence a contaminated sample or if your colleagues / a sequencing company have sent you a wrong file!
  
  
  ---
  
  ### Answer by username_8 on 2017-06-06. Score: 2
  
  The kind of QC you do routinely depends on what your lab's focus is. We do a lot of low-quality, multiplexed DNA and RNA. If you routinely do fresh frozen whole genomes, your QC will be different.
  
  Weighing in from the resequencing side of things (i.e. sequencing an organism that has a good reference), there are several things we are testing with quality control:
  
  * Did the library preparation succeed? Was the capture successful? Is our input material of high enough quality?
  * Did the sequencing work? For Illumina, did enough clusters pass filter? What's our Q30? Did we sequence into adapters?
  * Can we proceed with analysis? Are the indices on multiplexed samples correct?
    Do we have the correct reference sequence? Do we have the right targets?
  
  You may use a combination of tools:
  
  * [bcl2fastq](https://support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.html) : if you multiplexed, check the undetermined indices file after converting to fastq to make sure you caught all of the indexed reads
  * [FastQC](https://www.bioinformatics.babraham.ac.uk/projects/fastqc/): General sequencing quality, sequencing into adapters, handy flagging for obvious problems
  * [Illumina interop](https://github.com/Illumina/interop)/SAV: Cluster pass filter, Q30
  * Alignment (with [BWA](http://bio-bwa.sourceforge.net/), bowtie, etc), [bedtools](http://bedtools.readthedocs.io/en/latest/), [samtools](http://www.htslib.org/), [Picard](https://broadinstitute.github.io/picard/): check for contamination, make sure we're using the right reference, overlap with target regions, depth of coverage, soft-clipping for low-quality sequence
  
  
  ---
  
  ### Answer by username_9 on 2017-06-20. Score: 0
  
  If you are looking for NGS QC for your fastq, bam, bed and vcf files I would suggest a commercial tool called [omnomicsQ](https://om.euformatics.com/apex/f?p=118:1:).
  
  It automates coverage analysis and sample pass, warn fail according according to your defined SOPs. Unlike fastqc it also comes with a database so that you can compare protocols, samples, and your performance to peer organisations. It also charts performance over time, exceptions,and correlations between metrics
  
  
  ---
  
  