Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

Error in running SplitNCigarReads - htsjdk.samtools.util.RuntimeIOException

Answered
0

9 comments

  • Avatar
    KATHIE NGO

    Hi ,

    I also tried using the latest version GATK v4.1.8.1, and I still see the same issues. Any advice and guidance is appreciated.

    Thank you for your time and help.

    Best Regards,

    Kathie

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    Hi KATHIE NGO, this error is resulting from the program running out of space to write files. You can specify a temporary directory that has more space with the option --tmp-dir. 

    0
    Comment actions Permalink
  • Avatar
    KATHIE NGO

    Hi Genevieve, 

    Thank you for your time and help. Can you please tell me where does the program write files by default, is it where the command was executed? The output BAM files are empty and we have about 2.2TB of free space in our /home directory on the server where I executed the command.

    Best Regards,

    Kathie

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    Hi KATHIE NGO, it depends on your system. Some systems use an automatic temporary space that has less space than necessary for some jobs. If the --tmp-dir option does not work, please let us know so we can troubleshoot from there.

    0
    Comment actions Permalink
  • Avatar
    KATHIE NGO

    Hi Genevieve, 

    Thank you for your time and help. Specifying the option --tmp-dir solved my issue. 

    Best Regards,

    Kathie

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    Thanks for the update! Glad it worked. 

    0
    Comment actions Permalink
  • Avatar
    JAY SINGH

    Hello Genevieve,

    I am using GATK 4.2.5.0. I am getting the same issue.

    Unfortunately, my issue is not solved after specifying --tmp-dir.

    Any other suggestion?

    Thanks and regards,

    Jay

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    Hi JAY SINGH,

    Could you give more details about your issue, including the command and the error message?

    Thank you,

    Genevieve

    0
    Comment actions Permalink
  • Avatar
    Niharika

    Hi

    I am using GATK SplitNCigarReads version 4.0.5.1

    And getting an error (attached)

    gatk SplitNCigarReads -R SplitNCigar_Reads/human_genome.fna -I Mark_duplicates/PicardOutput/sorted_unique_reads_withRGAlignment/SRR1rtgtrr_Unique_reads_WithRG.bam -O SplitNCigar_Reads/output/SRR1rtgtrr_SplitNCigar_out.bam --sites-only-vcf-output true --TMP_DIR TMP_DIR
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:80: SyntaxWarning: "is" with a literal. Did you mean "=="?
      if len(args) is 0 or (len(args) is 1 and (args[0] == "--help" or args[0] == "-h")):
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:80: SyntaxWarning: "is" with a literal. Did you mean "=="?
      if len(args) is 0 or (len(args) is 1 and (args[0] == "--help" or args[0] == "-h")):
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:117: SyntaxWarning: "is" with a literal. Did you mean "=="?
      if len(args) is 1 and args[0] == "--list":
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:301: SyntaxWarning: "is" with a literal. Did you mean "=="?
      if call(["gsutil", "-q", "stat", gcsjar]) is 0:
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:305: SyntaxWarning: "is" with a literal. Did you mean "=="?
      if call(["gsutil", "cp", jar, gcsjar]) is 0:
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:458: SyntaxWarning: "is not" with a literal. Did you mean "!="?
      if not len(properties) is 0:
    /DATA/RamanDeep/Tools/anaconda3/bin/gatk:462: SyntaxWarning: "is not" with a literal. Did you mean "!="?
      if not len(filesToAdd) is 0:
    Using GATK jar /DATA/RamanDeep/Tools/anaconda3/share/gatk4-4.0.5.1-0/gatk-package-4.0.5.1-local.jar
    Running:
        java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -jar /DATA/RamanDeep/Tools/anaconda3/share/gatk4-4.0.5.1-0/gatk-package-4.0.5.1-local.jar SplitNCigarReads -R SplitNCigar_Reads/human_genome.fna -I Mark_duplicates/PicardOutput/sorted_unique_reads_withRGAlignment/SRR1rtgtrr_Unique_reads_WithRG.bam -O SplitNCigar_Reads/output/SRR1rtgtrr_SplitNCigar_out.bam --sites-only-vcf-output true --TMP_DIR TMP_DIR
    18:44:57.015 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/DATA/RamanDeep/Tools/anaconda3/share/gatk4-4.0.5.1-0/gatk-package-4.0.5.1-local.jar!/com/intel/gkl/native/libgkl_compression.so
    18:44:57.234 INFO  SplitNCigarReads - ------------------------------------------------------------
    18:44:57.234 INFO  SplitNCigarReads - The Genome Analysis Toolkit (GATK) v4.0.5.1
    18:44:57.235 INFO  SplitNCigarReads - For support and documentation go to https://software.broadinstitute.org/gatk/
    18:44:57.235 INFO  SplitNCigarReads - Executing as niharika@user on Linux v4.15.0-213-generic amd64
    18:44:57.235 INFO  SplitNCigarReads - Java runtime: OpenJDK 64-Bit Server VM v1.8.0_382-b05
    18:44:57.235 INFO  SplitNCigarReads - Start Date/Time: 30 November, 2023 6:44:57 PM IST
    18:44:57.235 INFO  SplitNCigarReads - ------------------------------------------------------------
    18:44:57.235 INFO  SplitNCigarReads - ------------------------------------------------------------
    18:44:57.235 INFO  SplitNCigarReads - HTSJDK Version: 2.15.1
    18:44:57.235 INFO  SplitNCigarReads - Picard Version: 2.18.2
    18:44:57.235 INFO  SplitNCigarReads - HTSJDK Defaults.COMPRESSION_LEVEL : 2
    18:44:57.235 INFO  SplitNCigarReads - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
    18:44:57.235 INFO  SplitNCigarReads - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
    18:44:57.235 INFO  SplitNCigarReads - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
    18:44:57.236 INFO  SplitNCigarReads - Deflater: IntelDeflater
    18:44:57.236 INFO  SplitNCigarReads - Inflater: IntelInflater
    18:44:57.236 INFO  SplitNCigarReads - GCS max retries/reopens: 20
    18:44:57.236 INFO  SplitNCigarReads - Using google-cloud-java patch 6d11bef1c81f885c26b2b56c8616b7a705171e4f from https://github.com/droazen/google-cloud-java/tree/dr_all_nio_fixes
    18:44:57.236 INFO  SplitNCigarReads - Initializing engine
    18:44:57.686 INFO  SplitNCigarReads - Done initializing engine
    18:44:57.729 INFO  SplitNCigarReads - Shutting down engine
    [30 November, 2023 6:44:57 PM IST] org.broadinstitute.hellbender.tools.walkers.rnaseq.SplitNCigarReads done. Elapsed time: 0.01 minutes.
    Runtime.totalMemory()=2175270912
    htsjdk.samtools.util.RuntimeIOException: Error opening file: file:///DATA1/Niharika/NSCLC_DATA_ANALYSIS/RE_ANALYSIS_RNA_SEQ/SplitNCigar_Reads/output/SRR1rtgtrr_SplitNCigar_out.bam
        at htsjdk.samtools.SAMFileWriterFactory.makeBAMWriter(SAMFileWriterFactory.java:302)
        at htsjdk.samtools.SAMFileWriterFactory.makeBAMWriter(SAMFileWriterFactory.java:259)
        at htsjdk.samtools.SAMFileWriterFactory.makeSAMOrBAMWriter(SAMFileWriterFactory.java:435)
        at htsjdk.samtools.SAMFileWriterFactory.makeWriter(SAMFileWriterFactory.java:474)
        at org.broadinstitute.hellbender.utils.read.ReadUtils.createCommonSAMWriterFromFactory(ReadUtils.java:1183)
        at org.broadinstitute.hellbender.utils.read.ReadUtils.createCommonSAMWriter(ReadUtils.java:1131)
        at org.broadinstitute.hellbender.engine.GATKTool.createSAMWriter(GATKTool.java:772)
        at org.broadinstitute.hellbender.tools.walkers.rnaseq.SplitNCigarReads.onTraversalStart(SplitNCigarReads.java:170)
        at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:992)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:135)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:180)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:199)
        at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160)
        at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203)
        at org.broadinstitute.hellbender.Main.main(Main.java:289)
    Caused by: java.nio.file.NoSuchFileException: SplitNCigar_Reads/output/SRR1rtgtrr_SplitNCigar_out.bam
        at sun.nio.fs.UnixException.translateToIOException(UnixException.java:86)
        at sun.nio.fs.UnixException.rethrowAsIOException(UnixException.java:102)
        at sun.nio.fs.UnixException.rethrowAsIOException(UnixException.java:107)
        at sun.nio.fs.UnixFileSystemProvider.newByteChannel(UnixFileSystemProvider.java:214)
        at java.nio.file.spi.FileSystemProvider.newOutputStream(FileSystemProvider.java:434)
        at java.nio.file.Files.newOutputStream(Files.java:216)
        at htsjdk.samtools.SAMFileWriterFactory.makeBAMWriter(SAMFileWriterFactory.java:290)
        ... 14 more
    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk