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 ??
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 ??
Comment