Question: How to create a bedGraph file for broadPeak peak caller?
0
gravatar for mrs.hope
3.4 years ago by
mrs.hope10
Russian Federation
mrs.hope10 wrote:

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!

 

ADD COMMENTlink modified 2.5 years ago by 5694977350 • written 3.4 years ago by mrs.hope10
0
gravatar for 569497735
2.5 years ago by
5694977350
5694977350 wrote:

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

ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by 5694977350
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: 2138 users visited in the last hour