We compute the spectra of open flux tubes formed between a static quark-antiquark pair for various
gauge groups in the large-$N_c$ limit, focusing on different symmetries manifested by the quantum numbers of angular momentum, parity, and charge conjugation. Specifically, we present spectra from $N_c=3$ up to $N_c=6$ and for eight different irreducible representations of the symmetry characterizing the flux tube. In this study, we employed an anisotropic Wilson action, a large number of suitable spacial operators, and solved the generalized eigenvalue problem (GEVP) to obtain a significant number of excitations for different combinations of flux tube quantum numbers. The spectra are compared with the Nambu-Goto string model, revealing novel phenomena such as the presence of massive axions in the flux tube spectrum. We find that the mass of the axion in the open flux tube spectra sector is consistent with the mass obtained in the closed flux tube sector.