We study the multi-boundary entanglement structure of the states prepared in (1+1) and (2+1) dimensional Chern-Simons theory with finite discrete gauge group G. The states in (1+1)-d are associated with Riemann surfaces of genus g with multiple S1 boundaries and we use replica trick to compute the entanglement entropy for such states. In (2+1)-d, we focus on the states associated with torus link complements which live in the tensor product of Hilbert spaces associated with multiple T2. We present a quantitative analysis of the entanglement structure for both abelian and non-abelian groups. For all the states considered in this work, we find that the entanglement entropy for direct product of groups is the sum of entropy for individual groups, i.e. EE(G1× G2) = EE(G1) + EE(G2). Moreover, the reduced density matrix obtained by tracing out a subset of the total Hilbert space has a positive semidefinite partial transpose on any bi-partition of the remaining Hilbert space. © 2020, The Author(s).