Rotational-vibrational spectroscopy of water in solid noble gas matrices has been studied for many decades. Despite that, discrepancies persist in the literature about the assignment of specific bands. We tackle the involved rotational-vibrational spectrum of the water isotopologues H216O, HD16O, and D216O with an unprecedented combination of experimental high-resolution matrix isolation infrared (MI-IR) spectroscopy and computational anharmonic vibrational spectroscopy by vibrational configuration interaction (VCI) on high-level ab initio potential energy surfaces. With VCI, the average deviation to gas-phase experiments is reduced from >100 to ≈1 cm-1 when compared to harmonic vibrational spectra. Discrepancies between MI-IR and VCI spectra are identified as matrix effects rather than missing anharmonicity in the theoretical approach. Matrix effects are small in Ne (≈1.5 cm-1) and a bit larger in Ar (≈10 cm-1). Controversial assignments in Ne MI-IR spectra are resolved, for example, concerning the ν3 triad in HDO. We identify new transitions, for example, the ν2 101 ← 110 transition in D2O and H2O or the ν3 000 ← 101 transition in D2O, and reassign bands, for example, the band at 3718.9 cm-1 that is newly assigned as the 110 ← 111 transition. The identification and solution of discrepancies for a well-studied benchmark system such as water prove the importance of an iterative and one-hand combination of theory and experiment in the field of high-resolution infrared spectroscopy of single molecules. As the computational costs involved in the VCI approach are reasonably low, such combined experimental/theoretical studies can be extended to molecules larger than triatomics.