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

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.8k views
ADD COMMENT
2
Entering edit mode
7.9 years ago
brentp 23k

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.

ADD COMMENT
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....

ADD REPLY
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
ADD REPLY
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.

ADD REPLY
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

ADD REPLY

Login before adding your answer.

Traffic: 2999 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6