The effect of Mach number on transonic flow past a circular cylinder is investigated numerically for the free-stream Mach number M ∞ from 0.85 to 0.98 and the Reynolds number 2×105 based on the diameter of the cylinder. The work provides an insight into several salient features, including unsteady and quasi-steady flow state, formation of local supersonic zone, and evolution of turbulent shear layer. Results show that there exist two flow states dependent of a critical Mach number M cr around 0.9. One is an unsteady flow state characterized by moving shock waves interacting with the turbulent flow in the near region of the cylinder for M ∞<M cr, and the other is a quasi-steady flow state with nearly stationary shock waves formed in the near wake for M ∞>M cr, suppressing vortex shedding from the cylinder. Some flow behaviors in the unsteady and quasi-steady flow states are revealed. From time evolution of flow structures, local supersonic zones are identified in the wake and generated by two typical processes, i.e. reverse flow behind the cylinder and shed vortices in the near wake. The convective Mach number M c of turbulent shear layers shed from the cylinder is identified nearly as M c<1 in the unsteady flow regime and M c>1 in the quasi-steady flow regime, resulting in different evolutions of the shear layers.