General Information
| Description | An application to short-read mapping toward accurate SNP/indel detection |
| Homepage | http://last.cbrc.jp |
Publications
Probabilistic alignments with quality scores: an application to short-read mapping toward accurate SNP/indel detection
Motivation:
Recent studies have revealed the importance of considering quality scores of reads generated by next-generation sequence (NGS) platforms in various downstream analyses. It is also known that probabilistic alignments based on marginal probabilities (e.g. aligned-column and/or gap probabilities) provide more accurate alignment than conventional maximum score-based alignment. There exists, however, no study about probabilistic alignment that considers quality scores explicitly, although the method is expected to be useful in SNP/indel callers and bisulfite mapping, because accurate estimation of aligned columns or gaps is important in those analyses.
Results:
In this study, we propose methods of probabilistic alignment that consider quality scores of (one of) the sequences as well as a usual score matrix. The method is based on posterior decoding techniques in which various marginal probabilities are computed from a probabilistic model of alignments with quality scores, and can arbitrarily trade-off sensitivity and positive predictive value (PPV) of prediction (aligned columns and gaps). The method is directly applicable to read mapping (alignment) toward accurate detection of SNPs and indels. Several computational experiments indicated that probabilistic alignments can estimate aligned columns and gaps accurately, compared with other mapping algorithms e.g. SHRiMP2, Stampy, BWA and Novoalign. The study also suggested that our approach yields favorable precision for SNP/indel calling.
Availability:
The method described in this article is implemented in LAST, which is freely available from: http://last.cbrc.jp.
Contact:
mhamada@k.u-tokyo.ac.jp
Supplementary Information:
Supplementary data are available at Bioinformatics online.
Related Tools
StatAlign 2.0
Published January 2013
Combining statistical alignment with RNA secondary structure prediction
PicXAA
Published May 2010
Greedy probabilistic construction of maximum expected accuracy alignment of multiple sequences
FAIR Evaluation
Sorry, you must be signed in to evaluate this object.
