Reaction-diffusion systems are of considerable importance in many areas of physical sciences. For many reasons; an external manipulation of the dynamics of these processes is desirable. Here we show numerically how spatiotemporal behavior like pattern formation and wave propagation in a two component nonlinear reaction-diffusion model of bacterial chemotaxis can be externally controlled. We formulate the control goal as an objective functional and apply numerical optimization for the solution of the resulting control problem.