E-MEM: efficient computation of maximal exact matches for very large genomes.
BACKGROUND
Alignment of similar whole genomes is often performed using anchors given by the maximal exact matches (MEMs) between their sequences. In spite of significant amount of research on this problem, the computation of MEMs for large genomes remains a challenging problem. The leading current algorithms employ full text indexes, the sparse suffix array giving the best results. Still, their memory requirements are high, the parallelization is not very efficient, and they cannot handle very large genomes.


RESULTS
We present a new algorithm, efficient computation of MEMs (E-MEM) that does not use full text indexes. Our algorithm uses much less space and is highly amenable to parallelization. It can compute all MEMs of minimum length 100 between the whole human and mouse genomes on a 12 core machine in 10 min and 2 GB of memory; the required memory can be as low as 600 MB. It can run efficiently genomes of any size. Extensive testing and comparison with currently best algorithms is provided.


METHODS
The source code of E-MEM is freely available at: http://www.csd.uwo.ca/∼ilie/E-MEM/ CONTACT: ilie@csd.uwo.ca


BACKGROUND
Supplementary data are available at Bioinformatics online.
Department of Computer Science, University of Western Ontario, London, Ontario, N6A 5B7, Canada.
Nilesh
N
Khiste
Department of Computer Science, University of Western Ontario, London, Ontario, N6A 5B7, Canada.
Lucian
L
Ilie
D000465:Algorithms; D000818:Animals; D016208:Databases, Factual; D016678:Genome; D006801:Humans; D051379:Mice; D011381:Programming Languages; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D012984:Software; D013045:Species Specificity; D049490:Systems Biology / Q000379:methods*; D014908:Triticum / Q000235:genetics