We present an evaluation of the glueball spectrum for configurations produced with $N_f=1$ dynamical fermions as a function of the $m_{\rm PCAC}$ mass. We obtained masses of states that fall into the irreducible representations of the octahedral group of rotations in combination with the quantum numbers of charge conjugation $C$ and parity $P$. Due to the low signal to noise ratio, practically, we can only extract masses for the irreducible representations $R^{PC}=$ $A_1^{++}$, $E^{++}$, $T_2^{++}$ as well as $A_1^{-+}$. We make use of the Generalized Eigenvalue Problem (GEVP) with an operator basis consisting only of gluonic operators. Throughout this work we are aiming towards the identification of the effects of light dynamical quarks on the glueball spectrum and how this compares to the statistically more precise spectrum of SU(3) pure gauge theory. We used large gauge ensembles which consist of ${\sim {~\cal O}}(10 {\rm K})$ configurations. Our findings demonstrate that the low-lying spectrum of the scalar, tensor as well as pseudo-scalar glueballs receive negligible contributions from the inclusion of $N_f=1$ dynamical fermions.