Depth of coverage of Bed file regions?
0
0
Entering edit mode
5.1 years ago
Jon17 ▴ 20

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 • 1.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 1960 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6