Two sections of the control region and the genes coding for NADH dehydrogenase subunits 5 and 6 (ND-5/6) of mitochondrial DNA (mtDNA) were amplified from Phoxinus eos with the polymerase chain reaction. Both sections of the control region were sequenced directly while the ND-5/6 fragment was sequenced in from each end only. Additionally, the entire ND-5/6 fragment was examined for sequence variation using RFLP analysis. No sequence variation was detected in the control region among 70 individuals sampled from 18 populations across three Ontario regions (Spanish River, Madawaska R. and Cataraqui R.). To examine ND-5/6 variation, a total of 75 individuals were sampled from five populations representing two of the three regions (Madawaska River and Cataraqui R.). Six haplotypes were detected by direct sequencing and four by RFLP analysis. Estimates of population subdivision from RFLP data, sequence analysis, and the two data sets combined for the ND-5/6 fragment, suggest that gene flow is restricted within and between regions. However, estimates of sequence divergence for both sequence and RFLP analysis of this fragment suggested that populations were either founded by already differentiated populations or that populations were founded by a single stock and divergence between regions occurred prior to isolation of populations within regions. These estimates of population structure are much greater than those obtained from allozyme analysis. Additionally, high levels of heterozygosity in nuclear DNA, but low mtDNA diversity suggests that populations have experienced reductions in population size sufficient to reduce only mtDNA variation. Random lineage extinction and limited time for the accumulation of new mutations are likely responsible for low levels of mtDNA variation in ND-5/6 and the control region, while functional constraints may limit variation more than expected in the control region in dace and other fishes.