Question: ATAC-seq macs2 peak splitting in sliding windows
gravatar for nancy
10 months ago by
nancy0 wrote:

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?

atac-seq • 337 views
ADD COMMENTlink modified 9 months ago by danvoronov10 • written 10 months ago by nancy0

What exactly is the problem?

Edit 26/9/19: Cross-posted =>

ADD REPLYlink modified 9 months ago • written 9 months ago by ATpoint36k
gravatar for danvoronov
9 months ago by
SZN, Naples, Italy
danvoronov10 wrote:

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

ADD COMMENTlink written 9 months ago by danvoronov10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1060 users visited in the last hour