Quality Control Software For Targeted Sequencing?
1
3
Entering edit mode
8.4 years ago
Luca Beltrame ▴ 240

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?

ngs analysis sequencing software • 2.9k views
2
Entering edit mode
8.4 years ago
brentp 24k

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.

0
Entering edit mode

Looks interesting. Now all that's "left" is to find a way to convert Illumina's manifest format to Picard interval....

0
Entering edit mode

to create that file, just do:

samtools view -H your.bam > intervals.list
awk 'BEGIN{FS=OFS="\t"}{ print $0"\t+\t"$1":"$2"-"$3 }' illumina.bed >> intervals.list

0
Entering edit mode

Also in the latest versions of Picard there is a command that helps one to make such intervals from Bed files. They are: picard BedToIntervalList. This command requires a sequence dictionary that can be made using: picard CreateSequenceDictionary.

0
Entering edit mode

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