The twinning of crystals, on the macroscopic scale, sometimes manifests as the presence of reentrant crystal angles.(S6) When the twinning interfaces are not parallel to each other then they develop as cyclic twins, which can resemble stars. However, the presence of morphological irregularities is in itself not a sufficient indicator. While concave crystal faces do not always indicate twinning, it is certainly a warning sign.(S5) On the atomic level, twinning, as with other crystal imperfections, might manifest as certain pathologies such as split reflections and deviation from the expected intensity statistics of diffraction data.(S7) However, some forms of twinning can be particularly difficult to determine and not evident until full data sets are analyzed. In particular, merohedral twinning is a type of twinning where the lattices of two or more distinct domains coincide exactly.(S8) As the real space lattices coincide, the reciprocal lattices of the domains will exactly overlap and the resulting diffraction pattern will appear normal. This type of twinning can occur in space groups of high symmetry, e.g., tetragonal, trigonal, hexagonal or cubic, where more than one Laue symmetry is possible. In such cases, the additional rotational symmetry operator of the twinning exactly superimposes parts of the diffraction pattern, causing the appearance of higher symmetry, i.e., the twinning adds symmetry. The true Laue symmetry (point symmetry) of the crystal is lower than that of the lattice.(S9) The hexagonal space group P622 can be a result of the merohedral twinning of the lower symmetry trigonal space group P321. Therefore, it is imperative for our analysis that this possibility be checked and definitively confirmed or denied. As explained above, since merohedral twinning will not be evident from the diffraction pattern, the intensity distribution of the complete data set needs to be examined. One such widely used analysis involves the identification of twins from intensity statistics. The values of the <|E2-1|> can be plotted and are known to differ for centrosymmetric vs. noncentrosymmetric structures. The theoretical statistical value for centric (S10) structures is 0.938, and that of acentric structures 0.736. Twinning results in a diffraction pattern of higher symmetry and yields a hypo-noncentrosymmetric intensity distribution.(S11) The cumulative distribution function N(z) gives the fraction of reflections whose relative intensities are less than z, where z is the normalized reflection intensity obtained by dividing each individual measurement by the average value for its resolution shell.(S9) When the plot of the cumulative intensity distribution N(Z) for acentric and centric data is plotted, the shape of the plot is a diagnostic. The plot for an un-twinned crystal appears exponential due to the existence of a small number of very weak or very strong reflections, while for a twinned crystal, the plot appears sigmoidal. The sigmoidal shape arises (S12) because a portion of the non-equivalent reflections in a twinned crystal will overlap exactly; some of the very strong reflection