In this article, an elaborated evaluation has been presented for the stratified MHD Jeffrey nanofluid flow towards a stretching surface in the presence of gyrotactic micro-organisms. The analysis of the flow by taking into consideration the viscous dissipation, Brownian motion, Joule heating and thermophoresis and the assumption of thermal, solutal and motile density stratification aspects. The governing nonlinear system of partial differential equations (PDEs), administering the flow, are turned to a group of the nonlinear ordinary differential equations (ODEs) via using an appropriate similarity transformation and later solved numerically by implementing the Keller-box approach with the help of algebraic software Matlab. The impacts of a variety of parameters on the fluid motion, heat, mass, density of the motile micro-organisms, local skin friction, local Nusselt number, local Sherwood number and local density number of the motile micro-organisms are examined with the aid of graphs as well as tables. The results suggest that the motile micro-organisms density is a decreasing function of the bioconvection Lewis number, bioconvection Pecket number and microorganisms concentration difference. The Sherwood number and density rate of motile micro-organisms are bigger in the case of magnetic parameter and Lewis number.