Image reconstruction and uncertainty quantification in photoacoustic tomography