  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 2011, 2012 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.
