The boracites with general formula M3B7O 13X (M = divalent metal, X = Cl, Br, I), shortly denoted as M-X, are among the first known multiferroic materials. They exhibit a, sequence of transitions from the high temperature paraelectric cubic phase to ferroelectric orthorhombic, monoclinic, trigonal phases, and finally to a monoclinic phase at low temperatures, where both ferroelectric and magnetic orders coexist. The lattice dynamics of boracites has been scarcely studied with the Raman spectroscopy, the main problem with non-cubic phases being the coexistence of twin variants with different crystallographic and polarization orientation. In this work, on the example of Co3B7O13Cl (Co-Cl), Co3B7O13Br (Co-Br) and Ni 3B7Oi13Br (Ni-Br) we demonstrate that using the Raman microscopy imaging one can visualize the twin variants, follow their transformation through the crystallographic transitions, obtain Raman spectra from untwinned domains in exact scattering configurations, determine the Raman mode symmetries, and assign Raman lines to definite atomic motions. The effects of structural transitions and elemental substitution are discussed in close comparison with results of ab initio calculations of the phonon structure of the cubic phase.