I do not know if you can obtain directly rates from BOLISG, even if I think it is possible. However it is not very complex to do by yourself.
When you have obtained the distribution, you download the cross sections from LXCAT and it is sufficient to calculate the integral (summation over the energy grid) of the quantity EEDF(E)*v(E)*sigma(E)dE. EEDF should be in eV^-1.
The collision rate coefficients ki(cm3 s−1) may be defined by where σ is the cross-section (cm2) of the collision i and v is the electronic velocity (cm s−1) related with electron temperature. The collision rate coefficients depend on the pressure through the electronic energy distribution function.
Moreover, the presence of inelastic collisions between electrons and neutrals can result
in either ionization or excitation of the neutral particle.