Wakes of bluﬀ bodies in a stratiﬁed environment are common in oceanic and atmospheric ﬂows. Some examples are marine swimmers, underwater submersibles and ﬂow over mountains and islands. Direct numerical simulations of ﬂow past a sphere in a stratiﬁed ﬂuid at a sub-critical Reynolds number (Re) of 3,700 and for a range of Froude numbers, Fr = U/ND ∈ [0.025,∞] are performed. The conservation equations are solved in a cylindrical coordinate system and an immersed boundary method is employed to represent the sphere. The prime objective of this investigation is to understand the statistical response of the near, intermediate and far wake of a sphere at sub-critical Re under the inﬂuence of buoyancy. It is observed that buoyancy leads to the inhibition of vertical motion resulting in faster decay of r.m.s. velocity in the vertical direction as compared to the horizontal r.m.s. velocity, collapse of the wake, propagation of internal gravity waves and the organization of the primarily horizontal ﬂow into coherent vortical structures. Unprecedented with respect to previous studies, the time averaged turbulent kinetic energy budget is closed for the unstratiﬁed and stratiﬁed cases. A novel ﬁnding of this research is the regeneration of turbulent ﬂuctuations in the near wake when the stratiﬁcation increases beyond a critical level (Fr decreases beyond a critical value) which is in contrast to the previous results at lower Re that suggest monotone suppression of turbulence with increasing stratiﬁcation. Vorticity evolution, energy spectra and the turbulence energy equation explain turbulence regeneration. Another objective of this study is to quantify the distinction between the body and turbulence generated internal waves, in terms of the amplitude, frequency, potential energy distribution and propagation angles. With a decrease in Fr, the body generation mechanism become stronger and waves exhibit upstream propagation.