Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

HaplotypeCaller raises an error with -A BaseCountsBySample

$
0
0

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.


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>