We show that if R n is equipped with a certain non-doubling metric and an Orlicz-Sobolev inequality holds for a special family of Young functions Φ , then weak solutions to quasilinear infinitely degenerate elliptic equations of the form − div A x , u ∇ u = ϕ 0 − div A ϕ → 1 are locally bounded. This is obtained by the implementation of a Moser iteration method, what constitutes the first instance of such technique applied to infinite degenerate equations. The results presented here partially extend previously known estimates for solutions of similar equations in which the right hand side does not have a drift term. We also obtain bounds for small negative powers of nonnegative solutions, which will be applied in a subsequent paper to prove continuity of solutions. We also provide examples of geometries in which our abstract theorem is applicable. We consider the family of functions f k , σ x = x ln k 1 x σ , k ∈ N , σ > 0 , − ∞ < x < ∞ , infinitely degenerate at the origin, and show that all weak solutions to − div A x , y , u ∇ u = ϕ x , y − div A ϕ → 1 x , y , A x , y , z ∼ 1 0 0 f k , σ x 2 , with rough data A , ϕ 0 , ϕ → 1 , are locally bounded when k = 1 and 0 < σ < 1 .