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!

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.
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.

#### ☞ 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.

Picard 2.10.4 has MAJOR CHANGES that impact throughput of pipelines. Default compression is now 1 instead of 5, and Picard now handles compressed data with the Intel Deflator/Inflator instead of JDK.
GATK version 4.beta.2 (i.e. the second beta release) is out. See the GATK4 BETA page for download and details.

# HaplotypeCaller 2.4

Member

I am getting the following error. What is the minimum read size to do assembly? 50 basepair too short?

##### ERROR stack trace

java.lang.IllegalStateException: Reads are too small for use in assembly.

##### ERROR ---------------------------------------------------------------------------------------

::::::::::::::

Tagged:

• Member

Actually the reads on this sample is 200 bases. It is from a MiSeq run and works fine with version 2.3.

Ooh, first bug report for 2.4 already! That was fast.

We'll take a look at it. Could you upload a snippet of your bam where the error occurs? We'll probably need to reproduce this locally to figure out what's going on here.

• Member

Let me check a few other things and let this run finish. I will be interested to see if the interval that I am running this on returns any snps/indels or if it exits out.

Speed has improved through by a factor of 4?

Sure, no problem. Ryan tells me it's probably a region where lots of reads are getting clipped and resulting in very short leftovers. Before 2.4 the HC would simply have skipped through the region silently, but now it's flipping out unnecessarily. I expect if you look at the interval in question that may confirm his supposition.

• Member

It would be nice if the stack trace would give me some idea as to which interval it is flipping out about.

• Member

Any Idea what this one means?

java.lang.ArrayIndexOutOfBoundsException: -1

##### ERROR ------------------------------------------------------------------------------------------
• Member

Off-by-one error while taking down the variant context info. Sounds like another edge case that's not being handled properly.

To see which interval HC is choking on, use --debug -- it will print out a lot of info including the interval it's processing.

See the instructions on how to file a bug report here:

• Member

For off-by-one error, the following seems to bypass the issue.
In GenotypingEngine.java:675, replace

 final byte refByte = ref[refPos-1]; if( BaseUtils.isRegularBase(refByte) ) { insertionAlleles.add( Allele.create(refByte, true) ); } 

with:

 // Only create a ref allele if it exists and is a regular base if (refPos > 0) { final byte refByte = ref[refPos-1]; if( BaseUtils.isRegularBase(refByte) ) { insertionAlleles.add( Allele.create(refByte, true) ); } } `

• Member

Is this file in the jar zip file? GenotypingEngine.java
Also I was using the Haplotype caller is it still this file?

This issue has been fixed in the latest version of GATK (2.4-7).

• Member

hrmm. I am using
(version 2.4-7-g5e89f01):

##### ERROR

I will work on getting an example for you guys. I did not have time in my last attempt unless you tell me this is an old version of 2.4-7

• Member

I have uploaded the snippet of the problem with associated scripts that I was using.
-rw-r--r-- 1 depristo wga 287606741 Mar 15 17:37 upload-minusoneError.zip

• Member

Geraldine,
Did you see I posted the snippet for the error above?
Thanks