Question: Empty VCF files produced by both Mutect2 and HaplotypeCaller?
gravatar for gprashant17
15 months ago by
gprashant1770 wrote:

I have tried using both Mutect2 and HaplotypeCaller (GATK) for calling somatic variants and generating VCF file from an aligned BAM file (RNA-seq). However, the generated VCF file did not have any information other than the header. I have seen this type of thread online, but I did not manage to find a solution. Is this because Base Recalibration was not done on the BAM file?

Meanwhile, when I try to use samtools mpileup on the same BAM file, the process aborts after a long time and it generates a VCF file of around 45GB.

ADD COMMENTlink modified 14 months ago by Biostar ♦♦ 20 • written 15 months ago by gprashant1770

Can you samtools view -h you bam file?

ADD REPLYlink written 15 months ago by Asaf8.4k

Yes, I am able to view the bam files using samtools. I should mention here that I have used AddOrReplaceReadGroups (GATK) to include the @RG tag in the header, as it was not present initially.

ADD REPLYlink written 15 months ago by gprashant1770

It would help if you could show us all the output (STDOUT and STDERR) from the HaplotypeCaller run.

ADD REPLYlink written 15 months ago by RamRS30k

I'm not sure which output should I post. The command ran with no errors. It showed ProgressMeter for all the traversals and ended like this:

17:39:21.868 INFO  ProgressMeter - Traversal complete. Processed 10840980 total regions in 19.1 minutes.
17:39:21.915 INFO  VectorLoglessPairHMM - Time spent in setup for JNI call : 0.0
17:39:21.915 INFO  PairHMM - Total compute time in PairHMM computeLogLikelihoods() : 0.0
17:39:21.916 INFO  SmithWatermanAligner - Total compute time in java Smith-Waterman : 0.00 sec
17:39:21.916 INFO  Mutect2 - Shutting down engine
[21 June, 2019 5:39:21 PM IST] done. Elapsed time: 19.87 minutes.
Tool returned:
ADD REPLYlink modified 15 months ago by RamRS30k • written 15 months ago by gprashant1770

Please give us the exact command and upload the entire log file to a GitHub Gist so we can work together and figure out what's happening.

Guide to GitHub Gist: A: How to Use Biostars Part-3: Formatting Text and Using GitHub Gists

Also, if you're pasting text here, please use the formatting bar (especially the code option) to present it better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.

ADD REPLYlink written 15 months ago by RamRS30k

This might the same problem/solution: Mutect2 (GATK4, Java8) outputs empty VCF file for tumor only mode

ADD REPLYlink written 14 months ago by igor11k
Please log in to add an answer.


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