How matrix dislocations, i.e. basal, prismatic and pyramidal, interact with {101¯2}101¯1¯ twin boundaries in hexagonal close-packed metals has been discussed extensively in the literature. However, so far no systematic investigation has been reported. In this work, we performed atomistic simulations to study interaction between matrix dislocations in pure Mg with a {101¯2} twin boundary. Our results show that for the basal and the prismatic slip, when the Burgers vector is parallel to the zone axis of the twins, a matrix basal dislocation can be transmuted to a twin prismatic dislocation and vice versa. However, when the Burgers vector of the matrix dislocation is non-parallel to the zone axis, no transmutation occurs and the dislocation is absorbed by the twin boundary which acts as a dislocation sink. For a matrix pyramidal dislocation, the dislocation is absorbed by the twin boundary and no transmutation occurs either. It appears that if the product dislocation is a real slip system, transmutation may occur during twin-slip interaction. Otherwise the matrix dislocation will be shredded by atomic shuffling and then absorbed by the twin boundary. If the core structure of the product dislocation is complex, transmutation may not occur as well and dislocation absorption will occur. Lattice correspondence in deformation twinning was applied in explaining the interaction mechanisms. Our results can be well correlated to macroscopic experimental observations which show twin-slip interaction only contributes negligibly to work hardening in deformation of hcp metals.