Following our recent work (arxiv:2311.17692), we consider 1/2 BPS supersymmetric circular Wilson loops in four-dimensional N = 2 SU(N) supersymmetric Yang-Mills theories with massless matter in a generic representation of the gauge group and a non-vanishing β-function. On the four-sphere S^4, we can employ supersymmetric localization to map these observable into a matrix model, provided that the one-loop determinants are consistently regularized. After constructing the regularized matrix model for these set-ups, we will demonstrate that the predictions for the Wilson loop at order g^4 match standard perturbative renormalization, based on the evaluation of Feynman diagrams, both on the sphere and also in flat space, despite conformal symmetry being broken at the quantum level. Moreover, I will revisit the difference theory approach, showing that when β-function is non-vanishing this method does not account for the presence of evanescent terms which are activated by the renormalization procedure and contribute to the renormalized observable at order g^6.