Hello, I am reletively new to bioinformatics and I'm hoping someone will be able to answer my question.
I am trying to use Rainbow to assemble RADseq (short paired-end) reads and eventually call for SNPs. However, I get exit status 139 and the error file says: "/var/spool/PBS/mom_priv/jobs/13307.launch.hpc.SC: line 15: 5705 Segmentation fault rainbow merge -a -i cat.rbdiv.out -o cat.rbasm.out -N 500"
This appears to be a problem with the merging step of rainbow. I get this error with both trimmed and raw reads, but other people are using rainbow on our HPC without any problems so I'm assuming it has something to do with my data...
My script:
module load app/rainbow
rainbow cluster -1 forward -2 reverse > cat.rbcluster.out 2> log
rainbow div -i cat.rbcluster.out -o cat.rbdiv.out
rainbow merge -a -i cat.rbdiv.out -o cat.rbasm.out -N 500
Has anyone seen this error or have any suggestions on how to fix this?
Thanks in advance
•
written
by
phair.nikki •
1 answer
This happened to me because my input file was formatted incorrectly.
Log in to answer this question.
-
mrbayes segmentation fault
written by dgrace999 •
-
IDBA-UD fails due to segmentation error
written by langziv •
-
Rename sample names in g.vcf files fails [Error Segmentation fault (core dumped) ]
written by Alewa
-
Segmentation fault using gemma
written by ggman
-
Samtools segmentation fault when extracting exons from bam file.
written by jstevenson2256 •
-
PBS scripts for CD-HIT
written by lucia •
-
Segmentation fault error with hera
written by P •
-
Bowtie2 throws error no. 139
written by anu014
-
Fastx quality filter
written by kriti.awasthi23 •
-
Cufflinks segmentation fault
written by snp87 •
Looks like you are using PBS scheduler. Are you assigning adequate memory for this job (assuming the command you are using for
rainboware correct)?My experience is Torque / PBS would explicit say if memory was insufficient. But why are you looking at
/var/spool/PBS/mom_priv? What are the messages from stdout and stderr (which, if memory servers me right, you set with#PBS -oand#PBS -e)?Today, when I was running the Rainbow on the cluster I noticed the similar issue. I have assigned 600 GB memory still it fails ("segmentation fault"). Now I am suspecting whether Rainbow requires few TBs of RAM at "merge" step.