We consider the statistical linear inverse problem of recovering the unknown initial heat state from noisy interior measurements over an inhomogeneous domain of the solution to the heat equation at a fixed time instant. We employ nonparametric Bayesian procedures with Gaussian series priors defined on the Dirichlet-Laplacian eigenbasis, yielding convenient conjugate posterior distributions with explicit expressions for posterior inference. We review recent theoretical results that provide asymptotic performance guarantees (in the large sample size limit) for the resulting posterior-based point estimation and uncertainty quantification. We further provide an implementation of the approach, and illustrate it via a numerical simulation study.