Many classes of extended scalar-tensor theories predict that dynamical instabilities can take place at high energies, leading to the formation of scalarized neutron stars. Depending on the theory parameters, stars in a scalarized state can form a solution-space branch that shares a lot of similarities with the so-called mass twins in general relativity appearing for equations of state containing first-order phase transitions. Members of this scalarized branch have a lower maximum mass and central energy density compared to Einstein ones. In such cases, a scalarized star could potentially overaccrete beyond the critical mass limit, thus triggering a gravitational phase transition where the star sheds its scalar hair and migrates over to its nonscalarized counterpart. Such an event resembles, but is distinct from, a nuclear or thermodynamic phase transition. We dynamically track a gravitational transition by first constructing hydrostatic, scalarized equilibria for realistic equations of state, and then allowing additional material to fall onto the stellar surface. The resulting bursts of monopolar radiation are dispersively stretched to form a quasicontinuous signal that persists for decades, carrying strains of order ≳10^{-22} (kpc/L)^{3/2} Hz^{-1/2} at frequencies of ≲300 Hz, detectable with the existing interferometer network out to distances of L≲10 kpc, and out to a few hundred kpc with the inclusion of the Einstein Telescope. Electromagnetic signatures of such events, involving gamma-ray and neutrino bursts, are also considered.