Global weighted L estimates are obtained for the gradient of solutions to nonlinear elliptic Dirichlet boundary value problems over a bounded nonsmooth domain. Morrey and Hölder regularity of solutions are also established, as a consequence. These results generalize various existing estimates for nonlinear equations. The nonlinearities are of at most linear growth and assumed to have a uniform small mean oscillation. The boundary of the domain, on the other hand, may exhibit roughness but assumed to be sufficiently flat in the sense of Reifenberg. Our approach uses maximal function estimates and Vitali covering lemma, and also known regularity results of solutions to nonlinear homogeneous equations.