Photoacoustic tomography is a hybrid imaging technique that has various applications in biomedicine. In a photoacoustic image reconstruction problem (inverse problem), an initial pressure distribution is reconstructed from measured ultrasound waves which are generated by the photoacoustic effect induced by an optical excitation. In this work, the image reconstruction problem is approached in the framework of Bayesian inversion. The approach is tested with three dimensional numerical simulations. The initial pressure distribution is reconstructed in full-view and limited-view setups. In addition, the reliability of the obtained estimates is assessed. The numerical studies show that accurate estimates of the initial pressure distribution and uncertainty information can be obtained utilizing Bayesian approach