SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
GATK "Code exception" memento Bioinformatics 2 02-07-2012 12:41 AM
GATK IndelRealigner error kga1978 Bioinformatics 7 12-16-2011 08:58 AM
GATK IndelRealigner error Heisman Bioinformatics 2 09-14-2011 03:25 PM
GATK IndelRealigner returns empty output file bjfan Bioinformatics 3 06-08-2011 10:43 AM
GATK IndelRealigner error carljo Bioinformatics 1 04-06-2011 11:25 AM

Reply
 
Thread Tools
Old 12-14-2010, 03:25 AM   #1
HongxiangZheng
Junior Member
 
Location: China

Join Date: Dec 2010
Posts: 5
Default GATK IndelRealigner Error:Code Exception

Hi, all

Sorry to interrupt again.

I did some work on Local realignment around indels and followed the steps online.

The step 1 was ok.When I ran the 2nd command,An error appeared.
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR stack trace
java.lang.ArrayIndexOutOfBoundsException
at java.lang.System.arraycopy(Native Method)
at org.broadinstitute.sting.utils.sam.AlignmentUtils.createIndelString(AlignmentUtils.java:669)
at org.broadinstitute.sting.utils.sam.AlignmentUtils.leftAlignIndel(AlignmentUtils.java:568)
at org.broadinstitute.sting.gatk.walkers.indels.IndelRealigner.determineReadsThatNeedCleaning(IndelRealigner.java:807)
at org.broadinstitute.sting.gatk.walkers.indels.IndelRealigner.clean(IndelRealigner.java:591)
at org.broadinstitute.sting.gatk.walkers.indels.IndelRealigner.cleanAndCallMap(IndelRealigner.java:463)
at org.broadinstitute.sting.gatk.walkers.indels.IndelRealigner.map(IndelRealigner.java:437)
at org.broadinstitute.sting.gatk.walkers.indels.IndelRealigner.map(IndelRealigner.java:68)
at org.broadinstitute.sting.gatk.traversals.TraverseReads.traverse(TraverseReads.java:101)
at org.broadinstitute.sting.gatk.traversals.TraverseReads.traverse(TraverseReads.java:49)
at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:67)
at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:103)
at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:97)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:244)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:87)
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A GATK RUNTIME ERROR has occurred (version 1.0.4705):
##### ERROR
##### ERROR Please visit to wiki to see if this is a known problem
##### ERROR If not, please post the error, with stack trace, to the GATK forum
##### ERROR Visit our wiki for extensive documentation http://www.broadinstitute.org/gsa/wiki
##### ERROR Visit our forum to view answers to commonly asked questions http://getsatisfaction.com/gsa
##### ERROR
##### ERROR MESSAGE: Code exception (see stack trace for error itself)
##### ERROR ------------------------------------------------------------------------------------------

The input file should be OK, for I ran validsamfile and only found some information" CIGAR should have zero elements for unmapped read", which was considered no problem in Pichard website. The input file is sorted.

Can anyone help?


Additional, I thought the problem might be:

1.those unmapped read have non zero elements in Cigar.
or

2 there might be invalid charactor such as "W" or "K" in reads.

Then, How can I remove those invalid reads or charactors?
HongxiangZheng is offline   Reply With Quote
Old 12-14-2010, 07:21 PM   #2
HongxiangZheng
Junior Member
 
Location: China

Join Date: Dec 2010
Posts: 5
Default

It seems solved.For those reads mapped off the end of the contig might cause the problem.

New version of GATK would solved it.
Thank Erik Banks @http://getsatisfaction.com/gsa/.
HongxiangZheng is offline   Reply With Quote
Reply

Tags
gatk, indels

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:48 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO