Quantcast
Viewing all articles
Browse latest Browse all 12345

Invalid intervals error

Hi, I'm currently working with Genomestrip (version 3.7.GS-r1748-0-g74bfe0b) to do some structural variant calling on some BAMs that were aligned to hg38. Naturally, I used the resources available from the hg38 prerelease bundle for Genomestrip available here: ftp://ftp.broadinstitute.org/pub/svtoolkit/hg38_pre_release/

So far, whenever I use Genomestrip with hg38, I get several invalid interval errors. These appears to be legitimate intervals, so I'm not sure how Genomestrip is failing on these. Here are the intervals it fails on:

chr22:24317068-24322453
chrY:10022097-10027150
chr21:20211411-20216694

And here's a stack trace that's representative of these errors:

##### ERROR --
##### ERROR stack trace 
java.lang.RuntimeException: Invalid interval: chrY:16566323-16571422
        at org.broadinstitute.sv.metadata.depth.ReadCountFileReader$ReadCountDataIterator.<init>(ReadCountFileReader.java:415)
        at org.broadinstitute.sv.metadata.depth.ReadCountFileReader.getReadCountDataIteratorInternal(ReadCountFileReader.java:313)
        at org.broadinstitute.sv.metadata.depth.ReadCountFileReader.getReadCountDataIterator(ReadCountFileReader.java:282)
        at org.broadinstitute.sv.metadata.depth.ReadCountFileReader.getReadCounts(ReadCountFileReader.java:257)
        at org.broadinstitute.sv.common.ReadCountCache.getReadCounts(ReadCountCache.java:100)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.computeRefReadCounts(GenotypingDepthModule.java:283)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.computeRefReadCounts(GenotypingDepthModule.java:240)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.getReadCounts(GenotypingDepthModule.java:225)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.getCnpReadCounts(GenotypingDepthModule.java:212)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.genotypeCnp(GenotypingDepthModule.java:136)
        at org.broadinstitute.sv.genotyping.GenotypingDepthModule.genotypeCnp(GenotypingDepthModule.java:57)
        at org.broadinstitute.sv.genotyping.GenotypingAlgorithm.genotypeCnpInternal(GenotypingAlgorithm.java:149)
        at org.broadinstitute.sv.genotyping.GenotypingAlgorithm.genotypeCnp(GenotypingAlgorithm.java:114)
        at org.broadinstitute.sv.genotyping.SVGenotyperWalker.processVCFFile(SVGenotyperWalker.java:276)
        at org.broadinstitute.sv.genotyping.SVGenotyperWalker.map(SVGenotyperWalker.java:220)
        at org.broadinstitute.sv.genotyping.SVGenotyperWalker.map(SVGenotyperWalker.java:57)
        at org.broadinstitute.gatk.engine.traversals.TraverseLociNano.traverse(TraverseLociNano.java:106)
        at org.broadinstitute.gatk.engine.traversals.TraverseLociNano.traverse(TraverseLociNano.java:48)
        at org.broadinstitute.gatk.engine.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:98)
        at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:316)
        at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:123)
        at org.broadinstitute.sv.main.SVCommandLine.execute(SVCommandLine.java:133)
        at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:256)
        at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:158)
        at org.broadinstitute.sv.main.SVCommandLine.main(SVCommandLine.java:87)
        at org.broadinstitute.sv.main.SVGenotyper.main(SVGenotyper.java:21)
##### ERROR ------------------------------------------------------------------------------------------

Any help would be greatly appreciated!


Viewing all articles
Browse latest Browse all 12345

Trending Articles



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