I'm trying to find out the average depth of coverage per target region in a bed file. I can easily use bedtools -depth to find the depth of coverage at each base, but I'd much rather have it as average coverage per region. The bedtools -hist option produces 110 reports on one of my features. I only want 1 number (average coverage) per target region. Any way to do this?
I found threads from a few years ago but wonder if there is anything current. One of the solutions doesn't compile on my Ubuntu server and the bash script gets a "fatal: division by zero attempted"
1 answer
Please take a look at my answer here: Compute mean depth coverage for exome data with paired end, overlapping, features I believe that you require BEDTools version > 2.0.
Also some simple code for plotting the profile: Histogram showing overlap between 2 bed files
Kevin
Log in to answer this question.
-
awk: cmd. line:1: fatal: division by zero attempted when calculating average depth when using subpr…
written by Manuel •
-
How to obtain zero-based coordinates read depth using bedtools coverage for a specific region?
written by vladimiralan •
-
Difference in .bed target region base count and mapped base count for different tools
written by gavehan •
-
GC Bias in Sequence Capture - Analysis & Quantification
written by jnowacki
-
How to calculate average depth for each Amplicon in targeted amplicon sequencing?
written by bioinfo89 •
-
calculate average target coverage from bam file
written by bioguy24
-
Which Read-Depth is Needed to Find Insertions?
written by michaela_boell
-
samtools bedcov vs. bedtools coverage
written by Dataman
-
calculation of coverage for a given region
written by tonja.r
-
Coverage plots of gvcf files
written by Floris Brenk
