From what I understand, I don’t think mk_angndx returns indices for water molecules (at least I couldn’t find a way to get it to do so, if someone knows how then please correct me), if you were to simulate a protein in your system, then the command should work for the atoms of that protein.
You should be able to find the data on the angles of the water molecules from the .itp file of the water model you used in the folder containing your forcefield.