samtools coverage
1
0
Entering edit mode
4 months ago
Candela • 0

Hi guys, I have mapped metagenomic sequencing reads to a collection of assembled sequences that represent several MAGs (previously binned from the metagenome sample) obtaining my bam file, and my idea is to get the average coverage values within 5 Kbp sliding windows sampled every 100 bp in each MAG draft genome. I found the function samtools depth that gives you depth of coverage per base pair, but anyone can tell me how to obtain the average coverage of the desire region that i want.

Thanks in advance.

samtools bedtools coverage • 259 views
ADD COMMENT
0
Entering edit mode
4 months ago
yhoogstrate ▴ 110

I believe you can do that using the -b <bed> list of positions or regions argument and provide it a BED file with the regions of interest.

ADD COMMENT

Login before adding your answer.

Traffic: 2687 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