If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra

Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
We will be out of the office on October 14, 2019, due to the U.S. holiday. We will return to monitoring the forum on October 15.

Allelic Coverage Pulldown at genomic loci

ncamardancamarda DFCIMember, Broadie

Is there a GATK tool that can pull down the allelic coverage at a site in the genome for all alleles? The output format I'm imaging is (A,C,T,G) counts for each genomic loci I specify. Can someone please advise?
Thank you so much.

Best Answer


  • ncamardancamarda DFCIMember, Broadie

    Hi Sheila, I've been using the DepthOfCoverage tool from GATK 3.8.0 and this has been satisfactory for my purposes. Because none of the somatic variant calling tools in the GATK perform joint calling on tumor samples, I've been trying to implement a force-caller, and for this, I need to be able to identify whether - across samples - a given base has been covered. I've found that it's not necessarily true that the alleles that M2 emits for a given site encompass all of the bases that have coverage at that location.

    Do you think the GATK will ever implement such joint calling for MuTect2?

  • SheilaSheila admin Broad InstituteMember, Broadie, Moderator admin


    I remember asking the team about this, and the response was no. I cannot find the exact thread that discusses this, but perhaps this thread will help.


Sign In or Register to comment.