The aim of this paper is to obtain a Bayesian estimator of stress-strength reliability based on generalized order statistics for Pareto distribution. The dependence of the Pareto distribution support on the parameter complicates the calculations. Hence, in literature, one of the parameters is assumed to be known. In this paper, for the first time, two parameters of Pareto distribution are considered unknown. In computing the Bayesian confidence interval for reliability based on generalized order statistics, the posterior distribution has a complex form that cannot be sampled by conventional methods. To solve this problem, we propose an acceptance-rejection algorithm to generate a sample of the posterior distribution. We also propose a particular case of this model and obtain the classical and Bayesian estimators for this particular case. In this case, to obtain the Bayesian estimator of stress-strength reliability, we propose a variable change method. Then, these confidence intervals are compared by simulation. Finally, a practical example of this study is provided.