We numerically study the Hamiltonian lattice formulation of the two-flavor Schwinger model using matrix product states. Keeping the mass of the first flavor at a fixed positive value, we tune the mass of the second flavor through a range of negative values, thus exploring a regime where conventional Monte Carlo methods suffer from the sign problem and may run into instabilities due to zero modes. Our results indicate a phase transition at the point where the absolute value of the second flavor mass approaches the first flavor mass. The phase transition is accompanied by the formation of a fermion condensate, a steep drop of the average electric field, and a peak in the bipartite entanglement entropy. Our data hints at a second order transition, which is the 1+1D analog of the CP-violating Dashen phase transition in QCD.