Abstract: The unquenched spectrum of the Wilson Dirac operator in the microscopic (epsilon) regime including order a^2-corrections is discussed. We explain how to obtain exact analytical expressions for the microscopic spectral density from the supersymmetric method and Wilson chiral perturbation theory. The results are derived for individual topological sectors, where topology at nonzero lattice spacing is defined by the index of the Wilson Dirac operator. The microscopic spectral density shows in detail the transition to the Aoki phase. Finally, we discuss constraints on the coefficients of Wilson chiral perturbation theory.