Fast and accurate long-read alignment with Burrows–Wheeler transform
Tóm tắt
Motivation: Many programs for aligning short sequencing reads to a reference genome have been developed in the last 2 years. Most of them are very efficient for short reads but inefficient or not applicable for reads >200 bp because the algorithms are heavily and specifically tuned for short queries with low sequencing error rate. However, some sequencing platforms already produce longer reads and others are expected to become available soon. For longer reads, hashing-based software such as BLAT and SSAHA2 remain the only choices. Nonetheless, these methods are substantially slower than short-read aligners in terms of aligned bases per unit time.
Results: We designed and implemented a new algorithm, Burrows-Wheeler Aligner's Smith-Waterman Alignment (BWA-SW), to align long sequences up to 1 Mb against a large sequence database (e.g. the human genome) with a few gigabytes of memory. The algorithm is as accurate as SSAHA2, more accurate than BLAT, and is several to tens of times faster than both.
Availability: http://bio-bwa.sourceforge.net
Contact: [email protected]
Từ khóa
Tài liệu tham khảo
Altschul, 1997, Gapped BLAST and PSI-BLAST: a new generation of protein database search programs, Nucleic Acids Res., 25, 3389, 10.1093/nar/25.17.3389
Blumer, 1985, The smallest automaton recognizing the subwords of a text, Theor. Comput. Sci., 40, 31, 10.1016/0304-3975(85)90157-4
Burrows, 1994, A block-sorting lossless data compression algorithm, Technical report 124
Eid, 2009, Real-time DNA sequencing from single polymerase molecules, Science, 323, 133, 10.1126/science.1162986
Ferragina, 2000, Opportunistic data structures with applications, Proceedings of the 41st Symposium on Foundations of Computer Science (FOCS 2000), 390
Jiang, 2008, SeqMap: mapping massive amount of oligonucleotides to the genome, Bioinformatics, 24, 2395, 10.1093/bioinformatics/btn429
Kent, 2002, BLAT–the BLAST-like alignment tool, Genome Res., 12, 656
Lam, 2008, Compressed indexing and local alignment of DNA, Bioinformatics, 24, 791, 10.1093/bioinformatics/btn032
Langmead, 2009, Ultrafast and memory-efficient alignment of short DNA sequences to the human genome, Genome Biol., 10, R25, 10.1186/gb-2009-10-3-r25
Li, 2009, Fast and accurate short read alignment with Burrows-Wheeler transform, Bioinformatics, 25, 1754, 10.1093/bioinformatics/btp324
Li, 2008, Mapping short DNA sequencing reads and calling variants using mapping quality scores, Genome Res., 18, 1851, 10.1101/gr.078212.108
Li, 2009, The Sequence Alignment/Map format and SAMtools, Bioinformatics, 25, 2078, 10.1093/bioinformatics/btp352
Li, 2008, SOAP: short oligonucleotide alignment program, Bioinformatics, 24, 713, 10.1093/bioinformatics/btn025
Ma, 2002, PatternHunter: faster and more sensitive homology search, Bioinformatics, 18, 440, 10.1093/bioinformatics/18.3.440
Meek, 2003, OASIS: an online and accurate technique for local-alignment searches on biological sequences, Proceedings of 29th International Conference on Very Large Data Bases (VLDB 2003), 910
Morgulis, 2008, Database indexing for production megablast searches, Bioinformatics, 24, 1757, 10.1093/bioinformatics/btn322
Ning, 2001, SSAHA: a fast search method for large DNA databases, Genome Res., 11, 1725, 10.1101/gr.194201
Pearson, 1988, Improved tools for biological sequence comparison, Proc. Natl Acad. Sci. USA, 85, 2444, 10.1073/pnas.85.8.2444
Rumble, 2009, SHRiMP: accurate mapping of short color-space reads, PLoS Comput. Biol., 5, e1000386, 10.1371/journal.pcbi.1000386
Smith, 2008, Using quality scores and longer reads improves accuracy of Solexa read mapping, BMC Bioinformatics, 9, 128, 10.1186/1471-2105-9-128
Weese, 2009, RazerS–fast read mapping with sensitivity control, Genome Res., 19, 1646, 10.1101/gr.088823.108