This paper develops a least-squares approach to the solution of the optimal control problem for the Navier-Stokes equations. We recast the optimality system as a first-order system by introducing velocity-flux variables and associated curl and trace equations. We show that a least-squares principle based on $L^2$ norms applied to this system yields optimal discretization error estimates in the $H^1$ norm in each variable.