How to find the total phase of this transfer function
There is no error in the phase response given by Wolfram Alpha. If you look at your transfer function, you can see the combination of three components:

a right-half-plane zero (RHPZ) lies in the numerator
a pole at the origin for \$s=0\$
two poles located at higher frequency

The first component will shift the phase down to -90°. This is because the zero is located in the right half-plane, indicated by the minus sign in the numerator. Should you have a plus sign instead, the zero would be in the left half-plane (a LHPZ) and the phase would lead up to 90°.
Then the pole at the origin (a 1st-order integrator) will permanently lag the phase by 90°. So together with the RHPZ, you already have -90-90=-180°. Then, the two poles are located in the left half-plane and contribute another 180° phase lag, totaling the phase response to -180-180=-360° or 0°.
A quick Mathcad sheet shows the response of your formula while going through different arrangements. As expected, at frequencies beyond 100 Hz, the phase is 0° as predicted: