The true model is
y = ßx + ε
Stata calculates
ß = Σxy / Σxx
Substitution shows
ß = ß + Σxε / Σxx
Assume E(xε) = 0
E(ß) = ß
V(ß) = V(Σxε) / (Σxx)2 = σ2 / Σxx
 
V(Σxε) = V(x1ε1 + ... + xnεn)
= x12V(ε1) + ... + xn2V(εn)
... + cov(xiεi,xjεj) + ... for all ij
= x12V(ε1) + ... + xn2V(εn)
= σ2 Σxx
The true model is
y = ßx + δz + ε
Stata calculates
ß = Σxy / Σxx
Substitution shows
ß = ß + Σxz / Σxx + Σxε / Σxx
E(ß) = ß + Σxz / Σxx