We introduce a numerical method for Brownian dynamics with position dependent diffusion tensor which is second order accurate for sampling the invariant measure while requiring only one force evaluation per timestep. Analysis of the sampling bias is performed using the algebraic framework of exotic aromatic Butcher-series. Numerical experiments confirm the theoretical order of convergence and illustrate the efficiency of the new method.