The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Picard 2.9.0 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

Mutect2 dont list Variant

EADGEADG KielMember Posts: 101 ✭✭

Hi,
Iam running mutect2 on capture data, with tumor normal. The data flows the "classic" way on a local wdl-pipeline with bwa =>merge=>markDup=>recal=>mutect2 (3.6).

I see the variant clearly in IGV in both tumor and normal sample and even in the bamout i created for it (see screenshot). But it dont get listed in the output vcf.

I know that mutect2 apply some hard filters to the variant, but I should see the variant in the output marked with an alt_allele_in_normal ?

Greetings

EADG

image
=>tumot
=>normal
=>bamout

Screenshot from 2017-02-13 14:27:03.png
607 x 1028 - 19K
Tagged:

Best Answer

Answers

  • shleeshlee CambridgeMember, Broadie, Moderator Posts: 524 admin

    Hi @EADG,

    I see the variant clearly in IGV in both tumor and normal sample

    MuTect2 is a somatic caller that only calls on differences in the tumor from the normal. However, as you say, the call should be present with an annotation that filters it. Your screenshots, are these the raw BAMs or the BAMOUTs from a MuTect2 run?

  • EADGEADG KielMember Posts: 101 ✭✭

    Hi @shlee,
    thank you for you answer...the screenshots are showing in order tumor.raw.bam, normal.raw.bam and tumor-normal.bamout.bam.

  • shleeshlee CambridgeMember, Broadie, Moderator Posts: 524 admin

    @EADG,

    The BAMOUT needs to by grouped by sample, colored by HC tag and sorted by HC tag. Also, can you show the reference sequence and genomic location for this locus? You've previously done a helpful IGV coloring at the end of this thread for a MuTect2 BAMOUT.

    In general, we need more information than what you've so far provided to effectively answer questions such as why a variant wasn't called. Please help us help you more effectively by provided the needed information up front so we can minimize this back and forth. We are a small team. Thanks.

  • EADGEADG KielMember Posts: 101 ✭✭

    @shlee
    sry:(...i made a new screenshot and hope it will satisfy and also I will provide some helpful information (hope so :).
    I run the pipline in the broadinstitute/genomes-in-the-cloud 2.2.4-1469632282 Docker container with cromwell version 21, so its gatk 3.6.
    Furthermore i use scatter/ghatter and merging the resulting vcf´s. My command for mutect2 is:

        java -Xmx4g -jar /usr/gitc/GATK36.jar \
        -T MuTect2 \
        --dbsnp ${dbSNP_vcf} \
        -dt NONE \
        -R ${ref_fasta} \
        -I:normal ${input_normal_bam} \
        -I:tumor ${input_tumor_bam} \
        -L ${interval_list} \
        -o ${vcf_basename}.MT.vcf
    

    Say a what you want and I deliver more information :).

    Thank you for your help!

    image

    bamout-2.JPG
    1361 x 596 - 136K

    Issue · Github
    by shlee

    Issue Number
    1748
    State
    closed
    Last Updated
    Assignee
    Array
    Milestone
    Array
    Closed By
    sooheelee
  • shleeshlee CambridgeMember, Broadie, Moderator Posts: 524 admin

    Thanks for supplying a better screenshot. I hope there is an easy explanation for this odd phenomena @EADG. Would you mind generating a snippet of the data needed to recreate this image, i.e. your tumor and normal raw BAMs as well as your BAMOUT and callset? Instructions for sending us data are at http://gatkforums.broadinstitute.org/gatk/discussion/1894/how-do-i-submit-a-detailed-bug-report.

  • EADGEADG KielMember Posts: 101 ✭✭

    Hi @shlee

    i manage to upload the Files to your server (gosh that 20 user limitation is quite an annoyance). The Folder is EADG-Mutect2.

    I have observe the same issue with an other sample which carries the same variant...Do you mean i it is useful when I check if i get the variant called with HC, since the relationship between HC and Mutect2 ?

    Thank you!

    Greetings EADG

  • EADGEADG KielMember Posts: 101 ✭✭

    Hi @shlee,

    just to be shure I checked the position an noticed that all of them have rs-numbers and are included in dbSNP. I know you have plenty off other stuff to do...argh but you have a little status update for me ?

    Greetings EADG

  • shleeshlee CambridgeMember, Broadie, Moderator Posts: 524 admin

    Hi @EADG,

    @Sheila has returned from her conference and she should be able to help you going forward, especially given your question requires in-depth probing. I've just been filling in for the team this last week since most are away.

  • SheilaSheila Broad InstituteMember, Broadie, Moderator Posts: 4,695 admin

    @EADG
    Hi,

    I will have a look and get back to you soon.

    -Sheila

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,647 admin
    As I recall M2 doesn't call obvious germline variants, or at least it's not supposed to, no? My impression was that we only output "germline_risk" records if they have a chance of being somatic and are borderline. If they're obviously germline, why bother? Would just make much bigger files for nothing... But I may have misunderstood the default behavior.

    Geraldine Van der Auwera, PhD

  • EADGEADG KielMember Posts: 101 ✭✭

    Hi everbody and thank for your effort,
    if its like @Geraldine_VdAuwera said it is ok, but I need an "official" confirm, because I find nothing about it in the doc´s.
    Maybe a verbose-mode were good were mutect2 puts every think out even it is a obviously germline variant. Just to be shure everything is allright.

    Thanks so far!

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,647 admin

    @EADG We're checking on this with the M2 devs to get you that official confirmation, and document this more clearly as well. Sorry it's taking so long -- part of our team just got back from an international workshop so we're playing catch-up.

    Geraldine Van der Auwera, PhD

  • EADGEADG KielMember Posts: 101 ✭✭

    @Geraldine_VdAuwera

    Good to know :) Thanks for your efforts! Hope Belgium was nice!

Sign In or Register to comment.