Abstract:Correct and bias-free interpretation of the deep sequencing data is inevitably dependent on the complete mapping of all mappable reads to the reference sequence, especially for quantitative RNA-seq applications. Seed-based algorithms are generally slow but robust, while Burrows-Wheeler Transform (BWT) based algorithms are fast but less robust. To have both advantages, we developed an algorithm FANSe2 with iterative mapping strategy based on the statistics of real-world sequencing error distribution to substantially accelerate the mapping without compromising the accuracy. Its sensitivity and accuracy are higher than the BWT-based algorithms in the tests using both prokaryotic and eukaryotic sequencing datasets. The gene identification results of FANSe2 is experimentally validated, while the previous algorithms have false positives and false negatives. FANSe2 showed remarkably better consistency to the microarray than most other algorithms in terms of gene expression quantifications. We implemented a scalable and almost maintenance-free parallelization method that can utilize the computational power of multiple office computers, a novel feature not present in any other mainstream algorithm. With three normal office computers, we demonstrated that FANSe2 mapped an RNA-seq dataset generated from an entire Illunima HiSeq 2000 flowcell (8 lanes, 608 M reads) to masked human genome within 4.1 hours with higher sensitivity than Bowtie/Bowtie2. FANSe2 thus provides robust accuracy, full indel sensitivity, fast speed, versatile compatibility and economical computational utilization, making it a useful and practical tool for deep sequencing applications. FANSe2 is freely available at http://bioinformatics.jnu.edu.cn/software/fanse2/.

A Re-Sequencing Tool For High Mismatch-Tolerant Short Read Alignment Based On Burrows-Wheeler Transform

An efficient Burrows-Wheeler transform-based aligner for short read mapping

Fast and accurate long-read alignment with Burrows–Wheeler transform

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

Sensitive Long-Indel-Aware Alignment of Sequencing Reads

Mapping short reads to a genome without using hash look-up table algorithm and Burrows Wheeler Transformation

Fast and Accurate Read Alignment for Resequencing.

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

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

Fast and accurate mapping of Complete Genomics reads

Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

DBG2OLC: Efficient Assembly of Large Genomes Using Long Erroneous Reads of the Third Generation Sequencing Technologies

A Novel Multi-Alignment Pipeline for High-Throughput Sequencing Data.

The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote

ReadsMap: a new tool for high precision mapping of DNAseq and RNAseq read sequences

A fast read alignment method based on seed-and-vote for next generation sequencing

Perm: Efficient Mapping of Short Sequencing Reads with Periodic Full Sensitive Spaced Seeds

ProbAlign: a re-alignment method for long sequencing reads

Read Annotation Pipeline for High-Throughput Sequencing Data.

Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM

FANSe2: a robust and cost-efficient alignment tool for quantitative next-generation sequencing applications.