Neural Operator induced Gaussian Process framework for probabilistic solution of parametric partial differential equations