Using hs37d5 genome to call SVs on data aligned with hg19

Steph_UCSteph_UC Member
edited February 2016 in GenomeSTRiP

Hi, I need help please.

SO my question here is does the message below has to do with the fact that I happened to call SVs on hg19 aligned data with the HS37D5 genome? If so how can I fix it without having to realign all my bam files; or run the whole SVPreprocess again (takes a while with my resources).
SVPreprocess worked without trouble, and SVDiscovery ran all the way to this point (almost the end).

Please help, Here is the trace error message.

INFO 22:09:02,218 31-Jan-2016 SVDiscovery - Processing locus: hs37d5:1-35477943:100-1000000
INFO 22:09:02,218 31-Jan-2016 SVDiscovery - Locus search window: hs37d5:1-35477943
WARNING: BAM index file /projects/Families/mybams/person7.bam.bai is older than BAM /projects/Families/mybams/person7.bam
INFO 22:09:03,388 31-Jan-2016 GATKRunReport - Uploaded run statistics report to AWS S3

ERROR ------------------------------------------------------------------------------------------
ERROR stack trace

java.lang.IllegalArgumentException: Invalid reference index -1
at htsjdk.samtools.QueryInterval.(QueryInterval.java:24)
at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.query(SamReader.java:473)
at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.queryContained(SamReader.java:356)
at org.broadinstitute.sv.discovery.DeletionDiscoveryAlgorithm.runTraversal(DeletionDiscoveryAlgorithm.java:141)
at org.broadinstitute.sv.discovery.SVDiscoveryWalker.onTraversalDone(SVDiscoveryWalker.java:107)
at org.broadinstitute.sv.discovery.SVDiscoveryWalker.onTraversalDone(SVDiscoveryWalker.java:42)
at org.broadinstitute.gatk.engine.executive.Accumulator$StandardAccumulator.finishTraversal(Accumulator.java:129)
at org.broadinstitute.gatk.engine.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:116)
at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:319)
at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:121)
at org.broadinstitute.sv.main.SVCommandLine.execute(SVCommandLine.java:124)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:248)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:155)
at org.broadinstitute.sv.main.SVCommandLine.main(SVCommandLine.java:78)
at org.broadinstitute.sv.main.SVDiscovery.main(SVDiscovery.java:21)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version ):
ERROR
ERROR This might be a bug. Please check the documentation guide to see if this is a known problem.
ERROR If not, please post the error message, with stack trace, to the GATK forum.
ERROR Visit our website and forum for extensive documentation and answers to
ERROR commonly asked questions http://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Invalid reference index -1
ERROR ------------------------------------------------------------------------------------------

INFO 22:09:26,161 QGraph - Writing incremental jobs reports...
INFO 22:09:26,162 QGraph - 2 Pend, 0 Run, 1 Fail, 85 Done
INFO 22:09:26,166 QCommandLine - Writing final jobs report...
INFO 22:09:26,166 QCommandLine - Done with errors

Post edited by Steph_UC on

Answers

Sign In or Register to comment.