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.

#### ☞ Did you remember to?

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.

GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

# ReduceReads problem in GATK v2.6-2-ge03a5e9, Compiled 2013/06/20 20:07

Oxford, UKPosts: 9

Hello,
While trying to run ReduceReads on my BAM outputs of PrintReads, using the latest build v2.6-2-ge03a5e9, I get the following error, kindly comment on this error trace:

Command:
java -Xmx8g -jar /path/to/my/server/ashish/tools/GenomeAnalysisTK-2.6-2-ge03a5e9/GenomeAnalysisTK.jar -T ReduceReads -R /path/to/my/server/pipeline/lib/hs37d5.fa -I /path/to/my/server/output/bamext.base.recal/41453_TTAGGC_L002_001.bam -o /path/to/my/server/output/bamext.reduced/41453_TTAGGC_L002_001.bam

##### ERROR stack trace

org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Removed too many insertions, header is now negative at position 151039408
at org.broadinstitute.sting.gatk.traversals.TraverseReadsNano$TraverseReadsReduce.apply(TraverseReadsNano.java:251) at org.broadinstitute.sting.gatk.traversals.TraverseReadsNano$TraverseReadsReduce.apply(TraverseReadsNano.java:240)

##### ERROR ------------------------------------------------------------------------------------------
Tagged:

• Oxford, UKPosts: 9
edited June 2013

I've uploaded the bam snippet named "ReduceRead_problem_data.tar.gz" and reproduced the error on this snippet locally.
Thanks for the prompt attention.

Thanks for the test data! We'll let you know when we have an update in this thread.

Geraldine Van der Auwera, PhD

• Posts: 8

Hi,I met the same problem.Would you tell me the way to resolve?

• Oxford, UKPosts: 9

I'm also waiting for an updated solution from the GATK team.
Hopefully, they should be able to replicate the error and then resolve this bug from my uploaded bam snippet and thereupon, release a patched up version of this version of GATK release.

Sorry about the updated "answer tag" with this question. I'd posted for the first time and didn't realise that opting yes for "Did this answer the question" would update the final status flag of the main question rather than for that specific assistance step.

Hi folks,

We have processed the test files (thanks Ashish!) and were able to replicate the error. However we haven't yet had time to solve it. We are a little short-staffed right now with several members on vacation, so it will take a bit longer than usual to get bugs fixed. But we'll do our best to minimize the wait time. Thanks for your patience!

(I updated the answer flag to mark this as unresolved)

Geraldine Van der Auwera, PhD

• Oxford, UKPosts: 9

Thanks Geraldine.
However, on the latest release notes, I see this mentioned ##Added fix for the "Removed too many insertions, header is now negative" bug.##

That case was addressed and resolved; we think that your case is a different issue causing the same symptom.

Geraldine Van der Auwera, PhD

• Oxford, UKPosts: 9

I see. Now one other twist to the scenario - with the GATK verion 2.6, I've been able to generate the BAM Print Reads for all BAMs (n~1100+), but the Reduce Reads affect some BAMs, while other BAMs do generate upon the reduction step.

However, the same error is not replicated when using the previous GATK version 2.5 on that BAM snippet for Reduce Reads. So what if I run the last step of BAM reduction using the GATK version 2.5 (assuming that it runs successfully on all BAMs) on the Print Reads from the GATK version 2.6? Would that be a reasonable hit to take, in interest of time to be able to run Unified Genotyper on this large no. of BAMs? Or you reckon there would be other ramifications on the quality of calls generated?

• Oxford, UKPosts: 9

That's wonderful, thanks.

• Posts: 6

Hi there, I have the same problem with ReduceReads.
Is the fix already available ?
Thanks.

• Oxford, UKPosts: 9

Hi,

I can confirm that the nightly build of 26/06/2013 did not give me an error on the BAM snippet that I'd provided to the GATK team.
I am now running the reduction step on all my BAMs with this new build and shall get back if I see any further issues.

GenomeAnalysisTK-nightly-2013-06-26-g9fa3206

• Oxford, UKPosts: 9

I've finished the reduction step for all 1100+ BAMs and encountered no problems in the runs.
I also noticed faster run times than previously seen.
Thanks GATK team :-)

• Posts: 8

Hi.I met the same error:

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

Would you tell me the way to resolve?

Thanks!

Hi Hong,

This bug is fixed in our development version, and the fix is available in the latest nightly build (see the Downloads page).

Geraldine Van der Auwera, PhD

• Posts: 1

Is there a nightly build of GATK Queue available?

No, development on Queue is not so fast as to need nightly builds.

Geraldine Van der Auwera, PhD

• Posts: 16

Hi I am now getting the same error message using the latest GATK:

##### ERROR ------------------------------------------------------------------------------------------
• Oxford, UKPosts: 9

Try the version in the nightly build.

• Posts: 16

problem solved with the nightly build.

Thank you