We compute the spectra of flux tubes formed between a static quark antiquark pair up to a significant number of excitations and for eight symmetries of the flux tubes, up to $\Delta_u$, using pure $SU(3)$ gauge lattice QCD in 3+1 dimensions. To accomplish this goal, we use a large set of appropriate operators, an anisotropic tadpole improved action, smearing techniques, and solve a generalized eigenvalue problem. Moreover, we compare our results with the Nambu-Goto string model to evaluate possible tensions which could be a signal for novel phenomena.