1. FBB: a fast Bayesian-bound tool to calibrate RNA-seq aligners
- Author
-
Jeff Hasty, Irene Rodriguez-Lujan, and Ramon Huerta
- Subjects
0301 basic medicine ,Statistics and Probability ,Computer science ,Bayesian probability ,RNA-Seq ,Saccharomyces cerevisiae ,Biochemistry ,Genome ,03 medical and health sciences ,0302 clinical medicine ,Escherichia coli ,Molecular Biology ,Base Sequence ,Sequence Analysis, RNA ,Bayes Theorem ,FBB ,Sequence Analysis, DNA ,Base (topology) ,Computer Science Applications ,Computational Mathematics ,030104 developmental biology ,Computational Theory and Mathematics ,Calibration ,RNA ,Algorithm ,Algorithms ,Software ,030217 neurology & neurosurgery - Abstract
Motivation Despite RNA-seq reads provide quality scores that represent the probability of calling a correct base, these values are not probabilistically integrated in most alignment algorithms. Based on the quality scores of the reads, we propose to calculate a lower bound of the probability of alignment of any fast alignment algorithm that generates SAM files. This bound is called Fast Bayesian Bound (FBB) and serves as a canonical reference to compare alignment results across different algorithms. This Bayesian Bound intends to provide additional support to the current state-of-the-art aligners, not to replace them. Results We propose a feasible Bayesian bound that uses quality scores of the reads to align them to a genome of reference. Two theorems are provided to efficiently calculate the Bayesian bound that under some conditions becomes the equality. The algorithm reads the SAM files generated by the alignment algorithms using multiple command option values. The program options are mapped into the FBB reference values, and all the aligners can be compared respect to the same accuracy values provided by the FBB. Stranded paired read RNA-seq data was used for evaluation purposes. The errors of the alignments can be calculated based on the information contained in the distance between the pairs given by Theorem 2, and the alignments to the incorrect strand. Most of the algorithms (Bowtie, Bowtie 2, SHRiMP2, Soap 2, Novoalign) provide similar results with subtle variations. Availability and Implementation Current version of the FBB software is provided at https://bitbucket.org/irenerodriguez/fbb. Supplementary information Supplementary data are available at Bioinformatics online.
- Published
- 2016
- Full Text
- View/download PDF