Annotating large genomes with exact word matches.
We have developed a tool for rapidly determining the number of exact matches of any word within large, internally repetitive genomes or sets of genomes. Thus we can readily annotate any sequence, including the entire human genome, with the counts of its constituent words. We create a Burrows-Wheeler transform of the genome, which together with auxiliary data structures facilitating counting, can reside in about one gigabyte of RAM. Our original interest was motivated by oligonucleotide probe design, and we describe a general protocol for defining unique hybridization probes. But our method also has applications for the analysis of genome structure and assembly. We demonstrate the identification of chromosome-specific repeats, and outline a general procedure for finding undiscovered repeats. We also illustrate the changing contents of the human genome assemblies by comparing the annotations built from different genome freezes.
Healy
John
J
Cold Spring Harbor Laboratory, Cold Spring Harbor, New York 11724, USA. healy@cshl.edu
Thomas
Elizabeth E
EE
Schwartz
Jacob T
JT
Wigler
Michael
M
D000465:Algorithms; D002888:Chromosomes, Human, Pair 19 / Q000737:chemistry / Q000235:genetics; D019295:Computational Biology / Q000379:methods* / Q000706:statistics & numerical data; D020451:Contig Mapping / Q000379:methods / Q000639:trends; D015342:DNA Probes / Q000235:genetics; D018076:DNA, Complementary / Q000235:genetics; D030541:Databases, Genetic / Q000706:statistics & numerical data; D015894:Genome, Human; D006801:Humans
D015342:DNA Probes; D018076:DNA, Complementary