The analysis of organic mercury (Hg) species in polluted soils is a necessary tool to assess the environmental risk(s) of mercury in contaminated legacy sites. The artificial formation of monomethylmercury (MeHg) during soil extraction and/or analysis is a well-known limitation and is especially relevant in highly polluted areas where MeHg/Hg ratios are notoriously low. Although this has been known for almost 30 years, the thorough characterisation of artificial formation rates is rarely a part of the method development in scientific literature. Here we present the application of two separate procedures (inorganic Hg (iHg) spiking and double-spike isotope dilution analyses (DSIDA)) to determine and correct for artificial Hg methylation in MeHg-selective acid-leaching/organic solvent extraction procedure. Subsequently, we combined corrected MeHg and ethylmercury (EtHg) measurements with PCR amplification of hgcA genes to distinguish between naturally formed MeHg from primary deposited MeHg in soils from a legacy site in a Swiss mountain valley. We found the DSIDA procedure incompatible with the organomercury selective extraction method due to the quantitative removal of iHg. Methylation factors from iHg spiking were in the range of (0.0075 ± 0.0001%) and were consistent across soils and sediment matrices. Further, we suggest that MeHg was deposited and not formed in-situ in two out of three studied locations. Our line of evidence consists of 1) the concomitant detection of EtHg, 2) the elevated MeHg concentrations (up to 4.84 μg kg^-1), and 3) the absence of hgcA genes at these locations. The combination of Hg speciation and methylation gene (hgcA) abundance analyses are tools suited to assess Hg pollution pathways at Hg legacy sites.