Genome indexing with hisat2 and samtools
2
2
Entering edit mode
3.2 years ago
MK ▴ 40

Hello, I found this code in Biostar handbook (RNA-Seq by Example) in section about genome indexing.

# The reference genome.
IDX=refs/genome.fa

# Build the genome index.
hisat2-build $IDX $IDX

# Index the reference genome with samtools.
samtools faidx $IDX

I was wondering about samtools indexing step. I've never seen anyone using it before. Isn't hisat2-build not enough to properly index genome?

samtools index Hisat2 RNA-seq genome • 1.5k views
ADD COMMENT
2
Entering edit mode
3.2 years ago
Rafael Soler ★ 1.3k

Yes, the HISAT2 indexing step is enough :)

ADD COMMENT
0
Entering edit mode
3.2 years ago
GenoMax 147k

samtools index fasta is required for querying fasta files, extracting regions from fasta file.

ADD COMMENT

Login before adding your answer.

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