The current GATK version is 3.6-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!

Powered by Vanilla. Made with Bootstrap.

Min mapping quality for UnifiedGenotyper

indapaindapa Member Posts: 1
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

  • katiefskatiefs Member Posts: 1

    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 Administrator, Dev Posts: 10,704 admin

    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.

    Geraldine Van der Auwera, PhD

  • marzimarzi netherlandsMember Posts: 5

    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 Posts: 5

    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 Administrator, Dev Posts: 10,704 admin

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

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.