This is a test version of Biostars. For the public version, visit https://www.biostars.org.
extract some feature from peak?

I did peak calling using MACS2. Now I like to calculate some features (like below) from them but I am not sure which tools can be helpful.

  • peak height
  • number of read / count for each peak -- > I think map / genomecoverage function of bedtools may be helpful but I am not sure!

Thanks for your help in advanced.

bedtools diffbind macs2 featurecount chip-seq

you may use some tools in bedtools suite- intersecting the bam with the intervals bed file by intersectBed, use coverage to compute depth and breadth, and even genomecov and an in house script.

0 answers

No answers yet.

Log in to answer this question.