We develop an exact map of all states and operators from 2D lattices of spins-1/2 into lattices of fermions and bosons with mutual semionic statistical interaction that goes beyond previous dualities of Z2 lattice gauge theories because it does not rely on imposing local conservation laws and captures the motion of "charges" and "fluxes" on equal footing. This map allows to explicitly compute the Berry phases for the transport of fluxes in symmetry enriched topologically ordered states that can be either chiral, non-chiral, abelian or non-abelian, and whose numerical complexity reduces to diagonalizing free-fermion Hamiltonians. Among other results, we establish numerically the conditions under which the Majorana-carrying flux excitation in Ising-Topologically-Ordered states enriched by translations acquires 0 or π phase when moving around a single plaquette.