problem with bowtie2 allignment
0
0
Entering edit mode
5.9 years ago
arshil • 0

Can any help me out with this problem. I am using bowtie2 to allign my fastq files to reference genome hg38. every time when I run the following Command i am getting and error (ERR): bowtie2-align died with signal 9 (KILL)".

command I am using :

bowtie2 -p 24 -x hg38 -X 2000 -1 R-CA-013-2-14_S3_R1.fastq.gz -2 R-CA-013-2-14_S3_R2.fastq.gz -S R-CA-013-2-14_S3 --no-unal

I would really appreciate if anyone could let me know whats the error about? sorry for the english. thanks

RNA-Seq bowtie2 • 2.3k views
ADD COMMENT
0
Entering edit mode

How much memory do you have available? And did you try with fewer threads, like -p 6?

ADD REPLY
0
Entering edit mode

Total of 50 gb memory and I tried using p as 6, still Getting the same error!!

ADD REPLY
0
Entering edit mode

Was bowtie2 already installed on this cluster (and know to work) or did you download the program (and compile?) yourself?

ADD REPLY
0
Entering edit mode

It was already installed on cluster! All I did was module load bowtie2 And wrote my command

ADD REPLY
0
Entering edit mode

Are you also using a job scheduling program then? If you are on a cluster you may not be allowed to run programs outside the scheduler. Is the above command only thing you typed in? You may need to talk with your local IT support to find out what is the proper way of submitting jobs on your cluster.

ADD REPLY
0
Entering edit mode

Hi! So I am running this command on HPC cluster which has 32 cores!!

ADD REPLY
0
Entering edit mode

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized.

@h.mon asked how much "memory" you have available?

ADD REPLY
0
Entering edit mode

You can try to download the latest binary and run the command with it instead of the preinstalled one.

ADD REPLY
0
Entering edit mode

Did you build a reference genome index?

It shouldn't matter I believe, but the argument to -S is the name of your sam file so would be better to use R-CA-013-2-14_S3.sam.

Even better would be to let bowtie write to stdout and convert to bam on the fly (requiring a reasonably recent version of samtools):

bowtie2 -p 24 -x hg38 -X 2000 -1 R-CA-013-2-14_S3_R1.fastq.gz -2 R-CA-013-2-14_S3_R2.fastq.gz --no-unal | satmtools sort -o R-CA-013-2-14_S3.bam
ADD REPLY

Login before adding your answer.

Traffic: 2350 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6