Both cold and flame jets find numerous applications in different fields, ranging from domestic applications to aerospace and space technology. Indeed, the applications of isothermal and non-isothermal jets in the flame heating industry fascinated the researchers to gain an in-depth understanding. Nevertheless, these benefits are not standalone, rather, they are associated with major disadvantages such as improper jet mixing and flame instabilities that require careful remedies. In the present investigation, three-inline jets, having methane jet at the center and two oxygen jets at the periphery, are studied computationally in a three-dimensional domain, with and without considering the effects of combustion. To study the mixing characteristics of cold jets, the radial velocity distributions at different streamwise locations have been analyzed at the jet inlet velocity of 27 m/s. However, for oxygen and methane flame jets, inlet velocities are varied as 27 m/s and 54 m/s. Moreover, to investigate the effects of temperature variation on mixing characteristics at a particular jet velocity, the inlet temperatures of reactants are varied as 300 K, 500 K, and 700 K, at the jet inlet velocity of 27 m/s. Combustion is found to have a profound impact on the mixing characteristics. At the inlet temperature of 300 K, a higher centerline velocity decay is observed for non-reactive jets as compared to reactive flame jets. Moreover, the turbulent kinetic energy distribution is relatively higher in the case of non-reactive jets, which is the direct evidence of an augmented mixing. As is obvious, the turbulent kinetic energy at the jet shear layer is maximum due to the formation of large-scale coherent eddies. The decay in centerline velocity is found to be increasing with an increase of inlet temperature. Additionally, with an increase of jet velocity, the radial velocity profiles become more asymmetrical, consequently yielding an unstable flame.