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.

differences between mutect1 and mutect2 output

tommy1984tommy1984 RomeMember

Hi I am trying to use mutect2 on my data.
Comparing mutect1.1.7 output and mutect2 output I noticed that the somatic calls are totally different.

Here is reported the mutect-1.1.7 command line:

java -jar /data01/Software/mutect/mutect-1.1.7.jar --analysis_type MuTect \
--reference_sequence human_hg19.fa \
--cosmic cosmic_chr.vcf \
--dbsnp dbSNP_144_00-All_chr.vcf \
--intervals SureSelect_All_Exon_V4_chr.list \
--input_file:tumor tumor.bam \
--input_file:normal normal.bam \
--out classical.out \
--coverage_file out.txt \
--vcf mutect1.vcf &

The input bam files are preprocessed using the GATK BaseRecalibrator and the GATK IndelRealigner.

Here is reported the mutect2 command line:

java -jar GenomeAnalysisTK_3.5.jar -T MuTect2 -R human_hg19.fa \
-I:tumor tumor.bam \
-I:normal normal.bam \
--dbsnp dbSNP_144_00-All_chr.vcf \
--cosmic cosmic_chr.vcf \
-L SureSelect_All_Exon_V4_chr.list \
-o mutect2.vcf

In this case the input bam files are preprocessed without the use of GATK IndelRealigner.

Comparing mutect1.vcf and mutect2.vcf there is only one common somatic variant call.


  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    We do expect to see differences in the results between these two versions. MuTect2 has some dramatic operational changes compared to the original MuTect. Have you looked at their respective calls to evaluate whether the results look better in one case or the other? Can you describe the differences in more detail?

Sign In or Register to comment.