Unveiling non-Abelian statistics of vortex Majorana bound states in iron-based superconductors using fermionic modes