JUM (v2.0.2 and up) stops and get killed at a certain step
Dear all, I'm using JUM pipeline (https://github.com/qqwang-berkeley/JUM/wiki/3.1.-Manual-running-JUM-(v2.0.2-and-up)) to analyze alternative splicing RNA-sequencing but the bash code stops and get killed at this step:
process_sortedbam () {
sortfile=$1;
bamfiletobed=${sortfile%_sorted.bam}.bed;
bamsortedbed=${sortfile%_sorted.bam}.sorted.bed;
intersect_bam_long_intron=${sortfile%_sorted.bam}"_"intersect_long_intron.txt;
intersect_bam_short_intron1=${sortfile%_sorted.bam}"_"intersect_short_intron1.txt;
intersect_bam_short_intron2=${sortfile%_sorted.bam}"_"intersect_short_intron2.txt;
bedtools bamtobed -i $sortfile > $bamfiletobed;
sort -k1,1 -k2,2n $bamfiletobed > $bamsortedbed;
intersectBed -a $bamsortedbed -b output_long_intron_sorted.gff -sorted -wa -u -f 1 > $intersect_bam_long_intron
intersectBed -a $bamsortedbed -b output_short_intron1_sorted.gff -sorted -wa -u -f 1 > $intersect_bam_short_intron1
intersectBed -a $bamsortedbed -b output_short_intron2_sorted.gff -sorted -wa -u -F 1 > $intersect_bam_short_intron2
}
Exact error:
Killed sort -k1,1 -k2,2n $bamfiletobed > $bamsortedbed
Has anyone have any experience with this error or know what should I do? I'm guessing it might be a CPU error as it gets killed but more than that, not sure. I'm working with a 16 threads system and once used all and other times just 14 of them, but it gets killed anyway. I should mention that JUM runs all the samples (6 in my case) in parallel.
Code I'm using:
bash /path/JUM_2.0.2_multi/JUM_A.sh --Folder /path/JUM_2.0.2_multi --JuncThreshold 5 --Condition1_fileNum_threshold 2 --Condition2_fileNum_threshold 2 --IRthreshold 5 --Readlength 101 --Thread 16 --Condition1SampleName cntrl1,cntrl2,cntrl3 --Condition2SampleName treat1,treat2,treat3
Thanks.
• 1,662 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Look into if you are running out of RAM. Using less cores may actually help in that case since RAM requirements may be reduced.
I reduced to 8, but it didn't work. Still get killed at the same step. The developers told me to run the script on non-parallel mode but there isn't any bash code or other guide that I can switch it to that mode.
Try
--Thread 1then.They sent a non parallel code as it was analyzing all my samples in parallel mode. My system wasn't able to do that.