Bioinformatics Advance Access originally published online on March 5, 2009
Bioinformatics 2009 25(9):1118-1124; doi:10.1093/bioinformatics/btp131
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
A consistency-based consensus algorithm for de novo and reference-guided sequence assembly of short reads
1International Max Planck Research School for Computational Biology and Scientific Computing, Ihnestr. 63 - 73, 2Algorithmische Bioinformatik, Institut für Informatik, Takustr. 9, 14195 Berlin, Germany and 3J. Craig Venter Institute, 9704 Medical Center Drive, Rockville, MD 20850, USA
*To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: Novel high-throughput sequencing technologies pose new algorithmic challenges in handling massive amounts of short-read, high-coverage data. A robust and versatile consensus tool is of particular interest for such data since a sound multi-read alignment is a prerequisite for variation analyses, accurate genome assemblies and insert sequencing.
Results: A multi-read alignment algorithm for de novo or reference-guided genome assembly is presented. The program identifies segments shared by multiple reads and then aligns these segments using a consistency-enhanced alignment graph. On real de novo sequencing data obtained from the newly established NCBI Short Read Archive, the program performs similarly in quality to other comparable programs. On more challenging simulated datasets for insert sequencing and variation analyses, our program outperforms the other tools.
Availability: The consensus program can be downloaded from http://www.seqan.de/projects/consensus.html. It can be used stand-alone or in conjunction with the Celera Assembler. Both application scenarios as well as the usage of the tool are described in the documentation.
Contact: rausch{at}inf.fu-berlin.de
Associate Editor: Limsoon Wong
Received on November 3, 2008; revised on January 23, 2009; accepted on March 2, 2009