Question: Gatk Threshold Parameter For Clipping Of Reads During Indel Calling Running Unitifed Genotyper?
1
gravatar for 14134125465346445
4.0 years ago by
United Kingdom
141341254653464453.3k wrote:

GATK 1.6 claims to apply a Q20 threshold to clip ends of the reads for the indel caller:

http://www.broadinstitute.org/gatk/guide/article?id=1237
3. Indel Calling with the Unified Genotyper

[...] while many of the parameters are common between indel and SNP calling, some parameters have different meaning or operate differently. For example, --min_base_quality_score has a fixed, well defined operation for SNPs (bases at a particular location with base quality lower than this threshold are ignored). However, indel calling is by definition delocalized and haplotype-based, so this parameter does not make sense. Instead, the indel caller will clip both ends of the reads if their quality is below a certain threshold (Q20), up to the point where there is a base in the read exceeding this threshold.

Also here:
http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_genotyper_UnifiedGenotyper.html

--min_base_quality_score / -mbq ( int with default value 17 )

Minimum base quality required to consider a base for calling. The minimum confidence needed in a given base for it to be used in variant calling. Note that the base quality of a base is capped by the mapping quality so that bases on reads with low mapping quality may get filtered out depending on this value. Note too that this argument is ignored in indel calling. In indel calling, low-quality ends of reads are clipped off (with fixed threshold of Q20).

Can one specify a different threshold other than this "fixed threshold of Q20" when running the Unified Genotyper?

indel gatk • 1.3k views
ADD COMMENTlink modified 4.0 years ago • written 4.0 years ago by 141341254653464453.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 601 users visited in the last hour