One often resorts to a non-minimal coupling of the electromagnetic field in order to generate magnetic fields during inflation. The coupling is expected to depend on a scalar field, possibly the same as the one driving inflation. At the level of three-point functions, such a coupling leads to a non-trivial cross-correlation between the perturbation in the scalar field and the magnetic field. This cross-correlation has been evaluated analytically earlier for the case of non-helical electromagnetic fields. In this work, we numerically compute the cross-correlation for helical magnetic fields. Non-Gaussianities are often generated as modes leave the Hubble radius. The helical electromagnetic modes evolve strongly (when compared to the non-helical case) around Hubble exit and one type of polarization is strongly amplified immediately after Hubble exit. We find that helicity considerably boosts the amplitude of the dimensionless non-Gaussianity parameter that characterizes the amplitude and shape of the cross-correlation between the perturbations in the scalar field and the magnetic field. We discuss the implications of the enhancement in the non-Gaussianity parameter due to parity violation. © 2018 IOP Publishing Ltd and Sissa Medialab.