We study a random matrix model which interpolates between the singular values of the Gaussian unitary ensemble (GUE) and of the chiral Gaussian unitary ensemble (chGUE). This symmetry crossover is analogous to the one realized by the Hermitian Wilson Dirac operator in lattice QCD, but our model preserves chiral symmetry of chGUE exactly unlike the Hermitian Wilson Dirac operator. This difference has a crucial impact on the statistics of near-zero eigenvalues, though both singular value statistics build a Pfaffian point process. The model in the present work is motivated by the Dirac operator of 3d staggered fermions, 3d QCD at finite isospin chemical potential, and 4d QCD at high temperature. We calculate the spectral statistics at finite matrix dimension. For this purpose, we derive the joint probability density of the singular values, the skew-orthogonal polynomials and the kernels for the k-point correlation functions. The skew-orthogonal polynomials are constructed by the method of mixing bi-orthogonal and skew-orthogonal polynomials, which is an alternative approach to Mehta's one. We compare our results with Monte Carlo simulations and study the limits to chGUE and GUE. As a side product, we also calculate a new type of a unitary group integral.PACS numbers: 02.10. Yn,05.50.+q,11.15.Ex,11.15.Ha,: 15B52, 33C45 ‡ The authors of [20,35] did not study the pion condensate ud + du . As such, they had no source term j and thus had a decoupling of the two quarks.