We investigate the supersymmetric Sachdev-Ye-Kitaev (SYK) model, $N$ Majorana fermions with infinite range interactions in $0+1$ dimensions. We have found that, close to the ground state $E approx 0$, discrete symmetries alter qualitatively the spectral properties with respect to the non-supersymmetric SYK model. The average spectral density at finite $N$, which we compute analytically and numerically, grows exponentially with $N$ for $E approx 0$. However the chiral condensate, which is normalized with respect the total number of eigenvalues, vanishes in the thermodynamic limit. Slightly above $E approx 0$, the spectral density grows exponential with the energy. Deep in the quantum regime, corresponding to the first $O(N)$ eigenvalues, the average spectral density is universal and well described by random matrix ensembles with chiral and superconducting discrete symmetries. The dynamics for $E approx 0$ is investigated by level fluctuations. Also in this case we find excellent agreement with the prediction of chiral and superconducting random matrix ensembles for eigenvalues separations smaller than the Thouless energy, which seems to scale linearly with $N$. Deviations beyond the Thouless energy, which describes how ergodicity is approached, are universality characterized by a quadratic growth of the number variance. In the time domain, we have found analytically that the spectral form factor $g(t)$, obtained from the connected two-level correlation function of the unfolded spectrum, decays as $1/t^2$ for times shorter but comparable to the Thouless time with $g(0)$ related to the coefficient of the quadratic growth of the number variance. Our results provide further support that quantum black holes are ergodic and therefore can be classified by random matrix theory.