The Coulomb excitation experiment to study electromagnetic properties of the heaviest stable Mo isotope, Mo, was performed using a 76 MeV S beam from the Warsaw cyclotron U-200P. Magnitudes and relative signs of 26 , , , and matrix elements coupling nine low-lying states in Mo were determined using the least-squares code gosia. Diagonal matrix elements (related to the spectroscopic quadrupole moments) of the 2, 2, and 2 states as well as the 4 state were extracted. The resulting set of reduced matrix elements was complete and precise enough to obtain, using the quadrupole sum rules approach, quadrupole deformation parameters of Mo in its two lowest 0 states: ground and excited. The overall deformation of the 0 and 0 states in Mo is of similar magnitude, in both cases larger compared to what was found for the neighboring isotopes Mo and Mo. At the same time, the asymetry parameters obtained for both states strongly differ, indicating a triaxial shape of the Mo nucleus in the ground state and a prolate shape in the excited 0 state. Low-energy quadrupole excitations of the Mo nucleus were studied in the frame of the general quadrupole collective Bohr Hamiltonian model (GBH). The potential energy and inertial functions were calculated using the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) method starting from two possible variants of the Skyrme effective interaction: SIII and Sly4. The overall quadrupole deformation parameters resulting from the GBH calculations with the SLy4 variant of the Skyrme interaction are slightly closer to the experimentally obtained values than those obtained using SIII.