TSS enrichment score serves as an important quality control metric for ATACseq data. I want to write a script for single cell ATACseq data.

From the Encode page:

Transcription Start Site (TSS) Enrichment Score - The TSS enrichment calculation is a signal to noise calculation. The reads around a reference set of TSSs are collected to form an aggregate distribution of reads centered on the TSSs and extending to 1000 bp in either direction (for a total of 2000bp). This …