The aim of the present study is to investigate the effect of uncertainty associated with soil friction angle (phi) and soil unit weight (gamma) on the stability of both unreinforced and reinforced cohesionless soil slopes subjected to strip loading. The magnitude of CoV of phi and gamma are varied to account uncertainties. The location of the footing on the top of the slope is also changed. Stability of both unreinforced and reinforced slopes is presented in terms of factor of safety (FoS). Deterministic FoS values are computed first by using a two-dimensional finite difference software FLAC. To perform probabilistic analyses, FLAC is combined with Monte Carlo simulations. The outcomes of the probabilistic analyses are presented in terms of probability of failure (pF) and reliability index (beta). The value of beta obtained from the present study is compared with the guidelines provided by USACE. It is found out that with the increase in the value of CoV, pF increases and beta decreases. As expected, the failure probability of slope is found to be maximum, when footing is placed on the edge of the unreinforced slope. With the inclusion of a single layer of geotextile in the slope for the same footing position, p(F) reduces drastically, and beta increases significantly. As footing position shifts from the slope edge, pF increases for a particular CoV value of phi and gamma. The effect of uncertainty related to phi is found to be more prominent with compared to the uncertainty related to gamma. The influence of cross-correlation between phi and gamma is also studied. It is found that there is no significant change in the value of pF with the change in the value of cross correlation coefficient. Though the present study is related to a simple slope stability problem, but using the same methodology, probabilistic analyses of complex slopes can also be performed.