To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.