Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Invalid sequence number 0 in index file

lucasllucasl Planet EarthMember

Hi, I've been using the RealignerTargetCreator tool with a pair of BAMs representing a tumor/normal pair to do joint indel realignment, and I get the following error message:

ERROR --
ERROR stack trace

org.broadinstitute.gatk.utils.exceptions.ReviewedGATKException: Invalid sequence number 0 in index file /sc/orga/scratch/xlochol01/liver/sam/SAMEA1324465_1_recal_fix_sorted_rg.bam.bai
at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.readReferenceSequence(GATKBAMIndex.java:127)
at org.broadinstitute.gatk.engine.datasources.reads.BAMSchedule.(BAMSchedule.java:105)
at org.broadinstitute.gatk.engine.datasources.reads.BAMScheduler.getNextOverlappingBAMScheduleEntry(BAMScheduler.java:296)
at org.broadinstitute.gatk.engine.datasources.reads.BAMScheduler.advance(BAMScheduler.java:185)
at org.broadinstitute.gatk.engine.datasources.reads.BAMScheduler.populateFilteredIntervalList(BAMScheduler.java:105)
at org.broadinstitute.gatk.engine.datasources.reads.BAMScheduler.createOverMappedReads(BAMScheduler.java:67)
at org.broadinstitute.gatk.engine.datasources.reads.IntervalSharder.shardOverMappedReads(IntervalSharder.java:55)
at org.broadinstitute.gatk.engine.datasources.reads.SAMDataSource.createShardIteratorOverMappedReads(SAMDataSource.java:1217)
at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.getShardStrategy(GenomeAnalysisEngine.java:649)
at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:307)
at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:113)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:255)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:157)
at org.broadinstitute.gatk.engine.CommandLineGATK.main(CommandLineGATK.java:108)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 3.6-0-g89b7209):
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 https://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Invalid sequence number 0 in index file /sc/orga/scratch/xlochol01/liver/sam/SAMEA1324465_1_recal_fix_sorted_rg.bam.bai
ERROR ------------------------------------------------------------------------------------------

Any help with this would be appreciated.

Answers

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @lucasl
    Hi,

    This error usually is due to incompatible reference and BAM files. Can you please post your FASTA.dict file and the BAM header? Please also try deleting the BAM index and re-generating it.

    Thanks,
    Sheila

Sign In or Register to comment.