Krylov and Safonov estimates for degenerate quasilinear elliptic PDEs