To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

Is it possible for GATK4 beta3 Mutect2 to output germline variants also?

By the default configuration, it seems that all variants are somatic. I could not find the parameter which would regulate that.

Tagged:

Best Answer

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    Mutect2 is not a germline variant caller; its modeling algorithm is tuned specifically for somatic mutations. If you're looking for germline variation as well you'll need to use a germline caller like HaplotypeCaller in a separate run.

  • @Geraldine_VdAuwera,
    thanks for the fast reply. In the previous versions, correct me if I'm wrong, the output VCF contained both germline and somatic variants. Was that changed? Or maybe some parameter can prevent outputting of the germline variants. Here is the command line:
    java -Xmx2048M -jar /opt/gatk-package-4.beta.3-local.jar Mutect2 --output Tumor_exomeSeq.vcf --tumorSampleName Tumor_exomeSeq --reference Homo_sapiens.GRCh38.dna.primary_assembly.fa --input tumor_recalibrated.bam --input normal_recalibrated.bam --normalSampleName Normal_exomeSeq --intervals 9.bed --dbsnp dbsnp_144.hg38.vcf

  • SheilaSheila Broad InstituteMember, Broadie, Moderator

    @Vladimir_Kovacevic
    Hi,

    I am not aware of previous versions of Mutect outputting germline variants.

    I am assuming you are asking about sites that are either in the PoN, germline resource or in the matched normal. For the sites in the PoN or germline resource, you can simply take a look at those in the resource files and compare them to the output of Mutect2. For the matched normal, it is best to run a germline variant caller like HaplotypeCaller. Note, FilterMutectCalls will also flag any sites called by Mutect2 that are in the PoN or germline resource or have a high LOD score in the normal.

    -Sheila

Sign In or Register to comment.