The current GATK version is 3.2-2

#### Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

# HaplotypeCaller 2.4

Posts: 13Member

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

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

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

Tagged:

• Posts: 13Member

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.

Geraldine Van der Auwera, PhD

• Posts: 13Member

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.

Geraldine Van der Auwera, PhD

• Posts: 13Member

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

• Posts: 13Member

Any Idea what this one means?

##### ERROR ------------------------------------------------------------------------------------------
• Posts: 13Member

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:

Geraldine Van der Auwera, PhD

• Posts: 1Member

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) ) {
}

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) ) {
}
}
• Posts: 13Member

Is this file in the jar zip file? GenotypingEngine.java or do I download a different source version of GATK? 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).

Geraldine Van der Auwera, PhD

• Posts: 13Member

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