Summary: ``Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme ($M1$) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes {\tt POSSIS} and {\tt Skynet}, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively.''