Bioinformatics Advance Access originally published online on January 19, 2007
Bioinformatics 2007 23(5):629-630; doi:10.1093/bioinformatics/btl681
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Compressed suffix tree—a basis for genome-scale sequence analysis
1Department of Computer Science, P.O. Box 68 (Gustaf Hällströmin katu 2b), FI-00014 University of Helsinki, Finland, 2Technische Fakultät, Universität Bielefeld, Germany and 3Indian Institute of Technology, Kanpur, India
*To whom correspondence should be addressed.
| Abstract |
|---|
Summary: Suffix tree is one of the most fundamental data structures in string algorithms and biological sequence analysis. Unfortunately, when it comes to implementing those algorithms and applying them to real genomic sequences, often the main memory size becomes the bottleneck. This is easily explained by the fact that while a DNA sequence of length n from alphabet
= {A, C, G, T } can be stored in n log |
|= 2n bits, its suffix tree occupies O(n log n) bits. In practice, the size difference easily reaches factor 50.
We provide an implementation of the compressed suffix tree very recently proposed by Sadakane (Theory of Computing Systems, in press). The compressed suffix tree occupies space proportional to the text size, i.e. O(n log} |
|) bits, and supports all typical suffix tree operations with at most log n factor slowdown. Our experiments show that, e.g. on a 10 MB DNA sequence, the compressed suffix tree takes 10% of the space of normal suffix tree. Typical operations are slowed down by factor 60.
Availability: The C++ implementation under GNU license is available at http://www.cs.helsinki.fi/group/suds/cst/. An example program implementing a typical pattern discovery task is included. Experimental results in this note correspond to version 0.95.
Contact: vmakinen{at}cs.helsinki.fi
Received on November 3, 2006; revised on January 5, 2007; accepted on January 5, 2007