I'm working with ATAC-Seq data and I need to extract the features in my bed files that are located in transcription start site regions of the mm10 genome. Can anyone recommend a way to do this?
I was thinking of generating a bed file of only tss regions from the mm10 genome and finding the intersect between the TSS file and the sample file. However, I'm not sure how to generate a bed file that only contains TSS sites.
Any help would be appreciated.