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!

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!


Surround blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block.
Powered by Vanilla. Made with Bootstrap.
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.

HaplotypeCaller throws ERROR: Stack trace on generating raw variants VCF

firasr83firasr83 New YorkMember Posts: 7
edited December 2013 in Ask the GATK team

Hello,

This is in regards to a small exome capture experiment that I'd recently talked about here: http://gatkforums.broadinstitute.org/discussion/3531/bam-index-files-from-1000g-not-recognized-other-questions#latest . There are 8 experimental BAMs and 22 BAMs from 1000G supplied as inputs to HaplotypeCaller. And it quits with a code exception error at 52.1%. I've attached the relevant output. Also attached is the input files list. My OS info is:

Linux 3.2.0-56-generic #86-Ubuntu SMP x86_64 GNU/Linux

Any help would be greatly appreciated!

Thank you,
Firas

txt
txt
gatkdebug1.txt
854K
txt
txt
inputbamlistforrawvcf.list.copy.txt
3K

Answers

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

    Hi Firas,

    This looks like it might be a bug. Can you confirm that all your BAM files are validated (using Picard tools)? Also, could you please check whether this also occurs with the latest version (2.8)?

    Geraldine Van der Auwera, PhD

  • firasr83firasr83 New YorkMember Posts: 7

    Thanks for the reply @Geraldine_VdAuwera. I'll try to report back with what happens with version 2.8. I don't have the original BAM files. Just the ones that have gone through the "reduced reads" step in the pipeline. I had validated the original BAMs before with no issues but trying to validate the post-"reduced reads" BAMs throws up errors, but I guess this is normal. I can send you the verification output if required.

    Stay tuned. Thanks!

  • firasr83firasr83 New YorkMember Posts: 7

    Hi @Geraldine_VdAuwera, so I ran the same input through HaplotypeCaller, GATK version 2.8 and it again threw a NullPointerException error, but this time at 29.5%. Attached is the output.

    Please let me know what you think. Would you suggest that I give up on using HaplotypeCaller and use UnifiedGenotyper instead?

    Thanks,
    Firas

    txt
    txt
    debug.txt
    558K
  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,388 admin

    This may be a problem of incompatibility between versions. If your bams were reduced with an older version of GATK, they should be called with the same version. More recent versions are not backwards compatible with older reduced files because the format evolved rapidly over the last few versions. UG may experience the same compatibility issues.

    I recommend doing whatever is necessary to obtain the original bams. Otherwise you'll need to downgrade to an older version of GATK, but there's no guarantee that will work, and unfortunately we can't provide support for older versions.

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.