We construct backreacted geometries dual to the supersymmetric mass deformation of the IKKT matrix model. They are Euclidean type IIB supergravity solutions given in terms of an electrostatic potential, having SO(7)\times SO(3) isometry and 16 supersymmetries. Quantizing the fluxes, we find that the supergravity solutions are in one-to-one correspondence with fuzzy sphere vacua of the matrix model.