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.

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

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

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

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

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

A Short Reads Alignment Algorithm Oriented To Massive Data

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

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

Read Annotation Pipeline for High-Throughput Sequencing Data.

Sensitive Long-Indel-Aware Alignment of Sequencing Reads

Fast and Accurate Read Alignment for Resequencing.

Short Read Alignment Based On Maximal Approximate Match Seeds

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

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

Aligning High Error Rate Reads Using Enhanced Sparse Suffix Array Index

An Efficient Filtration Method Based on Variable-Length Seeds for Sequence Alignment.

Re-alignment of the Unmapped Reads with Base Quality Score

smsMap: mapping single molecule sequencing reads by locating the alignment starting positions

The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads

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

Efficient Distributed Parallel Aligning Reads and Reference Genome with Many Repetitive Subsequences Using Compact de Bruijn Graph