Stress boundary conditions for the lattice Boltzmann equation that are consistent to Burnett order are proposed and imposed using a moment-based method. The accuracy of the method with complicated spatially-dependent boundary conditions for stress and velocity is investigated using the regularized lid-driven cavity flow. The complete set of boundary conditions, which involve gradients evaluated at the boundaries, are implemented locally. A recently-derived collision operator with modified equilibria and velocity-dependent collision rates to reduce the defect in Galilean invariance is also investigated. Numerical results are in excellent agreement with existing benchmark data and exhibit second-order convergence. The lattice Boltzmann stress field is studied and shown to depart significantly from the Newtonian viscous stress when the ratio of Mach to Reynolds numbers is not negligibly small.