The origin of the Martian moons is still uncertain, and knowledge about their interior could provide support to some of its leading theories. In preparation for the JAXA Martian Moons eXploration (MMX) mission, we review our current knowledge on the interior of Phobos, and provide synthetic tests showing how the gravity and rotation determination could allow the detection of specific interior-structure properties. The inversion of the geodetic observables for the retrieval of the internal mass distribution of a set of synthetic interior models is performed via non-linear least-squares, where the interior parameterization is based on the level-set method. We additionally provide simple expressions allowing to relate some of these interior models to the geodetic observables of Phobos. The results, based on realistic measurement resolution and noise scenarios, show good retrievals for most of the models at the data resolutions expected from MMX. Specifically, we find the gravity information is realistically sufficient for the detection of mass anomalies below the Stickney crater, as well as large scale heterogeneous regions within plausible rubble-pile structures. Libration helps retrieve the more degenerate models for gravity, such as those with concentric layers or with density varying linearly with depth. The incremental improvement from further adding an hypothetical mean obliquity measurement is marginal. Finally, we apply the level-set inversion and the analytical formulas to estimate possible interior characteristics of the ‘real’ Phobos from the currently-available scarce geodetic observables. The level-set solutions for the real-data inversion generally converge to a higher mass concentration towards the surface in the equatorial region. Markov chain Monte Carlo estimations of parameters relative to a simple 2-layer model or a radial density distribution similarly hint at a lighter region inside of Phobos.