We present first-principle calculations on the vertical ionization potentials (IPs), electron affinities (EAs), and singlet excitation energies on an aromatic-molecule test set (benzene, thiophene, 1,2,5-thiadiazole, naphthalene, benzothiazole, and tetrathiafulvalene) within the GW and BetheSalpeter equation (BSE) formalisms. Our computational framework, which employs a real-space basis for ground-state and a transition-space basis for excited-state calculations, is well-suited for high-accuracy calculations on molecules, as we show by comparing against G0W0 calculations within a plane-wave-basis formalism. We then generalize our framework to test variants of the GW approximation that include a local-density approximation (LDA)-derived vertex function (ΓLDA) and quasiparticle-self-consistent (QS) iterations. We find that ΓLDA and quasiparticle self-consistency shift IPs and EAs by roughly the same magnitude, but with opposite sign for IPs and same sign for EAs. G0W0 and QSGW ΓLDA are more accurate for IPs, while G0W0ΓLDA and QSGW are best for EAs. For optical excitations, we find that perturbative GW -BSE underestimates the singlet excitation energy, while self-consistent GW -BSE results in good agreement with previous best-estimate values for both valence and Rydberg excitations. Finally, our work suggests that a hybrid approach, where G0W0 energies are used for occupied orbitals and G0W0ΓLDA for unoccupied orbitals, also yields optical excitation energies in good agreement with experiment but at a smaller computational cost.