Question: Using Diffbind on ATAC data - adapter offset issue
gravatar for YaGalbi
13 months ago by
MRC Harwell Institute, Oxford, UK
YaGalbi20 wrote:



I'd like to use Diffbind on Atac data. However what is done about the adapter adjustment that occurs before peak calling? For peak calling all forward reads are offset by +4 bases and reverse reads by -5 so aligning peak positons to BAM file positions of the reads may need reversing of the offset.

Forward strand example

Position:     1   2   3   4   5   6   7   8   9   10 

Base:          A   G   A   T  A   C  C   A   A    T

BAM read:        G--------------C                                                  Aligner places G---C read at positions 2 to 6

Peak:                                   6----------------10                           Peak calling after offset calls a peak at 6-10


Does anyone have experience with this? Am I thinking correctly? It should be a simple fix, I'm just wondering if I am thinking correctly and this is needed.




ADD COMMENTlink modified 13 months ago by Rory Stark2.4k • written 13 months ago by YaGalbi20
gravatar for Rory Stark
13 months ago by
Rory Stark2.4k
CRUK, Cambridge, UK
Rory Stark2.4k wrote:

Let's consolidate discussion of this issue here:

DiffBind for ATAC-seq


ADD COMMENTlink written 13 months ago by Rory Stark2.4k
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.2.0
Traffic: 362 users visited in the last hour