Arbitrary-order mixed methods for heterogeneous anisotropic diffusion on general meshes