I'm mainly doing (amplicon-based) resequencing over small (50-100) number of genes with the Illumina technology. Aside from the normal metrics for quality control produced by Qualimap and FastQC, I'm wondering if there's something out there (aside TEQC, see below for while) that targets specifically targetd resequencing experiments.
I'm mostly interested in assessing coverage uniformity across a number of samples sequenced for my targets. I know about TEQC but it has a number of problems, in particular related to reading in paired reads (all of my experiments are paired-end), so I'm looking at alternatives.
Is there any software for QC that can be used for these types of experiments?
•
written
by
Luca Beltrame
picard's CollectTargetedPcrMetrics (http://picard.sourceforge.net/command-line-overview.shtml#CollectTargetedPcrMetrics) is a simple place to start.
It will tell you what % of your target bases are covered at 2X, 10X, 20X, 30X. What percent of your target bases are not covered and a few other useful metrics.
Looks interesting. Now all that's "left" is to find a way to convert Illumina's manifest format to Picard interval....
The link above is broken. Here is a currently unbroken link to CollectTargetedPcrMetrics: https://software.broadinstitute.org/gatk/documentation/tooldocs/4.0.5.0/picard_analysis_directed_CollectTargetedPcrMetrics.php
Log in to answer this question.
