KHN biquad: why does LTSPICE plot a wrong phase plot for bandpass filter?
What you see is called phase wrapping, and it's the result of calculating the phase with the four-quadrant version of atan(), atan2(). Its domain is from \$-\pi\$ to \$\pi\$, or, in degrees, from \$-180^\text{o}\$ to \$180^\text{o}\$.
In LTspice (or most other SPICEs), there is an option to unwrap the phase, which displays it continuously, without the jumps: if you're using version XVII then right-click, else left-click on the right side of the waveform window, on the Y-axis. You will see this little dialog come up (ignore the readings, it's a quick test run):

The Unravel Branch Wrap is checked, by default. Unchecking it will result in the phase being displayed as you see in Mathematica. If you want for Mathematica to also display it unwrapped, then you'll have to do a bit of calculation involving the derivative of the phase, searching for the jumps, and adding or subtracting \$2\pi\$, as needed. Or maybe there's a builtin function.