Hi GATK team , FYI I'm getting the following error with GATK 3.7:
<br />java -X -jar GenomeAnalysisTK.jar -T HaplotypeCaller -R human_g1k_v37.fasta --validation_strictness LENIENT -I .bam.list \
-o "out.vcf.gz" -l INFO -nct 10 --dbsnp "138.b37.vcf" --annotation PossibleDeNovo --annotation AS_FisherStrand --annotation AlleleBalance --annotation AlleleBalanceBySample --annotation BaseCountsBySample --annotation GCContent --annotation ClippingRankSumTest --pedigree in.ped -L:BED "in.bed"
##### ERROR stack trace
java.lang.IllegalStateException: Never found start -1 or stop -1 given cigar 100I
at org.broadinstitute.gatk.utils.sam.AlignmentUtils.getBasesCoveringRefInterval(AlignmentUtils.java:204)
at org.broadinstitute.gatk.utils.sam.AlignmentUtils.countBasesAtPileupPosition(AlignmentUtils.java:1418)
at org.broadinstitute.gatk.tools.walkers.annotator.BaseCountsBySample.getBaseCounts(BaseCountsBySample.java:135)
at org.broadinstitute.gatk.tools.walkers.annotator.BaseCountsBySample.annotate(BaseCountsBySample.java:113)
at org.broadinstitute.gatk.tools.walkers.annotator.VariantAnnotatorEngine.annotateGenotypes(VariantAnnotatorEngine.java:517)
at org.broadinstitute.gatk.tools.walkers.annotator.VariantAnnotatorEngine.annotateContextForActiveRegion(VariantAnnotatorEngine.java:278)
at org.broadinstitute.gatk.tools.walkers.annotator.VariantAnnotatorEngine.annotateContextForActiveRegion(VariantAnnotatorEngine.java:260)
at org.broadinstitute.gatk.tools.walkers.haplotypecaller.HaplotypeCallerGenotypingEngine.annotateCall(HaplotypeCallerGenotypingEngine.java:328)
at org.broadinstitute.gatk.tools.walkers.haplotypecaller.HaplotypeCallerGenotypingEngine.assignGenotypeLikelihoods(HaplotypeCallerGenotypingEngine.java:290)
at org.broadinstitute.gatk.tools.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:962)
at org.broadinstitute.gatk.tools.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:250)
at org.broadinstitute.gatk.engine.traversals.TraverseActiveRegions$TraverseActiveRegionMap.apply(TraverseActiveRegions.java:709)
at org.broadinstitute.gatk.engine.traversals.TraverseActiveRegions$TraverseActiveRegionMap.apply(TraverseActiveRegions.java:705)
at org.broadinstitute.gatk.utils.nanoScheduler.NanoScheduler$ReadMapReduceJob.run(NanoScheduler.java:471)
at java.util.concurrent.Executors$RunnableAdapter.call(Executors.java:511)
at java.util.concurrent.FutureTask.run(FutureTask.java:266)
at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1142)
at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:617)
at java.lang.Thread.run(Thread.java:745)
searching for '100I': the following command returns nothing:
find wordir -name "**.bam" | while read F; do samtools view $F "12:132806991-133841896" | grep 100I ; done
removing 'BaseCountsBySample' seems to fix the problem.
P.