Question: Error while using hisat2 (hisat2-align exited with value 137)
0
gravatar for szabo.marton
12 months ago by
szabo.marton0 wrote:

Dear All,

I've tried mapping my trimmed RNA-Seq reads against the indexed reference genome using hisat2 but I've got the following error message:

Killed
(ERR): hisat2-align exited with value 137

I coudn't find any answers related to this type of error. Here is the used command:

hisat2 -x Rno_index --known-splicesite-infile rno6_splice_sites.txt -1 S459_R1_paired.fastq.gz -2 S459_R2_paired.fastq.hisat2 -x Rno_index --known-splicesite-infile rno6_splice_sites.txt -1 S459_R1_paired.fastq.gz -2 S459_R2_paired.fastq.gz | samtools view -bS - > SIR_1.bam

I've followed the workflow provided by the Babraham Bioinformatics

For the mapping I've biuld my own indexed genome using hisat2-build from the Rattus_norvegicus.Rnor_6.0.dna.toplevel.fa reference genome, obtained from Ensembl. Also I've tried with the already indexed rat genome from Ensembl, but without any succes.

Before mapping the very first step was merging the individual sequencing runs from the same samples with cat file_1R1.fastq file_2R1.fastq.gz > file.fastq.gz then I've performed adapter trimming:

java -jar /home/Documents/Trimmomatic-0.38/trimmomatic-0.38.jar PE S458_R1.fastq.gz S459_R2.fastq.gz S459_R1_paired.fastq.gz S459_R1_unpaired.fastq.gz S459_R2_paired.fastq.gz S459_R2_unpaired.fastq.gz HEADCROP:15 CROP:55 ILLUMINACLIP:TruSeq3-PE.fa:2:30:10:8:true TRAILING:10 SLIDINGWINDOW:4:20 MINLEN:36

Any idea about what did I wrong are appreciated.

rna-seq alignment • 855 views
ADD COMMENTlink written 12 months ago by szabo.marton0
1

How much memory do you have available? Are you running this under a job scheduler? Your job is getting killed so you are likely exceeding some resource allocation.

ADD REPLYlink modified 12 months ago • written 12 months ago by genomax92k

Thank you for your reply! Would it be that simple? I have only 4 Gb RAM, it is not to much but it was enough for previous analysis. I don't use any job scheduler, would it be useful?

ADD REPLYlink written 12 months ago by szabo.marton0
1

4 GB RAM is not enough to be able to align rat genome with pretty much any NGS aligner.

ADD REPLYlink written 12 months ago by genomax92k

Could it help if I double the RAM or should I use a better computer?

ADD REPLYlink written 12 months ago by szabo.marton0
1

Find better hardware if you can. If you can't then doubling RAM may allow you to use bwa.

ADD REPLYlink modified 12 months ago • written 12 months ago by genomax92k

I've double the RAM and now it's working perfectly

ADD REPLYlink written 12 months ago by szabo.marton0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1529 users visited in the last hour