HISAT2 Error: Encountered internal HISAT2 exception (#1)
0
0
Entering edit mode
4.3 years ago
ja4123 ▴ 20

Hello, I am using HISAT2 for aligning transcripts to the reference genome. I have build the indexes first. And try the command as follows:

hisat2-build -p 8 --dta -x name_tran -1 input1.fastq.gz -2 input2.fastq.gz -S out.sam

However, I am getting the following error:

Error: Encountered internal HISAT2 exception (#1)

Kindly help.

hisat2 • 5.6k views
ADD COMMENT
1
Entering edit mode

Command you posted has a typo/is incorrect (hisat2-bulit -p 8 --dta -x name_tran -1 input1.fastq.gz -2 input2.fastq.gz -S out.sam). It is not a command that you would use for the alignment.

ADD REPLY
0
Entering edit mode

I am useing this procedure:

Align the RNA-seq reads to the genome
Map the reads for each sample to the reference genome:
$ hisat2 -p 8 --dta -x chrX_data/indexes/chrX_tran -1 
chrX_data/samples/ERR188044_chrX_1.fastq.gz -2 
chrX_data/samples/ERR188044_chrX_2.fastq.gz -S ERR188044_chrX.sam

That is an example

From Nature protocol "Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown"

ADD REPLY
0
Entering edit mode

Are you using the latest hisat version available?

ADD REPLY
0
Entering edit mode

You know, I have change the command:

hisat2-align-s -p 8 --dta -x name_tran -1 input1.fastq.gz -2 input2.fastq.gz -S out.sam

And now I have error like this:

Error: reads file does not look like a FASTQ file
terminate called after throwing an instance of 'int'
Aborted (core dumped)

But they are fastq files ..

Answering youre question: No, I have hisat2-2.0.1 beta

ADD REPLY
0
Entering edit mode

What do you get when you run these commands file input1.fastq.gz and file input2.fastq.gz?

ADD REPLY
0
Entering edit mode
input1.fastq.gz: gzip compressed data, was "input1.fastq", last modified: Fri Jan  3 10:46:39 2020, from Unix

input2.fastq.gz: gzip compressed data, was "input2.fastq", last modified: Fri Jan  3 10:45:58 2020, from Unix
ADD REPLY
0
Entering edit mode

Ok. So these files appear to be in the right format. How much memory do you have available on your machine?

Are you able to successfully run:

hisat2 --dta -x name_tran -1 input1.fastq.gz -2 input2.fastq.gz -S out.sam
ADD REPLY
0
Entering edit mode

Output of youre command:

Command 'hisat2' not found

After mine change:

hisat2-align-s --dta -x name_tran -1 input1.fastq.gz -2 input2.fastq.gz -S out.sam

I get the same error:

Error: reads file does not look like a FASTQ file

Memory is about 80 GB

ADD REPLY
0
Entering edit mode

It is puzzling that you are not able to find the hisat2 wrapper. You are not supposed to run hisat-align directly.

ADD REPLY

Login before adding your answer.

Traffic: 3803 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