Question: Can anyone help me with normal-only calling for panel of normals(PON) with MuTect2.
0
gravatar for ShahiRB
3.3 years ago by
ShahiRB0
Brussel
ShahiRB0 wrote:

Can anyone help me with normal-only calling for panel of normals with MuTect2.

The code below does not spit any error as well as variants in vcf, except vcf header.

Here is the code:

java -Xmx45g -XX:+AggressiveOpts -jar /nexus/GenomeAnalysisTK-3.5/GenomeAnalysisTK.jar \
-T MuTect2 \
-R /home/ngs/data/ucsc.hg19.fasta \
-I:tumor /nexus/normal.bam \
--dbsnp /home/ngs/data/dbsnp_137.hg19.vcf \
--cosmic /nexus/rshahi/data/CosmicCodingMuts.vcf5 \
--artifact_detection_mode \
--intervals /nexus/NimbV3_target.gatk.interval_list \
-o $HOME/jobout/normal_only.vcf

Thanx,
Rajendra
Brussels

next-gen software error • 4.6k views
ADD COMMENTlink modified 3.0 years ago by 2nelly130 • written 3.3 years ago by ShahiRB0
2
gravatar for Dan Gaston
3.3 years ago by
Dan Gaston7.1k
Canada
Dan Gaston7.1k wrote:

The most likely source of error is your intervals list. If this is incomplete for instance, or not matching up with the regions you have data for, you would get no variants and also no error messages (unless the intervals list file was actually malformed). Try running the same command but without passing it an intervals file and see what happens.

You should also probably check your BAM file to make sure it is actually containing aligned data as expected. Visualize it with IGV for instance and check some genomic coordinates of your favourite genes as a quick check.

ADD COMMENTlink written 3.3 years ago by Dan Gaston7.1k

Hi Dan,

1.Even without passing intervals list I got the same results i.e no error, no variants in vcf, header only.

2.The bam file is perfect.

Any further suggestions?

Regards,

ADD REPLYlink written 3.3 years ago by ShahiRB0

Maybe try changing your java calling? It's possible something is going on with the JVM. You are passing 45G of RAM and Aggressive Opts to Java. If the BAM file is for a single exome or targeted resequencing experiment (assuming since you're passing a Nimblegen targets list) then 45G is about 10x-20x what is actually needed. If GATK runs out of memory it will throw an error, but if the problem is earlier on with the JVM then it could potentially die silently, especially with the -XX:+AggressiveOpts also being passed. That will turn on a whole suite of options and I'm not sure what effect that may have as I have never used it.

If neither of those seem to be the case it is hard to diagnose without trying on your same BAM file (or a portion of the BAM file) to diagnose. There doesn't seem to be anything wrong with the command line. When you look at the BAM file in IGV or something can you find obvious variants bye eye that should be called? If you throw the BAM at FreeBayes or another variant caller do you get results then?

ADD REPLYlink written 3.3 years ago by Dan Gaston7.1k

@Dan, I tried without XX:+AggressiveOpts argument; still the same problem. But the HaplotypeCaller runs perfect with the same bam and intervals list.

ADD REPLYlink written 3.3 years ago by ShahiRB0

Well I'm stumped. I would suggest you post with the GATK directly. Since MuTect2 is still in Beta you may have come across some sort of bug or edge-case.

ADD REPLYlink written 3.3 years ago by Dan Gaston7.1k

@Dan, Thanx for your effort. I will update you with some tweaks from Geraldine_VdAuwera(GATK). Cheers!

ADD REPLYlink written 3.3 years ago by ShahiRB0
1
gravatar for 2nelly
3.0 years ago by
2nelly130
Athens
2nelly130 wrote:

Don t worry your command is correct, it s not your fault. There is a bug in MuTect2 when try to create a PON

[http://gatkforums.broadinstitute.org/gatk/discussion/6789/mutect2-panel-of-normals-empty-vcf-fileption][1]

I am facing the same issue.

Some people claim that the problem was solved but I am still getting an empty vcf file.

Try the nightly build

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by 2nelly130
0
gravatar for igor
3.3 years ago by
igor7.6k
United States
igor7.6k wrote:

It looks like the --artifact_detection_mode flag is for HaplotypeCaller (the GATK recommendation is to use HaplotypeCaller for generating the panel of normals): https://www.broadinstitute.org/gatk/guide/tooldocs/org_broadinstitute_gatk_tools_walkers_cancer_m2_MuTect2.php

 

ADD COMMENTlink written 3.3 years ago by igor7.6k
1

I believe that is a mistake in the documentation. I don't think HC has the artifact_detection_mode as an option.

ADD REPLYlink written 3.3 years ago by Dan Gaston7.1k

I haven't tried MuTect2 yet. I did think that was odd, which is why I thought it might be causing the error.

ADD REPLYlink written 3.3 years ago by igor7.6k
0
gravatar for ShahiRB
3.3 years ago by
ShahiRB0
Brussel
ShahiRB0 wrote:
Thanx igor, but that doc abt mutect2 looks misleading . Actually HC does not have that argument.
ADD COMMENTlink written 3.3 years ago by ShahiRB0

Hi biotechnepal, this should really be added as a comment to igor's answer as opposed to an answer itself.

ADD REPLYlink written 3.3 years ago by Dan Gaston7.1k
0
gravatar for winter_li
3.3 years ago by
winter_li40
winter_li40 wrote:

HI,  I think maybe you can find GATK4 to have a try, I heard it has been published.

ADD COMMENTlink written 3.3 years ago by winter_li40

Still in development. 3.5 is still the last official release

ADD REPLYlink written 3.3 years ago by Dan Gaston7.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 925 users visited in the last hour