Pleistocene events played an important role in the differentiation of North American vertebrate populations. Michigan, in particular, and the Great Lakes region, in general, were greatly influenced by the last glaciation. While several hypotheses regarding the recolonization of this region have been advanced, none have been strongly supported. We generated 148 complete ND2 mitochondrial DNA (mtDNA) sequences from common gartersnake (Thamnophis sirtalis) populations throughout the Great Lakes region to evaluate phylogeographic patterns and population structure and to determine whether the distribution of haplotypic variants is related to the post-Pleistocene retreat of the Wisconsinan glacier. The common gartersnake was utilized, as it is believed to have been one of the primary vertebrate invaders of the Great Lakes region following the most recent period of glacial retreat and because it has been a model species for a variety of evolutionary, ecological, behavioral, and physiological studies. Several genetically distinct evolutionary lineages were supported by both genealogical and molecular population genetic analyses, although to different degrees. The geographic distribution of the majority of these lineages is interpreted as reflecting post-glacial recolonization dynamics during the late Pleistocene. These findings generally support previous hypotheses of range expansion in this region.