View Single Post
Old 12-18-2012, 11:18 PM   #1
ersgupta
Member
 
Location: India

Join Date: Jun 2011
Posts: 26
Default Error while running ContEst

I am trying to run ContEst(http://www.broadinstitute.org/cancer/cga/contest)

Below is the command and options I used:

java -Xmx1200m -jar ContEst.jar \
-T Contamination \
-I test.bam \
-R hg19_GATK_Ref.fa \
-L target_bed_file.bed \
-Bop,vcf hg19_population_stratified_af_hapmap_3.3.with.chr.vcf \
-B:genotypes,vcf test.vcf \
-BTI genotypes \
-o test.vcf.contest.out.txt

Here is the error report:

INFO 17:18:47,638 TraversalEngine - chr15:65621857 1.48e+06 5.2 m 3.5 m 72.4% 7.2 m 119.5 s
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR stack trace
java.lang.ArrayIndexOutOfBoundsException: 0
at org.broadinstitute.cga.tools.gatk.walkers.cancer.contamination.ContaminationWalker.map(Unknown Source)
at org.broadinstitute.cga.tools.gatk.walkers.cancer.contamination.ContaminationWalker.map(Unknown Source)
at org.broadinstitute.sting.gatk.traversals.TraverseLoci.traverse(TraverseLoci.java:78)
at org.broadinstitute.sting.gatk.traversals.TraverseLoci.traverse(TraverseLoci.java:19)
at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:62)
at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:234)
at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:113)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:221)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:87)
##### ERROR ------------------------------------------------------------------------------------------

It seems to be some issue with chr15 because I tried it with 3 different bams files and the reference used to generate those bams is also same(same as the one I used with ContEst).

I am unable to find any such issue reported by anyone else.

Anybody else having this issue or knows a way around ??

Last edited by ersgupta; 12-19-2012 at 07:30 AM.
ersgupta is offline   Reply With Quote