The equation of motion for the Skyrme model with a pion mass term is studied in the framework of the Painlevé analysis, and information is obtained about the singularity structure of its solutions. As in the massless case, singularities exist, consisting of a first-order pole term superposed to a logarithmic branch point. For the solitonic solutions, the singularities form an infinite sequence of points located on the negative real axis of the squared radial variable, accumulating at the origin. Based on this property and on the asymptotic behavior of the solitonic solutions, an approximate representation is built for the baryonic soliton, which is extremely accurate.