The current GATK version is 3.8-0
Examples: Monday, today, last week, Mar 26, 3/26/04

#### Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

#### ☞ Got a problem?

1. Search using the upper-right search box, e.g. using the error message.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

#### ☞ Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ` ) each to make a code block as demonstrated here.

GATK version 4.beta.3 (i.e. the third beta release) is out. See the GATK4 beta page for download and details.

# Min mapping quality for UnifiedGenotyper

Member
edited August 2012

HI GATK - I am still using the GenomeAnalysisTK-1.6-5-g557da77 version for UnifiedGenotyper. This is probably a silly question, but is there a way to set a parameter for minimum mapping quality score for reads, in deciding whether to evaluate them for variant detection. I know there is a --min_base_quality_score parameter, but I don't see on for mapping quality. http://www.broadinstitute.org/gsa/gatkdocs/release/org_broadinstitute_sting_gatk_walkers_genotyper_UnifiedGenotyper.html

Tagged:

• Member

Hello,

I have this same question. I'm also curious what the minimum mapping quality currently is for Unified Genotyper to consider a read. I see that it will not consider reads with unknown mapping quality, but so far I have not been able to find a minimum value written anywhere. Apologies if I've missed it.

Thanks a lot!

By default by the UG does not look at reads that have unknown or zero mapping quality (using read filters as detailed here, but it does not apply a minimum mapping quality threshold as such. If you want to do that (although we prefer to filter after calling rather than impose pre-calling thresholds) you can use the MappingQuality filter.

• netherlandsMember

if I am right, MQ in VCF file is for filtering mapping quality, But MQ means Root Mean Square of the mapping quality, not mapping quality itself. So, if we use e.g 20 does it mean that I filtered mapping quality less than 20?

• netherlandsMember

Hi all,

I used UG to call my snps. I got so many snps which is not expected. e.g for chr1 of layer chicken which is the largest chr I got 44 million snps. It means almost every positions was called as SNP. I could not find the problem. any idea will be appreciated.