We calculate global (unified) wind models of main-sequence, giant, and supergiant O stars from our Galaxy. The models are calculated by solving hydrodynamic, kinetic equilibrium (also known as NLTE) and comoving frame (CMF) radiative transfer equations from the (nearly) hydrostatic photosphere to the supersonic wind. For given stellar parameters, our models predict the photosphere and wind structure and in particular the wind mass-loss rates without any free parameters. Our predicted mass-loss rates are by a factor of 2-5 lower than the commonly used predictions. A possible cause of the difference is abandoning of the Sobolev approximation for the calculation of the radiative force, because our models agree with predictions of CMF NLTE radiative transfer codes. Our predicted mass-loss rates agree nicely with the mass-loss rates derived from observed near-infrared and X-ray line profiles and are slightly lower than mass-loss rates derived from combined UV and Hα diagnostics. The empirical mass-loss rate estimates corrected for clumping may therefore be reconciled with theoretical predictions in such a way that the average ratio between individual massloss rate estimates is not higher than about 1.6. On the other hand, our predictions are by factor of 4.7 lower than pure Hα mass-loss rate estimates and can be reconciled with these values only assuming a microclumping factor of at least eight.