Mutect2 resources guide

Posted by shlee on 31 Jan 2018 (1)


A new tutorial for somatic calling

We have a new tutorial, Tutorial#11136, that outlines how to call somatic short variants, i.e. SNVs and indels, with GATK4 Mutect2. The tutorial provides small example data to follow along with.

Mutect2-compatible germline resources

Full-length Mutect2-compatible human germline resources are available on our [FTP server]( https://software.broadinstitute.org/gatk/download/bundle) and at gs://gatk-best-practices/. The resources are simplified from the gnomAD resource and retain population allele frequencies. Mutect2 and GetPileupSummaries are the two tools in the workflow that each require a germline resource.

Working WDL scripts

If you want to run the Somatic Short Variant Discovery Best Practices workflow using WDL, be sure to check out the official Mutect2 WDL script in the gatk-workflows repository. @bshifaw and other engineers optimize the scripts in the repository to run efficiently in the cloud. Furthermore, the scripts come with example JSON format inputs files filled out with publically-accessible cloud data.

For other Mutect2-related scripts, e.g. towards panel of normals generation, check out the gatk repository's scripts/mutect2_wdl directory. Our developers update these scripts on a continual basis.

For background information

If you are new to somatic calling, be sure to read Article#11127. It gives an overview of what traditional somatic calling entails. For one, somatic calling is NOT just a difference between two callsets in that germline variant sites are excluded from consideration.

For those switching from GATK3 MuTect2, Blog#10911 will bring you up to speed on the differences.

An off-label tutorial for simple difference calling

If you are interested in simply calling differences between two samples, Blog#11315 outlines an off-label two-pass Mutect2 workflow. Off-label means the workflow is not a part of the Best Practices and is therefore unsupported. However, if given enough community interest, we may be convinced to further flesh out the workflow. Please do post to the forum to express interest.



Return to top

Wed 31 Jan 2018



- Recent posts


- Upcoming events

See Events calendar for full list and dates


- Recent events

See Events calendar for full list and dates



- Follow us on Twitter

GATK Dev Team

@gatk_dev

@nparmalee @raymondpirouz Hah no our pipelines are rather less controversial! Be sure to check out our pre-built WD… https://t.co/nAV3JvTiE5
11 Nov 18
@RubenVanPaemel @bruler0111 @strnr M2 definitely better and has a tumor-only workflow now, which we would expect to… https://t.co/SiCMMKg53z
8 Nov 18
Our very own Andrey Smirnov is delighted to answer questions about the #GATK4 pipeline for common and rare germline… https://t.co/cpCe37W4Ul
17 Oct 18
@cpj3131 Depending on the tool you're using, some filters get applied that may modify what ends up getting counted.… https://t.co/yTGPKrsB9A
15 Oct 18
@cpj3131 It's a variant context annotation that counts the number of A, C, G, and T bases across all samples, in th… https://t.co/gVhAdpqOIv
13 Oct 18

- Our favorite tweets from others

Amazing talk by @dgmacarthur about the expansion of gnomAD and how size and diversity increase filtering power https://t.co/vQAULIMwlw
30 Oct 18
@geoffjentry Who doesn't love a Warp Pig? @WDL_dev and @gatk_dev are on the ball getting stickers out. Was happy to… https://t.co/91OODRpFOC
22 Oct 18
#ASHG18 VA: call with GATK @gatk_dev. Look for pathogenic / likely pathogenic. leverage ClinVar.
17 Oct 18
If you think your fascination with #GATK hit the roof wait until you meet @gatk_dev team! Has been a wonderful week… https://t.co/KwfHm1SzNh
21 Sep 18
@xdopazo @gatk_dev @ClinicalBioinfo @FProgresoysalud @INB_Official @CIBERER @jpflorido Thank you very much for such amazing time!
21 Sep 18

See more of our favorite tweets...