ATAC-seq macs2 peak splitting in sliding windows
I used macs2 to call peaks for atac-seq data. now my goal is to split the peaks into 50 bp windows with 25 bp steps and then calculate the Tn5 integration frequency in each window.
How should I proceed with that?
• 1,971 views
•
link
1 answer
Hello, you can use bedtools makewindows -b <your narrowpeak files from macs2> -w 50 -s 25 > <output file> this will give you the sliding windows you want.
The Tn5 integration frequency is unclear, because you should have the fragments generated by the transposase "cutting" the genome, so I think you need 2 transposase reactions to generate one fragment. Could you elaborate on what you need?
Best, Dan
• 0 views
•
link
Log in to answer this question.
What exactly is the problem?
Edit 26/9/19: Cross-posted => https://bioinformatics.stackexchange.com/questions/9365/atac-seq-macs2-peak-splitting-in-sliding-windows