The paper is concerned to find the distribution of the chemically reactant solute in the MHD flow of an electrically conducting viscous incompressible fluid over a stretching surface. The first order chemical reaction and the variable solute distribution along the surface are taken into consideration. The governing partial differential equations along with appropriate boundary conditions for flow field and reactive solute are transformed into a set of non-linear self-similar ordinary differential equations by using scaling group of transformations. An exact analytic solution is obtained for the velocity field. Using this velocity field, we obtain numerical solution for the reactant concentration field. It reveals from the study that the values of concentration profile enhances with the increase of the magnetic field and decreases with increase of Schmidt number as well as the reaction rate parameter. Most importantly, when the solute distribution along the surface increases then the concentration profile decreases
