Let $Γ_{g}$ be the fundamental group of a closed connected orientable surface of genus $g\geq2$. We develop a new method for integrating over the representation space $\mathbb{X}_{g,n}=\mathrm{Hom}(Γ_{g},S_{n})$ where $S_{n}$ is the symmetric group of permutations of $\{1,\ldots,n\}$. Equivalently, this is the space of all vertex-labeled, $n$-sheeted covering spaces of the the closed surface of genus $g$. Given $φ\in\mathbb{X}_{g,n}$ and $γ\inΓ_{g}$, we let $\mathsf{fix}_γ(φ)$ be the number of fixed points of the permutation $φ(γ)$. The function $\mathsf{fix}_γ$ is a special case of a natural family of functions on $\mathbb{X}_{g,n}$ called Wilson loops. Our new methodology leads to an asymptotic formula, as $n\to\infty$, for the expectation of $\mathsf{fix}_γ$ with respect to the uniform probability measure on $\mathbb{X}_{g,n}$, which is denoted by $\mathbb{E}_{g,n}[\mathsf{fix}_γ]$. We prove that if $γ\inΓ_{g}$ is not the identity, and $q$ is maximal such that $γ$ is a $q$th power in $Γ_{g}$, then \[ \mathbb{E}_{g,n}[\mathsf{fix}_γ]=d(q)+O(n^{-1}) \] as $n\to\infty$, where $d\left(q\right)$ is the number of divisors of $q$. Even the weaker corollary that $\mathbb{E}_{g,n}[\mathsf{fix}_γ]=o(n)$ as $n\to\infty$ is a new result of this paper. We also prove that if $γ$ is not the identity then $\mathbb{E}_{g,n}[\mathsf{fix}_γ]$ can be approximated to any order $O(n^{-M})$ by a polynomial in $n^{-1}$.