Bioremediation is becoming an increasingly popular approach for the remediation of sites contaminated with the explosive hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). Multiple lines of evidence are often needed to assess the success of such approaches, with molecular studies frequently providing important information on the abundance of key biodegrading species. Towards this goal, the current study utilized shotgun sequencing to determine the abundance and diversity of functional genes (xenA, xenB, xplA, diaA, pnrB, nfsI) and species previously associated with RDX biodegradation in groundwater before and after biostimulation at an RDX-contaminated Navy Site. For this, DNA was extracted from four and seven groundwater wells pre- and post-biostimulation, respectively. From a set of 65 previously identified RDX degraders, 31 were found within the groundwater samples, with the most abundant species being Variovorax sp. JS1663, Pseudomonas fluorescens, Pseudomonas putida, and Stenotrophomonas maltophilia. Further, 9 RDX-degrading species significantly (p<0.05) increased in abundance following biostimulation. Both the sequencing data and qPCR indicated that xenA and xenB exhibited the highest relative abundance among the six genes. Several genes (diaA, nsfI, xenA, and pnrB) exhibited higher relative abundance values in some wells following biostimulation. The study provides a comprehensive approach for assessing biomarkers during RDX bioremediation and provides evidence that biostimulation generated a positive impact on a set of key species and genes.