Hi,

I have started using GSNAP to align some RNA-Seq reads to a reference genome, and am having some trouble getting the program to run properly. The alignment works fine if I do not use the -s parameter to provide splicing information, but if I do try and use it, the program almost immediately quits for a reason I cannot particularly understand.

Normally I run the command like this:

• gsnap --gunzip -D <directory to="" database=""> -d <database> -t 12 --orientation=FR -A sam > output.sam -N 1 Read1.fastq.gz Read2.fastq.gz

However, when I add splicing (the file is called splicesites.iit and is located in the maps directory of the genome database):

• gsnap --gunzip -D <directory to="" database=""> -d <database> -t 12 --orientation=FR -A sam > output.sam -N 1 -s splicesites Read1.fastq.gz Read2.fastq.gz

the program runs only to the point of allocating memory to everything and then aborts as follows before even starting the alignment:

Allocating memory for ref positions, kmer 8, interval 1...Attached new memory for GSNAP_copy/GSNAP/X.laevis_database/X.laevis_database.ref081locpositions...done (4,891,849,244 bytes, 402.54 sec)

Not expecting version 1

Signal received: SIGABRT

Calling Access_emergency_cleanup

Removed existing memory for shmid 950287

Removed existing memory for shmid 917518

Removed existing memory for shmid 884749

Removed existing memory for shmid 851980

Removed existing memory for shmid 819211

Removed existing memory for shmid 786442

Removed existing memory for shmid 753673

Removed existing memory for shmid 720897

Removed existing memory for shmid 688128

any advice on what might be going on? I am unsure how adding splicing would lead to an abort error like this, but have been unable to find any information in the documentation. I also am fairly amateur when it comes to CS, so its possible I am screwing up something simple with my commands and do not realize. Any help would be appreciated!

Thanks!