We provide the full asymptotic description of the quasi-normal modes (resonances) in any strip of fixed width for Dirac fields in slowly rotating Kerr-Newman-de Sitter black holes. The resonances split in a way similar to the Zeeman effect. The method is based on the extension to Dirac operators of techniques applied by Dyatlov in [17], [18] to the (uncharged) Kerr-de Sitter black holes. We show that the mass of the Dirac field does not have effect on the two leading terms in the expansions of resonances.