Hi, I have a bed file containing the coordinates of certain regions of the genome. I want to get the genome coverage of these regions at each base position (i,e, similar to using -d option in bedtools genomecov).
example bed file
chr1 100 200 chr5 250 300
I want the output in the following manner:
chr1 100 "coverage at this position" chr1 101 "coverage" chr1 102 "coverage" . . . chr1 199 "coverage" chr1 200 "coverage" chr5 250 "coverage" chr5 251 "coverage" ...and so on
The genomecov of bedtools does this for the entire genome. Is there any way to get the coverage within the regions specified only?
•
written
by
ilovesuperheroes1993 •
1 answer
Mosdepth for the win (https://github.com/brentp/mosdepth ):
-b --by <bed|window> optional BED file or (integer) window-sizes.
•
written
by
ATpoint
Log in to answer this question.
-
Identifying RBP motifs overlapping exons
written by RK •
-
Efficiently find regions in bed file that are within start and end regions of another bed file.
written by rxs1018 •
-
partition single bed file into non-overlapping regions
written by jrose •
-
Merge regions in bedtools genomecov/bedgraph file
written by eoneill627 •
-
Split bed file by unique first three columns
written by graeme.thorn
-
Consensus of multiple BED files
written by bruce.moran
-
Merge bed files, keep the peak count
written by rioualen
-
compute count of position coverage with bedtools
written by sacha
-
about INCLUSIVE intersections of GRanges
written by Bogdan
-
How To Melt A Bed File (With Bedops?)
written by Irsan
https://bedtools.readthedocs.io/en/latest/content/tools/coverage.html