A virtual element method for the quasilinear equation −div(κ(u)gradu)=f using general polygonal and polyhedral meshes is presented and analysed. The nonlinear coefficient is evaluated with the piecewise polynomial projection of the virtual element ansatz. Well posedness of the discrete problem and optimal-order a priori error estimates in the H1- and L2-norm are proven. In addition, the convergence of fixed-point iterations for the resulting nonlinear system is established. Numerical tests confirm the optimal convergence properties of the method on general meshes.
