Question: Merging peak files for unionized set
0
gravatar for aam458
3 months ago by
aam4580
aam4580 wrote:

Hello,

I'm relatively new to atac and pretty much all of bioinformatics work but I've been picking up a lot of experience as I've gone and figured most things out on my own.

I'm beginning some work on some bed files I've (finally) attained for atac-seq analysis, I've used bedtools intersect to check how well the replicates overlap for downstream analysis (decently well at ~74%) and now I want to merge them (each replicate into a union set) so I can compare different timepoints.

Does anyone have any thoughts/experiences/procedures as to how to correctly merge two replicate peak bed files and maybe some hints as to common pitfalls to avoid? It occurs to me that it would be very hard to tell on first glance whether one has made a large mistake.

I was thinking about using bedtools merge function to do it without specifying any options (because none of them seem to apply, looking at the bedtools documentation).

Thanks!

atac-seq next-gen bedtools • 212 views
ADD COMMENTlink modified 3 months ago by trausch1.2k • written 3 months ago by aam4580
1
gravatar for trausch
3 months ago by
trausch1.2k
Germany
trausch1.2k wrote:

I am not sure if that's necessarily the best solution but we indeed use bedtools merge to get a unified peak list, then count fragments for each sample at these peaks and then run differential peak calling.Here is the link to our ATAC-Seq workflow:

https://github.com/tobiasrausch/ATACseq

ADD COMMENTlink written 3 months ago by trausch1.2k
Please log in to add an answer.

Help
Access

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