We calculate the cross section of the gamma-ray production from neutral-current neutrino-oxygen quasi-elastic interaction, $ u+^{16}O rightarrow u +p+^{15}N*$, or $ u+^{16}O rightarrow u+n+^{15}O*$, in which the residual nuclei (15N* or 15O*) lead to the gamma-ray emission with gamma-ray energy >6 MeV at the branching ratio of 41%. Above 200 MeV, this cross section dominates over that of gamma-ray production from the inelastic reaction, $ u+^{16}O-> u+^{16}O*$. In the present calculation, spectral function and the spectroscopic factors of $1p_{1/2}, 1p_{3/2} and 1s_{1/2}$ states are essential. The gamma-ray production is dominated by the deexcitation of $1p_{3/2}$ state of the residual nucleus.