Compute chromVAR deviations for groups of peaks. The goal of this function is similar to that of AddModuleScore except that it is designed for single-cell chromatin data. The chromVAR deviations for each group of peaks will be added to the object metadata.
features: A named list of features to include in each module. The name of each element in the list will be used to name the modules computed, which will be stored in the object metadata.
genome: A BSgenome object
assay: Name of assay to use. If NULL, use the default assay.