Resting state functional magnetic resonance imaging (rs-fMRI) dynamic functional network connectivity (dFNC) analysis has illuminated brain network interactions across many neuropsychiatric disorders. A common analysis approach involves using hard clustering methods to identify transitory states of brain activity, and in response to this, other methods have been developed to quantify the importance of specific dFNC interactions to identified states. Some of these methods involve perturbing individual features and examining the number of samples that switch states. However, only a minority of samples switch states. As such, these methods actually identify the importance of dFNC features to the clustering of a subset of samples rather than the overall clustering. In this study, we present a novel approach that more capably identifies the importance of each feature to the overall clustering. Our approach uses fuzzy clustering to output probabilities of each sample belonging to states and then measures their Kullback-Leibler divergence after perturbation. We show the viability of our approach in the context of schizophrenia (SZ) default mode network analysis, identifying significant differences in state dynamics between individuals with SZ and healthy controls. We further compare our approach with an existing approach, showing that it captures the effects of perturbation upon most samples. We also find that interactions between the posterior cingulate cortex (PCC) and the anterior cingulate cortex and the PCC and precuneus are important across methods. We expect that our novel explainable clustering approach will enable further progress in rs-fMRI analysis and to other clustering applications.