Probabilistic neural operators for functional uncertainty quantification