We study numerically the integrable turbulence developing from strongly nonlinear partially coherent waves, in the framework of the focusing one-dimensional nonlinear Schrodinger equation. We find that shortly after the beginning of motion the turbulence enters a state characterized by a very slow evolution of statistics (the quasi-stationary state - QSS), and we concentrate on the detailed examination of the basic statistical functions in this state depending on the shape and the width of the initial spectrum. In particular, we show that the probability density function (PDF) of wavefield intensity is nearly independent of the initial spectrum and is very well approximated by a certain Bessel function representing an integral of the product of two exponential distributions. The PDF corresponds to the value of the second-order moment of intensity equal to 4, indicating enhanced generation of rogue waves. All waves of large amplitude that we have studied are very well approximated - both in space and in time - by the rational breather solutions of either the first (the Peregrine breather), or the second orders.