<p>It's impossible to find nucleosome free and such regions without a BAM file. Your only options are:</p>
<ol>
<li>Filter the BAM files to contain only the size-range of interest and rerun the counting and statistics.</li>
<li>Find or write a counting program that can track fragment sizes.</li>
</ol>
<p>BTW, those size ranges could use some adjustments, you're throwing away all fragments between 100 and 180 bases using their ranges.</p>

CC BY-SA 4.0
2020-12-31T07:07:02.233Z
Thank you devon . "BTW, those size ranges could use some adjustments, you're throwing away all fragments between 100 and 180 bases using their ranges." how to do it add some region?
kcm
CC BY-SA 4.0
2020-12-31T10:52:29.557Z
Looks at a size distribution plot and adjust accordingly.
Devon Ryan
CC BY-SA 4.0
2020-12-30T07:33:37.917Z
2020-12-30T07:33:37.917Z
Devon Ryan
<p>I m following this <a href="https://rockefelleruniversity.github.io/RU_ATAC_Workshop.html" rel="nofollow noreferrer">tutorial</a> for my ATAC seq analysis. Here in of the section they show profile of open chromatin region plots.</p>
<pre><code>**Plotting insert sizes with Greenleaf open, mono- and di-nucleosome profiles
This looks very similar to the image from the Greenleaf paper.

We can now annotate our nucleosome free (&lt; 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) as in the Greenleaf study.**
</code></pre>
<p>Now Since i have done the differential accessibility analysis. So from that information how do i find out region which are &quot;nucleosome free (&lt; 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) &quot;.</p>
<p>My objective</p>
<ul>
<li>I would like to see this &quot;nucleosome free (&lt; 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) &quot; from my differential accessible region which was done from readcount and then I ran deseq2.</li>
<li>My goal is not to plot the way they have done rather find those regions mentioned above and see it between different conditions.</li>
</ul>
<p>Now based on the tutorial of that section where they do for each bam file.</p>
<p>Do i need to run the analysis on each bam file or each condition or Can I find out the above from my deseq2 analysis which has regions and those regions I have annotated it .</p>

<atac-seq>
Open chromatin region : nucleosome , mono-nucleosome , di-nucleosome
# Open chromatin region : nucleosome , mono-nucleosome , di-nucleosome

(Asked by: username_0 on 2020-12-29)

I m following this tutorial for my ATAC seq analysis. Here in of the section they show profile of open chromatin region plots.
**Plotting insert sizes with Greenleaf open, mono- and di-nucleosome profiles
This looks very similar to the image from the Greenleaf paper.

We can now annotate our nucleosome free (< 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) as in the Greenleaf study.**

Now Since i have done the differential accessibility analysis. So from that information how do i find out region which are "nucleosome free (< 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) ".
My objective

I would like to see this "nucleosome free (< 100bp), mono-nucleosome (180bp-247bp) and di-nucleosome (315-437) " from my differential accessible region which was done from readcount and then I ran deseq2.
My goal is not to plot the way they have done rather find those regions mentioned above and see it between different conditions.

Now based on the tutorial of that section where they do for each bam file.
Do i need to run the analysis on each bam file or each condition or Can I find out the above from my deseq2 analysis which has regions and those regions I have annotated it .

**Tags:** atac-seq

---

## 1 answers

### Answer by username_1 on 2020-12-30. Score: 1 | Accepted answer

It's impossible to find nucleosome free and such regions without a BAM file. Your only options are:

1. Filter the BAM files to contain only the size-range of interest and rerun the counting and statistics.
2. Find or write a counting program that can track fragment sizes.

BTW, those size ranges could use some adjustments, you're throwing away all fragments between 100 and 180 bases using their ranges.

> **Comments:**
> **username_0**: Thank you devon . "BTW, those size ranges could use some adjustments, you're throwing away all fragments between 100 and 180 bases using their ranges." how to do it add some region?  
> **username_1**: Looks at a size distribution plot and adjust accordingly.  

---

