MATLAB: $\tt freqz$ vs $\tt bode$
In fact both plots are the same. 

Your freqz plot uses a linear discrete frequency, $\omega$ axis scaled between $0$ and $1$ with units of $\pi$ rad per sample. For 200 Hz sampling frequency, a 60 Hz notch would be shown at $w_c = 0.6 \pi$ rads, which is put at 0.6.

The bode plot, however, uses a logarithmic scale for the continuous-time radian frequency $\Omega$. Given that you have a $f_c=60$ Hz of notch frequency, its radian value is $\Omega_c = 2\pi \times  60 \approx 377$ radians per second.

So both plots are showing the notch at the correct place.