Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

GetPileupSummaries all outputs zero

0

4 comments

  • Avatar
    Gökalp Çelik

    Hi Steven Strong

    Your mapping qualities are all 1 for these 2 reads. GetPileupSummaries tool requires a default minimum of 50 as mapping quality for counting reads. 

    --min-mapping-quality,-mmq <Integer>
                                  Minimum read mapping quality  Default value: 50. 

    Can you adjust this value to 1 or 0 and see if it works. Also your aligner may not be an optimal one to provide proper mapping qualities therefore we recommend you review your tools for mapping. 

    I hope this helps. 

    0
    Comment actions Permalink
  • Avatar
    Steven Strong

    That fixed it. Thanks so much! It might be useful if this filtering step was also logged like the other filtering steps, to make it easier to debug issues like this.

    0
    Comment actions Permalink
  • Avatar
    Gökalp Çelik

    Hi Steven Strong

    Thank you for the feedback. A PR is ready to add this log to the default output and will be available in the next release. 

    https://github.com/broadinstitute/gatk/pull/8781 

    Regards. 

    0
    Comment actions Permalink
  • Avatar
    Steven Strong

    thanks!

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk