The bioconvection flow of an incompressible micropolar fluid containing microorganisms between two infinite stretchable parallel plates is considered. A mathematical model, with a fully coupled nonlinear system of equations describing the total mass, momentum, thermal energy, mass diffusion, and microorganisms is presented. The governing equations are reduced to a set of nonlinear ordinary differential equations with the help of suitable transformations. The resulting nonlinear ordinary differential equations are linearized using successive linearization method, and the resulting system of linear equations is solved using the Chebyshev collocation method. The detailed analysis illustrating the influences of various physical parameters, such as the micropolar coupling number, squeezing parameter, the bioconvection Schmidt number, Prandtl numbers, Lewis number, and bioconvection Peclet number on the velocity, microrotation, temperature, concentration and motile microorganism distributions, skin friction coefficient, Nusselt number, Sherwood number, and density number of motile microorganism, is examined. The influence of the squeezing parameter is to increase the dimensionless velocities and temperature and to decrease the local Nusselt number and local Sherwood number. The density number of motile microorganism is decreasing with squeezing parameter, bioconvection Lewis number, bioconvection Peclet number, and bioconvection Schmidt number.