For analyzing QC of fastq files I know fastQC; It generates lots of statistics; Also, there are a set of tools (GATK, samtools, bedtools) to check the quality of reads after mapping.(bam files)
My question is, given different exome sequencing batches, and say different samples per batch; what parameters/figures would you make to present the quality of data in different steps ?
•
written
by
User 1933 •
1 answer
There is no single tool or even data level that you can track that will catch all problems, but a good high-level metric is the % of bases with coverage > 30x (30x is a minimum) for targeted bases. This number will get you a rough sense of how effective your variant calling will be. You can get this information from picard CalculateHsMetrics.
Others that are less important for gauging the success on a per-sample basis and can be tracked over time to assess systematic sequencing problems might include:
- Insert sizes
- Error rate
- Duplication rate/effective library size
- % aligned
- Overrepresented sequences
- Adapter contamination
Log in to answer this question.
-
RRBS fastq - biased per base sequence content
written by mbk0asis
-
Finding Overlapping Variants (I.E. Indels, Snps) Using Annovar Format.
written by youngcsong •
-
Non-Coding Region On Hg19
written by User 1933 •
-
How To Count The Reads Mapped To Exon And Exon-Exon Junctions?
written by lwc628
-
How To Fix A Bam-Readcount Sm Error: "Couldn'T Grab Single-End Mapping Quality For Read"
written by amnonim •
-
Converting Gff To Bed With Bedtools?
written by user
-
Reducing I/O Burden Of Bamtools C++ Api
written by Mitch Bekritsky
-
Dssp And 6 Different Hydrogen Bonds - What Are They Exactly
written by User 1933 •
-
Softwares For Calculation Of Genome-Wide Alignments Coverage
written by Bioscientist
-
Getting Unmapped Reads: Comparing Fastq To Bam
written by User 9996
check this post: fastQC html report to PDF (with a script)
the post you have refereed to, is about converting html to pdf; I don't want to make 300 pdf from 300 samples, ...
in my ans in that post i selected few images out of FASTQC ; those will be useful to you...
I'm a little unclear on what you are asking. Is your question "Of all the data that I can get out of a quality control assessment, what would be the most important criteria to communicate in regards to my quality control pipeline/algorithm?"