We present a high-order accurate discontinuous Galerkin method for evolving the spherically reduced Baumgarte-Shapiro-Shibata-Nakamura (BSSN) system expressed in terms of second-order spatial operators. Our multidomain method achieves global spectral accuracy and longtime stability on short computational domains. We discuss in detail both our scheme for the BSSN system and its implementation. After a theoretical and computational verification of the proposed scheme, we conclude with a brief discussion of issues likely to arise when one considers the full BSSN system.