Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Gsnap error - possibly very basic

    When I put in a basic input, it terminates fairly quickly:

    gsnap -d hg19 -N 1 InputR1.fq InputR2.fq
    GSNAP version 2015-06-12 called with args: gsnap -d hg19 -N 1 InputR1.fq InputR2.fq
    Novel splicing (-N) turned on => assume reads are RNA-Seq
    Checking compiler assumptions for popcnt: 000041A7 __builtin_clz=17 __builtin_ctz=0 _mm_popcnt_u32=7 __builtin_popcount=7
    Checking compiler assumptions for SSE2: 000041A7 10D63AF1 xor=10D67B56
    Checking compiler assumptions for SSE4.1: -89 -15 max=241 => compiler zero extends
    Finished checking compiler assumptions
    Note: >1 sequence detected, so index files are being memory mapped.
    GSNAP can run slowly at first while the computer starts to accumulate
    pages from the hard disk into its cache. To copy index files into RAM
    instead of memory mapping, use -B 3, -B 4, or -B 5, if you have enough RAM.
    For more speed, also try multiple threads (-t <int>), if you have multiple processors or cores.
    Allocating memory for compressed genome (oligos)...Error with shmget. Error 2: No such file or directory
    Signal received: SIGABRT
    Calling Access_emergency_cleanup
    Unable to retrieve queryseq for request

    When I add -B 3, it gets a little farther:

    gsnap -d hg19 -N 1 -B 3 InputR1.fq InputR2.fq
    GSNAP version 2015-06-12 called with args: gsnap -d hg19 -N 1 -B 3 InputR1.fq InputR2.fq
    Novel splicing (-N) turned on => assume reads are RNA-Seq
    Checking compiler assumptions for popcnt: 000041A7 __builtin_clz=17 __builtin_ctz=0 _mm_popcnt_u32=7 __builtin_popcount=7
    Checking compiler assumptions for SSE2: 000041A7 10D63AF1 xor=10D67B56
    Checking compiler assumptions for SSE4.1: -89 -15 max=241 => compiler zero extends
    Finished checking compiler assumptions
    Note: >1 sequence detected, so index files are being memory mapped.
    GSNAP can run slowly at first while the computer starts to accumulate
    pages from the hard disk into its cache. To copy index files into RAM
    instead of memory mapping, use -B 3, -B 4, or -B 5, if you have enough RAM.
    For more speed, also try multiple threads (-t <int>), if you have multiple processors or cores.
    Pre-loading compressed genome (oligos)...,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,.,..,.,..,.,..,.,....,.,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,......,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,.,.,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,..done (1,093,035,816 bytes, 266855 pages, 0.10 sec)
    Pre-loading compressed genome (bits)...,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,.,..,.,..,.,..,.,....,.,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,......,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,..,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,.,.,.,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,...,..done (1,093,035,840 bytes, 266855 pages, 0.10 sec)
    Allocating memory for indexij ptrs...Error with shmget. Error 2: No such file or directory
    Signal received: SIGABRT
    Calling Access_emergency_cleanup
    Unable to retrieve queryseq for request

    I have little idea what's going wrong. I'm going to keep tinkering, but I was hoping someone here would have some idea of what to do.

  • #2
    Looks like the program is not able to find/read the input sequence file(s). Can you try providing full paths to InputR1.fq InputR2.fq?

    You may also be getting an error (Error with shmget) with memory allocation in the second try.

    Comment


    • #3
      I just had a similar error using the same version.

      gsnap -d mydatabase --gunzip -B 4 -t 20 -m 8 -N 1 --force-single-end ConditionAlpha.fastq.gz

      GSNAP version 2015-06-12 called with args: gsnap -d mydatabase --gunzip -B 4 -t 20 -m 8 -N 1 --force-single-end ConditionAlpha.fastq.gz
      Novel splicing (-N) turned on => assume reads are RNA-Seq
      Checking compiler assumptions for popcnt: 6B8B4567 __builtin_clz=1 __builtin_ctz=0 _mm_popcnt_u32=17 __builtin_popcount=17
      Checking compiler assumptions for SSE2: 6B8B4567 327B23C6 xor=59F066A1
      Checking compiler assumptions for SSE4.1: -103 -58 max=198 => compiler zero extends
      Finished checking compiler assumptions
      Note: >1 sequence detected, so index files are being memory mapped.
      GSNAP can run slowly at first while the computer starts to accumulate
      pages from the hard disk into its cache. To copy index files into RAM
      instead of memory mapping, use -B 3, -B 4, or -B 5, if you have enough RAM.
      For more speed, also try multiple threads (-t <int>), if you have multiple processors or cores.
      Allocating memory for compressed genome (oligos)...Error with shmget. Error 2: No such file or directory
      Signal received: SIGABRT
      Calling Access_emergency_cleanup


      Would this be a configuration or access rights error? I don't have admin rights on that server and it is a local installation.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Recent Advances in Sequencing Analysis Tools
        by seqadmin


        The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
        Today, 07:48 AM
      • seqadmin
        Essential Discoveries and Tools in Epitranscriptomics
        by seqadmin




        The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
        04-22-2024, 07:01 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Today, 07:17 AM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-02-2024, 08:06 AM
      0 responses
      19 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-30-2024, 12:17 PM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-29-2024, 10:49 AM
      0 responses
      28 views
      0 likes
      Last Post seqadmin  
      Working...
      X