Let D be an open subset of R n with finite measure, and let x0 ∈ D. We introduce the p-Gauss gap of D w.r.t. x0 to measure how far are the averages over D of the harmonic functions u ∈ L p (D) from u(x0). We estimate from below this gap in terms of the ball gap of D w.r.t. x0, i.e., the normalized Lebesgue measure of D\B, being B the biggest ball centered at x0 contained in D.From these stability estimates of the mean value formula for harmonic functions in L p -spaces, we straightforwardly obtain rigidity properties of the Euclidean balls. We also prove a continuity result of the p-Gauss gap in the Sobolev space W 1,p ′, where p ′ is the conjugate exponent of p.