Many sources of uncertainty exist when emissions are modelled for a gas turbine combustion system. They originate from uncertain inputs, boundary conditions, calibration, or lack of sufficient fidelity in the model. In this paper, a non-intrusive polynomial chaos expansion (NIPCE) method is coupled with a chemical reactor network (CRN) model using Python to rigorously quantify uncertainties of NOx emission in a premixed burner. The first objective of the uncertainty quantification (UQ) in this study is development of a global sensitivity analysis method based on NIPCE to capture aleatory uncertainty due to the variation of operating conditions and input parameters. The second objective is uncertainty analysis of Arrhenius parameters in the chemical kinetic mechanism to study the epistemic uncertainty in the modelling of NO emission. A two-reactor CRN consisting of a perfectly stirred reactor (PSR) and a plug flow reactor (PFR) is constructed in this study using Cantera to model NOx for natural gas at the relevant operating conditions for a benchmark premixed burner. UQ is performed through the use of a number of packages in Python. The results of uncertainty and sensitivity analysis using NIPCE based on point collocation method (PCM) are then compared with the results of advanced Monte Carlo simulation (MCS). Surrogate models are also developed based on the NIPCE approach and compared with the forward model in Cantera to predict NOx emissions. The results show the capability of NIPCE approach for UQ using a limited number of evaluations to develop a UQ-enabled emission prediction tool for gas turbine combustion systems.