The generalization of the Bernstein polynomials based on Polya distribution is considered in the present article. Here, we introduce a mixed summation-integral type operators having Polya and Bernstein basis functions in summation and integration respectively. We establish some direct results which include an asymptotic formula, local and global approximation results for these operators in terms of modulus of continuity.