Notice:
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.

MuTect2 to compare 2 tumor samples

I have capture data of 2 tumor biopsies from 1 patient, and would like to compare the two for mutational similarities and differences. I don't have a normal control.
I tested if MuTect2 (4.0.1.2) could help me in my analysis by comparing the 2, with one of the tumor samples as "normal" and the other as "tumor". Analysis results gave me a nice list of variants that are not present in one tumor but present in the other, however i'm missing the variants that are overlapping (but differ in vaf levels, for example 20% in one tumor and 80% in the other).
Is there a way to tweak MuTect2 in such a way that also the overlapping variants (that are now obviously regarded as 'germline') are outputted?

Tagged:

Answers

  • SheilaSheila Broad InstituteMember, Broadie admin

    @abrink
    Hi,

    The team lead just said he would work on allowing Mutect2 to output germline sites as well as somatic sites. That should be a new feature soon.

    In the meantime, you may consider running Mutect2 in tumor-only mode on both of your samples and using SelectVariants with `--concordance.

    Also, you may find this post helpful. Although it illustrates the opposite of what you want, you may find some helpful tips in there.

    -Sheila

Sign In or Register to comment.