The focus of this study is to analyzed the 2-dimensional bioconvective flow of Eyring-Powell nanofluid over an exponentially stretching sheet is investigated numerically. Connective boundary conditions for both heat and mass transfer are employed. The governing highly nonlinear partial differential equations are converted into ordinary differential equations by using a similarity transformation. Numerical solutions of the nonlinear ordinary differential equations are found by bvp4c method in MATLAB software. Effects of Eyring-Powell fluid parameter ϵ and δ, Magnetic parameter M, Thermophoresis variable N_t, Lewis number L_b, Peclet number Pe and concentration difference of microorganismsω on velocity, temperature, concentration and motile density profiles are discussed. The nondimensional velocity of the nanofluid is increased as the significance of Eyring-Powell fluid parameter ϵ increases. By increasing the thermophoresis parameter N_t results in increasing profiles of temperature, concentration and motile. The motile profile decreases as the values of Peclet number Pe increases. The motile profile decreases with the rising values of microorganism’s concentration difference ω. Numerical evaluations of the skin friction coefficient, Nusselt numbers and Sherwood numbers are turbulated.