Robust Uncertainty Quantification for Non-Negative Matrix Factorization with Applications to Mutational Signatures Analysis