The current GATK version is 3.8-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.

Got a problem?


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
Download the latest Picard release at https://github.com/broadinstitute/picard/releases.
GATK version 4.beta.3 (i.e. the third beta release) is out. See the GATK4 beta page for download and details.

HaplotypeCaller throws ERROR: Stack trace on generating raw variants VCF

firasr83firasr83 New YorkMember
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

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    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)?

  • firasr83firasr83 New YorkMember

    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

    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

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    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.

Sign In or Register to comment.