Subsea production pipelines in deep water oil and gas fields are susceptible to a phenomenon called lateral buckling which is a major pipeline design concern. Accurate estimation of lateral buckle formation and the additional stresses generated are of particular importance for High-Temperature High-Pressure (HTHP) pipelines resting on soft clay and is the subject of pipe-soil interaction studies (PSI). Variability in pipe embedment and soil resistances constitute the largest uncertainty in PSI analyses. Best practice methods for PSI studies constitute multi-step strategies that involve finding solutions to nonlinear and implicit equations. In the present study, the impact of soil strength variability on embedment and subsequently peak (Hpeak/V) and residual (Hres/V) lateral friction factors are evaluated. An R code, PSI-Lateral, is developed and verified against published case studies. It is used to calculate the pipeline embedment and horizontal resistances mobilized in a parametric study on three pipe diameters (0.6, 0.8 m and 1.0 m) with three different wall thicknesses. For each pipe diameter-weight combination, embedment and lateral friction factors (peak and residual) are calculated for a series of mean shear strengths ranging from 2 kPa to 8 kPa and COV ranging from 5% to 37.5%. The high estimate and low estimate shear strengths, suHE and suLE, correspond to values that are +/- two standard deviations from the mean. Hence, the range of strengths encompassed between these two extremes represents approximately 95% of the shear strength distribution. The results of the parametric runs conducted are presented in the form of graphical charts where high and low embedment and lateral friction factors can be read-off for any selected pair of mean shear strength and COV. Similar trends are observed in the charts for all pipelines evaluated. It is found that the range of normalized embedment, z/D, increases with increasing COV and decreases with increasing mean shear strength for all pipelines tested. A similar trend is observed for the residual horizontal resistance. The expected range of peak horizontal resistance Hpeak/V shows the same trend but only for soils with mean strengths less than 4 kPa. For soils with higher mean strengths, the range of Hpeak/V is found to be mostly insensitive to variations in COV and further increases in the mean strength.