The phase transition in the q -state Potts model with homogeneous ferromagnetic couplings is strongly first order for large q, while it is rounded in the presence of quenched disorder. Here we study this phenomenon on different two-dimensional lattices by using the fact that the partition function of the model is dominated by a single diagram of the high-temperature expansion, which is calculated by an efficient combinatorial optimization algorithm. For a given finite sample with discrete randomness the free energy is a piecewise linear function of the temperature, which is rounded after averaging, however, the discontinuity of the internal energy at the transition point (i.e., the latent heat) stays finite even in the thermodynamic limit. For a continuous disorder, instead, the latent heat vanishes. At the phase transition point the dominant diagram percolates and the total magnetic moment is related to the size of the percolating cluster. Its fractal dimension is found d(f) = ( 5 + square root of 5)/4 and it is independent of the type of the lattice and the form of disorder. We argue that the critical behavior is exclusively determined by disorder and the corresponding fixed point is the isotropic version of the so-called infinite randomness fixed point, which is realized in random quantum spin chains. From this mapping we conjecture the values of the critical exponents as beta=2- d(f), beta(s) =1/2, and nu=1.