Clusters of correlated TF presence. Clusters (k = 12) of TF co-binding were identified using non-negative matrix factorization (NMF) applied to TF presence matrix at the (a) TSSs and (b) enhancers separately. The TF membership (H matrix) for top three clusters important for prediction is shown, next to their SHAP values