Abstract:BACKGROUND:Many genetic variants have been reported from sequencing projects due to decreasing experimental costs. Compared to the current typical paradigm, read mapping incorporating existing variants can improve the performance of subsequent analysis. This method is supposed to map sequencing reads efficiently to a graphical index with a reference genome and known variation to increase alignment quality and variant calling accuracy. However, storing and indexing various types of variation require costly RAM space.METHODS:Aligning reads to a graph model-based index including the whole set of variants is ultimately an NP-hard problem in theory. Here, we propose a variation-aware read alignment algorithm (VARA), which generates the alignment between read and multiple genomic sequences simultaneously utilizing the schema of the Landau-Vishkin algorithm. VARA dynamically extracts regional variants to construct a pseudo tree-based structure on-the-fly for seed extension without loading the whole genome variation into memory space.RESULTS:We developed the novel high-throughput sequencing read aligner deBGA-VARA by integrating VARA into deBGA. The deBGA-VARA is benchmarked both on simulated reads and the NA12878 sequencing dataset. The experimental results demonstrate that read alignment incorporating genetic variation knowledge can achieve high sensitivity and accuracy.CONCLUSIONS:Due to its efficiency, VARA provides a promising solution for further improvement of variant calling while maintaining small memory footprints. The deBGA-VARA is available at: https://github.com/hitbc/deBGA-VARA.

AliBI: An Alignment-Based Index for Genomic Datasets

Gene Sequence Alignment on a Public Computing Platform

Indexing Finite Language Representation of Population Genotypes

Indexing large genome collections on a PC

DNA sequences alignment method using sparse index on pan-genome graph

BSAlign: A Library for Nucleotide Sequence Alignment

Hybrid Indexes for Repetitive Datasets

Building a pangenome alignment index via recursive prefix-free parsing

Research on Multi-genome Index and Its Improved Sequence Alignment Algorithm

A Short Reads Alignment Algorithm Oriented To Massive Data

SigAlign: an alignment algorithm guided by explicit similarity criteria

A Compressed Self-Index for Genomic Databases

CGAP-align: a high performance DNA short read alignment tool.

Nucleotide String Indexing using Range Matching

Fast and accurate short read alignment with hybrid hash-tree data structure

Fast Read Alignment with Incorporation of Known Genomic Variants.

Integer programming framework for pangenome-based genome inference

A Bidirectional Fuzzy Index and Approximate Search Algorithm for Next Generation Sequencing

MUSIC: A Hybrid Computing Environment for Burrows-Wheeler Alignment for Massive Amount of Short Read Sequence Data

AliView: a fast and lightweight alignment viewer and editor for large datasets