Brownian motion is a central scientific paradigm. Recently, due to increasing efforts and interests towards miniaturization and small-scale physics or biology, the effects of confinement on such a motion have become a key topic of investigation. Essentially, when confined near a wall, a particle moves much slower than in the bulk due to friction at the boundaries. The mobility is therefore locally hindered and space-dependent, which in turn leads to the apparition of so-called multiplicative noises, and associated non-Gaussianities which remain difficult to resolve at all times. Here, we exploit simple, optimized and efficient numerical simulations to address Brownian motion in confinement in a broadrange and quantitative way. To do so, we integrate the overdamped Langevin equation governing the thermal dynamics of a negatively-buoyant single spherical colloid within a viscous fluid confined by two rigid walls, including surface charges. From the produced large set of long random trajectories, we perform a complete statistical analysis and extract all the key quantities, such as the probability distributions in displacements and their main moments. In particular, we propose a novel method to compute high-order cumulants by reducing convergence problems, and employ it to efficiently characterize the inherent non-Gaussianity of the confined process.Read less <