We develop a general framework to (numerically) study adiabatic braiding of quasiholes in fractional quantum Hall systems. Specifically, we investigate the Moore-Read (MR) state at nu=1/2 filling factor, a known candidate for non-Abelian statistics, which appears to actually occur in nature. The non-Abelian statistics of MR quasiholes is demonstrated explicitly for the first time, confirming the results predicted by conformal field theories.