Abstract:MOTIVATION:Based on the next generation genome sequencing technologies, a variety of biological applications are developed, while alignment is the first step once the sequencing reads are obtained. In recent years, many software tools have been developed to efficiently and accurately align short reads to the reference genome. However, there are still many reads that can't be mapped to the reference genome, due to the exceeding of allowable mismatches. Moreover, besides the unmapped reads, the reads with low mapping qualities are also excluded from the downstream analysis, such as variance calling. If we can take advantages of the confident segments of these reads, not only can the alignment rates be improved, but also more information will be provided for the downstream analysis.RESULTS:This paper proposes a method, called RAUR (Re-align the Unmapped Reads), to re-align the reads that can not be mapped by alignment tools. Firstly, it takes advantages of the base quality scores (reported by the sequencer) to figure out the most confident and informative segments of the unmapped reads by controlling the number of possible mismatches in the alignment. Then, combined with an alignment tool, RAUR re-align these segments of the reads. We run RAUR on both simulated data and real data with different read lengths. The results show that many reads which fail to be aligned by the most popular alignment tools (BWA and Bowtie2) can be correctly re-aligned by RAUR, with a similar Precision. Even compared with the BWA-MEM and the local mode of Bowtie2, which perform local alignment for long reads to improve the alignment rate, RAUR also shows advantages on the Alignment rate and Precision in some cases. Therefore, the trimming strategy used in RAUR is useful to improve the Alignment rate of alignment tools for the next-generation genome sequencing.AVAILABILITY:All source code are available at http://netlab.csu.edu.cn/bioinformatics/RAUR.html.

Aligning High Error Rate Reads Using Enhanced Sparse Suffix Array Index

ResSeq: Enhancing Short-Read Sequencing Alignment by Rescuing Error-Containing Reads

A Short Reads Alignment Algorithm Oriented To Massive Data

Short Read Alignment Based On Maximal Approximate Match Seeds

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

Sensitive Long-Indel-Aware Alignment of Sequencing Reads

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

Accelerating spliced alignment of long RNA sequencing reads using parallel maximal exact match retrieval

Efficient Frequency-Based De Novo Short-Read Clustering for Error Trimming in Next-Generation Sequencing

Fast and Accurate Read Alignment for Resequencing.

An algorithm for rapid noncoding RNA sequence-structure alignment

Rhat: Fast Alignment of Noisy Long Reads with Regional Hashing.

Re-alignment of the Unmapped Reads with Base Quality Score

Maximum Match Subsequence Alignment Algorithm Finely Grained (MMSAA FG)

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

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

Umap: Use Unique Sequence for Alignment of Short Sequence Reads and SNP Detection

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

DEEP-LONG: A Fast and Accurate Aligner for Long RNA-Seq

Accel-Align: a fast sequence mapper and aligner based on the seed-embed-extend method

BSAlign: A Library for Nucleotide Sequence Alignment