In this work, we study the damage in crystalline molybdenum material samples due to neutron bombardment in a primary knock-on atom (PKA) range of 0.5-10 keV at room temperature. We perform classical molecular dynamics (MD) simulations using a previously derived machine learning (ML) interatomic potential based on the Gaussian approximation potential (GAP) framework. We utilize a recently developed software workflow for fingerprinting and visualizing defects in damaged crystal structures to analyze the Mo samples with respect to the formation of point defects during and after a collision cascade. As a benchmark, we report results for the total number of Frenkel pairs (a self-interstitial atom and a single vacancy) formed and atom displacements as a function of the PKA energy. A comparison to results obtained using an embedded atom method (EAM) potential is presented to discuss the advantages and limits of the MD simulations utilizing ML-based potentials. The formation of Frenkel pairs follows a sublinear scaling law as ξ b where b is a fitting parameter