HISAT2
HISAT2 is a fast and sensitive alignment program for mapping next-generation sequencing reads (both DNA and RNA) to a population of human genomes (as well as against a single reference genome). Based on an extension of BWT for graphs [Sirén et al. 2014], the HISAT2 developers designed and implemented a graph FM index (GFM). In addition to using one global GFM index that represents a population of human genomes, HISAT2 uses a large set of small GFM indexes that collectively cover the whole genome (each index representing a genomic region of 56 Kbp, with 55,000 indexes needed to cover the human population). These small indexes (called local indexes), combined with several alignment strategies, enable rapid and accurate alignment of sequencing reads.
Documentation
- HISAT2 web site
- Kim D, Langmead B and Salzberg SL.HISAT: a fast spliced aligner with low memory requirements. Nature Methods2015
Usage on Bridges-2
To see what versions of HISAT2 are available and if there is more than one, which is the default, along with some help, type
module spider hisat2
To use HISAT2, include a command like this in your batch script or interactive session to load the HISAT2 module: (note ‘module load’ is case-sensitive):
module load hisat2