In this paper we analyze an optimal distributed control problem where the state equations are given by a stationary chemotaxis model coupled with the Navier-Stokes equations. We consider that the movement and the interaction of cells are occurring in a smooth bounded domain of $\mathbb{R}^n,n = 2,3,$ subject to homogeneous boundary conditions. We control the system through a distributed force and a coefficient of chemotactic sensitivity, leading the chemical concentration, the cell density, and the velocity field towards a given target concentration, density and velocity, respectively. In addition to the existence of optimal solution, we derive some optimality conditions.