MTfit is a Python module for Bayesian moment tensor source inversion of earthquake seismic data using polarities, amplitudes or amplitude ratios. It can solve for double couple or full moment tensor solutions, taking into account uncertainties in polarities, take-off angles of the rays from the source to the receiver, and amplitudes. It provides an easily accessible and extendable approach to earthquake source inversion which is particularly useful for local and regional events.