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!

Get notifications!


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.
2. Try the latest version of tools.
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.

Did we ask for a bug report?


Then follow instructions in Article#1894.

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.

Jump to another community
Download the latest Picard release at https://github.com/broadinstitute/picard/releases.
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

indapaindapa Member
edited August 2012 in Ask the GATK team

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

Answers

  • 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!

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    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.

  • marzimarzi 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?

  • marzimarzi 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.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    @marzi, it sounds like you may be using the wrong reference genome.

Sign In or Register to comment.