In this paper, we present the global $L^q$ estimates for the weak solution of the steady p-Stokes equations describing the motion of the shear-thinning flow under the nonslip boundary condition. Our interest is a non-smooth domain whose boundary can go beyond the Lipschitz category, and the coefficients belong to BMO(Bounded Mean Oscillation) space with sufficiently small BMO semi-norm.