Revealing long-term effects of contaminants on the genetic structure of organisms inhabiting polluted environments should encompass analyses at the population, molecular, and cellular level. Following this concept, we studied the genetic constitution of zebra mussel populations from a polluted (Dp) and reference sites (Cl) at the river Drava, Croatia, and applied microsatellite and DNA damage analyses (Comet assay, micronucleus test (MNT)). Additionally, mussels from both populations were exposed to polluted wastewater in the laboratory for three days, and DNA damage was analyzed to evaluate acclimatization and genetic adaptation of the investigated populations to the polluted environment. The two populations differed in their genetic constitution. Microsatellite analysis suggested that Dp had undergone a genetic bottleneck. Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that Dp had an increased percentage of micronuclei in hemocytes in comparison to Cl. The laboratory experiment revealed that Dp had a lower percentage of tail DNA and a higher percentage of micronuclei than Cl. These differences between populations were possibly caused by an overall decreased fitness of Dp due to genetic drift and by an enhanced DNA repair mechanism due to acclimatization to pollution in the source habitat.