Peak calling was performed separately for cells from each of the 16 cell types (excluding glutamatergic neurons; as the caudate is not known to possess excitatory neurons, the presence of this cell type could indicate inadvertent inclusion of another brain region at the time of dissection) using the CallPeaks function in Signac with default parameters. The function uses MACS286 for peak calling.
