HISAT2 alignment issue
1
0
Entering edit mode
18 days ago
eleven11 • 0

Hello! I am trying to do a pair end read alignment. When I run the following code --

~/hisat2/hisat2  hisat2 -x ~/hisat2/grch38 -1 ~/hisat2/astra_r1.fastq.gz -2 ~/hisat2/astra_r2.fastq.gz -S ~hisat2/ouput.sam

I get that

(ERR): "/Users/simonlab/hisat2/grch38" does not exist Exiting now ...

I know I tried many different iterations of defining the index sequence but keep getting the same error. I know I have to define the path for hisat2 to run. Any help appreciated! Thanks!

RNAseq Alignment HISAT2 • 235 views
ADD COMMENT
0
Entering edit mode

Many alignment tools need the full path to their alignment index files with prefix. What is the output to the following commands:

file ~/hisat2/grch38
ls -lRh ~/hisat2/grch38*
ADD REPLY
0
Entering edit mode

Hello! output to file command is

/Users/simonlab/hisat2/grch38: directory

Output to ls command is

total 8797648
-rw-r--r--@ 1 lab  staff   940M Mar 12  2016 genome.1.ht2
-rw-r--r--@ 1 lab  staff   702M Mar 12  2016 genome.2.ht2
-rw-r--r--@ 1 lab  staff    11K Mar 12  2016 genome.3.ht2
-rw-r--r--@ 1 lab  staff   702M Mar 12  2016 genome.4.ht2
-rw-r--r--@ 1 lab  staff   1.2G Mar 12  2016 genome.5.ht2
-rw-r--r--@ 1 lab  staff   715M Mar 12  2016 genome.6.ht2
-rw-r--r--@ 1 lab  staff     8B Mar 12  2016 genome.7.ht2
-rw-r--r--@ 1 lab  staff     8B Mar 12  2016 genome.8.ht2
-rwxr-xr-x@ 1 lab  staff   1.5K Mar 17  2016 make_grch38.sh
(base) simonlab-imac:~ simonlab$ 

I have also tried running the orginal tar.gz file and get the same error. Thank you!

ADD REPLY
0
Entering edit mode
18 days ago
Ram 43k

You need -x ~/hisat2/grch38/genome. The manual says:

Main arguments

-x <hisat2-idx> The basename of the index for the reference genome. The basename is the name of any of the index files up to but not including the final .1.ht2 / etc. hisat2 looks for the specified index first in the current directory, then in the directory specified in the HISAT2_INDEXES environment variable.

ADD COMMENT

Login before adding your answer.

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