genomecov VS bamCoverage
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
• 2,214 views
•
link
updated
by
Friederike
•
written
by
c.heininger •
•
written
by
c.heininger •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Selection Pressure Analysis
written by upasana •
-
File formats for ATACseq data
written by c.eskiw •
-
Remove reads from FASTQ file based on missing fixed base
written by c.heininger •
-
How to obtain zero-based coordinates read depth using bedtools coverage for a specific region?
written by vladimiralan •
-
noisy chip-seq peaks
written by kolikolmakani14 •
-
ATAC-seq coverage tracks using bamCoverage from Deeptools
written by mickey_95 •
-
Predictive models from RNA (co)expression data
written by mattia.bosio85 •
-
how to process spike-in with bamCoverage and plotHeatmap (with deeptools)
written by alessandro.pastore •
-
quantification of PolII Ser-5 enrichment per gene
written by Sentinel156
-
define regions of coverage and count snps
written by rob234king
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.
Did you use a scaling factor with
genomecov? What was the bin size you used forbamCoverage? 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.