Hi,
I have a bed file of annotated regions. I want to find genomic regions which are particularly dense with these annotations. I thought perhaps a way to do this is standard peak calling (Poisson model) and I was wondering if there are any tools out there that do this? Using MACS2 with the input bed file does not work because it keeps predicting that the read size (annotation feature size in my case) is 200. Any suggestions?
•
written
by
jrose •
1 answer
Seems like potential use-case for bedtools merge assuming these regions overlap. Merge and count regions.
Otherwise, a binning approach makes sense. Break genome up into 10 kb chunks, or whatever chunk size makes sense, intersect your regions with them, spit out counts for each bin, grab the top X of them or set an arbitrary threshold, etc.
Log in to answer this question.
-
MACS2 CallPeak Producing Empty Files When Specifying Genome size
written by Marcos •
-
How to annotate a MACS2 peak output file that is in tabular format?
written by Jen •
-
partition single bed file into non-overlapping regions
written by jrose •
-
Sorting .bed file by genes size
written by Marco Pannone •
-
MACS2 option for ATAC seq?
written by star
-
shift direction in ATAC/DNase data?
written by star
-
use Homer to annotation region of ChIP-seq peaks.
written by mikysyc2016
-
MACS2 MODEL VS NO MODEL FOR CHIRP SEQ
written by prathyushakonda •
-
From bed file to enrichment score
written by GK1610
-
Algorithmic Difference between MACS and MACS2
written by taraeicher •
