Genome indexing with hisat2 and samtools
2
1
Entering edit mode
13 days ago
MK ▴ 30

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 • 142 views
ADD COMMENT
2
Entering edit mode
12 days ago

Yes, the HISAT2 indexing step is enough :)

ADD COMMENT
0
Entering edit mode
12 days ago
GenoMax 107k

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

ADD COMMENT

Login before adding your answer.

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