This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Depth of coverage of Bed file regions?

Best way to find average depth of coverage for very specific regions via a bed file?

I've tried this: bedtools genomecov -i $BED -ibam $BAM

But my bam file says this:

chr1 69090 69190
chr1 69190 69290
chr1 69290 69390
chr1 69390 69490

and the output is this:

chr1 0 235026752 248956422 0.944048
chr1 1 4040485 248956422 0.0162297
chr1 2 983706 248956422 0.00395132
chr1 3 565978 248956422 0.0022734

which isn't helpful.

I'm hoping to get something similar to this:

chromosome Start stop Depth_of_coverage
chr1 69090 69190 45
chr1 69190 69290 46
chr1 69290 69390 45
chr1 69390 69490 30

bedtools bed genome coverage

0 answers

No answers yet.

Log in to answer this question.