GMAP: a genomic mapping and alignment program for mRNA and EST sequences.
BACKGROUND
We introduce GMAP, a standalone program for mapping and aligning cDNA sequences to a genome. The program maps and aligns a single sequence with minimal startup time and memory requirements, and provides fast batch processing of large sequence sets. The program generates accurate gene structures, even in the presence of substantial polymorphisms and sequence errors, without using probabilistic splice site models. Methodology underlying the program includes a minimal sampling strategy for genomic mapping, oligomer chaining for approximate alignment, sandwich DP for splice site detection, and microexon identification with statistical significance testing.


RESULTS
On a set of human messenger RNAs with random mutations at a 1 and 3% rate, GMAP identified all splice sites accurately in over 99.3% of the sequences, which was one-tenth the error rate of existing programs. On a large set of human expressed sequence tags, GMAP provided higher-quality alignments more often than blat did. On a set of Arabidopsis cDNAs, GMAP performed comparably with GeneSeqer. In these experiments, GMAP demonstrated a several-fold increase in speed over existing programs.


BACKGROUND
Source code for gmap and associated programs is available at http://www.gene.com/share/gmap


BACKGROUND
http://www.gene.com/share/gmap.
Wu
Thomas D
TD
Department of Bioinformatics Genentech, Inc., South San Francisco, CA 94080, USA. twu@gene.com
Watanabe
Colin K
CK
D000465:Algorithms; D002874:Chromosome Mapping / Q000379:methods*; D020224:Expressed Sequence Tags; D015894:Genome, Human; D006801:Humans; D022821:RNA Splice Sites / Q000235:genetics; D012333:RNA, Messenger / Q000235:genetics*; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D017423:Sequence Analysis, RNA / Q000379:methods*; D012689:Sequence Homology, Nucleic Acid; D012984:Software
D022821:RNA Splice Sites; D012333:RNA, Messenger