genomecov VS bamCoverage

1

Entering edit mode

2.4 years ago

c.heininger ▴ 10

Hello,

can someone help me understand the differences between bedtools' genomecov and deepTools' bamCoverage? Are they used interchangeably, are there key differences, can I generate the same results with both tools? I have sequencing data from a paired-end PRO-seq experiment and for generating the coverage tracks I am not sure which tool to "trust", as I have not been able to reproduce their respective results with each other.

Any experience, thoughts, directions? I am grateful for any input!

Cheers

coverage NGS PRO-seq • 1.2k views

ADD COMMENT • link updated 2.4 years ago by Friederike 8.9k • written 2.4 years ago by c.heininger ▴ 10

1

Entering edit mode

If you need to visualize your data in a Genome Browser (such as IGV or UCSC) I would suggest you go with bamCoverage. I have never used genomecov but I never had problems in using bamCoverage.

ADD REPLY • link 2.4 years ago by Marco Pannone ▴ 790

0

Entering edit mode

I have not been able to reproduce their respective results with each other

Did you use a scaling factor with genomecov? What was the bin size you used for bamCoverage? Unless you share both lines of codes it's difficult to help you pin point the issue. Apart from the two factors "scaling" and "bin size", there are more details that may differ, such as the counting of fragments rather than reads for paired-end data.

ADD REPLY • link 2.4 years ago by Friederike 8.9k

Login before adding your answer.