We then quantified the number of fragments overlapping each peak for each cell using the Signac FeatureMatrix function. We retained all cells that were retained in the analysis performed by the original authors of the BICCN dataset53. We reduced dimensionality using LSI and UMAP as described above using LSI components 2 to 100.