With a massive explosion of datasets across all areas of science and engineering, the central questions are: How do we optimally learn from data through the lens of models? And how do we account for uncertainties in both data and models? These questions can be mathematically framed as Bayesian inverse problems. While powerful and sophisticated approaches have been developed to tackle these problems, such methods are often challenging to implement and there is no available software that easily facilitates the analysis of Bayesian inverse problems. In this talk, we present an extensible FEniCS-based software framework hIPPYlib-MUQ that overcomes these challenges by providing access to state-of-the-art algorithms that offer the capability to solve complex large-scale Bayesian inverse problems across a broad spectrum of scientific and engineering areas.