jsh58 / Genrich

Detecting sites of genomic enrichment
MIT License
185 stars 27 forks source link

HowTo call peaks from bigwig files? #103

Closed malcook closed 7 months ago

malcook commented 1 year ago

I am involved in a collaboration in which .bam files are apparently no longer available and only bigwig files can be found. Do you have any advice/strategy as to whether Genrich can be used to call peaks from them? I fully understand there may be multiple reasons why one should not want to do this, but still, here we are. Thanks!

jsh58 commented 1 year ago

Thanks for the question. No, Genrich cannot call peaks from bigwig files. Even in -P mode, Genrich requires that the input log file have p-/q-values, which a bigwig file won't have.

keenhl commented 10 months ago

You can convert bigwig to bedgraph, then do peak calling with MACS.