How to create a bedGraph file for broadPeak peak caller?
1
0
Entering edit mode
6.4 years ago
mrs.hope ▴ 10

I am trying to use broadPeak package (http://jordan.biology.gatech.edu/page/software/broadpeak/) on my ChIP-seq data. As an input file I should provide "...The input file of the sorted ChIP-seq read-mapping profile in the genome needs to be in bedGraph format. The four tab-delimited columns are chromosome, start, stop and tag count. ...... The size of each record (the bin) should be equal (e.g. 200bp, if you scan the genome by dividing it to 200bp non-overlapping bins)..".

I wonder if there is an established and straightforward way to create such bedGraph file with windows of a fixed size from a sorted bam file? I failed to find some ready-to-use scripts a la bamtoBedGraph conversion. As far as I understand the one from the bedtools genomeCoverage doesn't allow you to choose windows where you want to count reads.

Does anyone know the way how to do it relatively easy and fast? My only guess is to create a bed file with a genome divided into regions of required size. Convert it into GFF and then annotate with the reads count. But I think that it will be very slow and I guess not very elegant way to solve the problem.

Any piece of advice will be appreciated!

 

broadPeak bedGraph ChIP-Seq peak calling • 2.9k views
ADD COMMENT
0
Entering edit mode
5.6 years ago
569497735 • 0

It seems that bedtools can do this. Check this website, https://gencore.bio.nyu.edu/variant-calling-pipeline/

ADD COMMENT

Login before adding your answer.

Traffic: 1918 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6