The energy is converged to less then 1 microHartree, and the OH bond length is converged to 0.0002 angstroms. You can also see that the final geometry parameters obtained with the default optimization criteria depend somewhat on the initial starting geometry. VeryTight all but eliminates these differences. VeryTight alone is not necessarily enough to converge the geometry to the point where the low frequencies are as close to zero as you would like. The results are in Table 4 .
More accurate grids are necessary for a truly converged optimization.
The frequencies are sorted by increasing absolute value, so that it's easier to distinguish rotational modes from vibrational modes. Given the straight forward convergence seen with Hartree-Fock theory, this might not seem to make sense. However, it does make sense if you recall that DFT is done using a numerical integration on a grid of points. The accuracy of the default grid is not high enough for computing low frequency modes very precisely. The solution is to use a more numerically accurate grid.