Could not locate a HISAT2 index to basename

Could not locate a HISAT2 index to basename

2

First time trying out HISAT2 and I’m having a problem here, even with the pre-made indices for GRCH38.

$ hisat2 -x /share/projects/RNASeq/data/reference/GRCh38/grch38_tran -1 /home/echang/PANCANCER-030817-JE3-35880845/KTP-10-43736695/KTP-10_S3_L001_R1_001.fastq.gz -2 /home/echang/PANCANCER-030817-JE3-35880845/KTP-10-43736695/KTP-10_S3_L001_R2_001.fastq.gz -S tmp.sam

Error follows

Could not locate a HISAT2 index corresponding to basename "/share/projects/RNASeq/data/reference/GRCh38/grch38_tran"
Error: Encountered internal HISAT2 exception (#1)
Command: /home/echang/software/hisat2-2.1.0/hisat2-align-s --wrapper basic-0 -S tmp.sam -1 /tmp/31969.inpipe1 -2 /tmp/31969.inpipe2 /share/projects/RNASeq/data/reference/GRCh38/grch38_tran
(ERR): hisat2-align exited with value 1

The contents of the grch38_tran

$ ls /share/projects/RNASeq/data/reference/GRCh38/grch38_tran
genome_tran.1.ht2  genome_tran.3.ht2  genome_tran.5.ht2  genome_tran.7.ht2  make_grch38_tran.sh
genome_tran.2.ht2  genome_tran.4.ht2  genome_tran.6.ht2  genome_tran.8.ht2

I thought perhaps the genome_tran.*.ht2 needed to change to grch38_tran.*.ht2, but it still does not work after renaming. I’ve since untarred the original file again without change. Was there a step I missed? Where can I find HISTAT2_INDEXES?

-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.


RNA-Seq


Hisat2

• 9.7k views

Read more here: Source link