We are trying to identify somatic mutations by comparing DNAs of normal and disease tissues. Because of the small size of the disease tissues, the DNA concentration or amount for some samples are very small. We sent 10 samples to a company to try, and 2 out of 10 (each sample generates 2 datasets) data sets failed in some of quality check. We used FastQC for QC. For these samples, we asked average coverage of 60.
The qualities that failed are:
- Per base sequence quality
- Per base GC content
- Per base sequence content
Should we increase the coverage requirement to 90 to have better quality? In general, how much coverage would give decent quality of data? We use Illumina HiSeq platform.
•
written
by
haiying.kong
2 answers
First, average coverage is not a really useful number. You'll want to think in terms of the proportion of bases (in the exome) covered at X (where X is the target coverage). For example, you might aim for 85% of bases covered at 100x or higher coverage.
Now, for somatic sequencing, you'll probably want to aim for significantly higher for the tumor than for the normal sample. Folks routinely sequence to 150X (per base target) and higher to allow one to discover somatic variants in samples with normal tissue admixture, lower variant allele frequency. What you choose will depend on your biological questions, your budget, and your sample purity.
•
written
by
Sean Davis
Log in to answer this question.
-
Per base sequence content of Metagenomic samples
written by Savindu •
-
why my fastqc on trimmed data failed?
written by kamran.kouchakzadeh.kk •
-
Read length for single -ended RNA seq data analysis
written by vinaya •
-
Biased GC content distribution in miRNA Illumina sequencing data
written by madhuratathode •
-
how to remove warning and failures observed inQuality report by Fastqc for RNA-Seq data
written by aiswaryabioinfo •
-
Do I have to check every fastqc to make good quality?
written by khaeuk
-
FastQC Results -- Per Sequence GC Content and Kmer Contents
written by haiying.kong
-
How many samples would be enough to estimate the quality of sequencing
written by haiying.kong
-
Trimmomatics parameters for WGS data quality improvement
written by ravi.uhdnis
-
Ngs-Qc-Per Base Sequence Content Failed
written by williamjohn360
Your figures / links are not working, the correct way to show images on biostars is posting them on a site like imgur and posting the link here.
or try postimage.org to post your figures here