Using Matrix Product Operators (MPO) the Schwinger model is simulated in thermal equilibrium. The variational manifold of gauge-invariant MPO is constructed to represent Gibbs states. As a first application, the chiral condensate in thermal equilibrium is computed and agreement with earlier studies is found. Furthermore, as a new application the Schwinger model is probed with a fractional charged static quark-antiquark pair separated infinitely far from each other. A critical temperature beyond which the string tension is exponentially suppressed is found, and is in qualitative agreement with analytical studies in the strong coupling limit. Finally, the CT symmetry breaking is investigated and our results strongly suggest that the symmetry is restored at any nonzero temperature.