Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Mike_Brown
    Member
    • Jan 2016
    • 13

    #1

    A GATK RUNTIME ERROR has occurred

    I am new to GATK. Today, I run
    Code:
    java -jar /media/zwang10/Elements/UK10K/GenomeAnalysisTK-3.5/GenomeAnalysisTK.jar -T CountLoci -R /media/zwang10/Elements/UK10K/human_g1k_v37.fasta -I _EGAR00001038931_36843.pe.raw.sorted.bam
    But I got
    Code:
    INFO  14:59:39,061 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  14:59:39,064 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.5-0-g36282e4, Compiled 2015/11/25 04:03:56 
    INFO  14:59:39,064 HelpFormatter - Copyright (c) 2010 The Broad Institute 
    INFO  14:59:39,064 HelpFormatter - For support and documentation go to http://www.broadinstitute.org/gatk 
    INFO  14:59:39,068 HelpFormatter - Program Args: -T CountReads -R /media/zwang10/Elements/UK10K/human_g1k_v37.fasta -I _EGAR00001038931_36843.pe.raw.sorted.bam 
    INFO  14:59:39,076 HelpFormatter - Executing as zwang10@office on Linux 3.13.0-74-generic amd64; OpenJDK 64-Bit Server VM 1.7.0_91-b02. 
    INFO  14:59:39,077 HelpFormatter - Date/Time: 2016/01/03 14:59:39 
    INFO  14:59:39,077 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  14:59:39,077 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  14:59:39,138 GenomeAnalysisEngine - Strictness is SILENT 
    INFO  14:59:39,342 GenomeAnalysisEngine - Downsampling Settings: No downsampling 
    INFO  14:59:39,354 SAMDataSource$SAMReaders - Initializing SAMRecords in serial 
    INFO  14:59:39,416 SAMDataSource$SAMReaders - Done initializing BAM readers: total time 0.06 
    INFO  14:59:42,897 GATKRunReport - Uploaded run statistics report to AWS S3 
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR stack trace 
    org.broadinstitute.gatk.utils.exceptions.ReviewedGATKException: Unable to open index file (mark/reset not supported)_EGAR00001038931_36843.pe.raw.sorted.bam.bai
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.openIndexFile(GATKBAMIndex.java:319)
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.<init>(GATKBAMIndex.java:97)
    	at org.broadinstitute.gatk.engine.datasources.reads.SAMDataSource.<init>(SAMDataSource.java:388)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.createReadsDataSource(GenomeAnalysisEngine.java:931)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.initializeDataSources(GenomeAnalysisEngine.java:822)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:286)
    	at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:121)
    	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.gatk.engine.CommandLineGATK.main(CommandLineGATK.java:106)
    Caused by: java.io.IOException: mark/reset not supported
    	at java.io.InputStream.reset(InputStream.java:347)
    	at htsjdk.samtools.SamIndexes.asBaiSeekableStreamOrNull(SamIndexes.java:80)
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.openIndexFile(GATKBAMIndex.java:314)
    	... 9 more
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR A GATK RUNTIME ERROR has occurred (version 3.5-0-g36282e4):
    ##### 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: Unable to open index file (mark/reset not supported)_EGAR00001038931_36843.pe.raw.sorted.bam.bai
    ##### ERROR ------------------------------------------------------------------------------------------
    zwang10@office:/media/zwang10/Elements/UK10K$ java -jar /media/zwang10/Elements/UK10K/GenomeAnalysisTK-3.5/GenomeAnalysisTK.jar -T CountReads -R /media/zwang10/Elements/UK10K/human_g1k_v37.fasta -I _EGAR00001038931_36843.pe.raw.sorted.bam > temp
    INFO  15:08:54,928 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  15:08:54,931 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.5-0-g36282e4, Compiled 2015/11/25 04:03:56 
    INFO  15:08:54,931 HelpFormatter - Copyright (c) 2010 The Broad Institute 
    INFO  15:08:54,931 HelpFormatter - For support and documentation go to http://www.broadinstitute.org/gatk 
    INFO  15:08:54,935 HelpFormatter - Program Args: -T CountReads -R /media/zwang10/Elements/UK10K/human_g1k_v37.fasta -I _EGAR00001038931_36843.pe.raw.sorted.bam 
    INFO  15:08:54,943 HelpFormatter - Executing as zwang10@office on Linux 3.13.0-74-generic amd64; OpenJDK 64-Bit Server VM 1.7.0_91-b02. 
    INFO  15:08:54,943 HelpFormatter - Date/Time: 2016/01/03 15:08:54 
    INFO  15:08:54,943 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  15:08:54,944 HelpFormatter - -------------------------------------------------------------------------------- 
    INFO  15:08:55,005 GenomeAnalysisEngine - Strictness is SILENT 
    INFO  15:08:55,171 GenomeAnalysisEngine - Downsampling Settings: No downsampling 
    INFO  15:08:55,181 SAMDataSource$SAMReaders - Initializing SAMRecords in serial 
    INFO  15:08:55,237 SAMDataSource$SAMReaders - Done initializing BAM readers: total time 0.06 
    INFO  15:08:57,411 GATKRunReport - Uploaded run statistics report to AWS S3 
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR stack trace 
    org.broadinstitute.gatk.utils.exceptions.ReviewedGATKException: Unable to open index file (mark/reset not supported)_EGAR00001038931_36843.pe.raw.sorted.bam.bai
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.openIndexFile(GATKBAMIndex.java:319)
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.<init>(GATKBAMIndex.java:97)
    	at org.broadinstitute.gatk.engine.datasources.reads.SAMDataSource.<init>(SAMDataSource.java:388)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.createReadsDataSource(GenomeAnalysisEngine.java:931)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.initializeDataSources(GenomeAnalysisEngine.java:822)
    	at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:286)
    	at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:121)
    	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.gatk.engine.CommandLineGATK.main(CommandLineGATK.java:106)
    Caused by: java.io.IOException: mark/reset not supported
    	at java.io.InputStream.reset(InputStream.java:347)
    	at htsjdk.samtools.SamIndexes.asBaiSeekableStreamOrNull(SamIndexes.java:80)
    	at org.broadinstitute.gatk.engine.datasources.reads.GATKBAMIndex.openIndexFile(GATKBAMIndex.java:314)
    	... 9 more
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR A GATK RUNTIME ERROR has occurred (version 3.5-0-g36282e4):
    ##### 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: Unable to open index file (mark/reset not supported)_EGAR00001038931_36843.pe.raw.sorted.bam.bai
    ##### ERROR ------------------------------------------------------------------------------------------
    I did check _EGAR00001038931_36843.pe.raw.sorted.bam.bai is in the same folder.
    What happened?
  • blancha
    Senior Member
    • May 2013
    • 367

    #2
    I have no idea, but I wouldn't start a filename with an underscore.

    Comment

    • Mike_Brown
      Member
      • Jan 2016
      • 13

      #3
      Originally posted by blancha View Post
      I have no idea, but I wouldn't start a filename with an underscore.
      It is the original name of the file. Do I need to change it?

      Comment

      • GenoMax
        Senior Member
        • Feb 2008
        • 7142

        #4
        Certainly worth trying.

        Remember to rename both bam and index files. Not sure if the file name is embedded in the index, otherwise you may need to re-index the bam after changing the name.

        Comment

        • Mike_Brown
          Member
          • Jan 2016
          • 13

          #5
          Originally posted by GenoMax View Post
          Certainly worth trying.

          Remember to rename both bam and index files. Not sure if the file name is embedded in the index, otherwise you may need to re-index the bam after changing the name.
          I changed name as EGAR00001038931.36843.pe.raw.sorted.bam and EGAR00001038931.36843.pe.raw.sorted.bam.bai


          I still get the same problem.

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
            by SEQadmin2



            CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

            Despite this, “CRISPR helped turn genome editing from a specialized technique into
            ...
            07-31-2026, 11:01 AM
          • SEQadmin2
            Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
            by SEQadmin2


            Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

            The systematic characterization of the human proteome has
            ...
            07-20-2026, 11:48 AM
          • SEQadmin2
            Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
            by SEQadmin2



            Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
            ...
            07-09-2026, 11:10 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Today, 07:41 AM
          0 responses
          9 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 08-03-2026, 10:13 AM
          0 responses
          25 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-31-2026, 02:55 AM
          0 responses
          38 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-24-2026, 12:17 PM
          0 responses
          25 views
          0 reactions
          Last Post SEQadmin2  
          Working...