Question

Per Base Coverage

0

Entering edit mode

13.1 years ago

Info_Nowise ▴ 10

Is there a way to obtain per-base coverage for a define chromosome interval using a bam file generated from Illumina single-end reads? genomeCoverageBed in Bedtools does not seem to have an option for it.

bedtools coverage • 4.4k views

ADD COMMENT • link updated 13.1 years ago by Aaronquinlan 12k • written 13.1 years ago by Info_Nowise ▴ 10

Ram · Answer 1 · 2012-03-09

4

Entering edit mode

13.1 years ago

Zev.Kronenberg 12k

samtools depth your.bam

ADD COMMENT • link updated 5.1 years ago by Ram 45k • written 13.1 years ago by Zev.Kronenberg 12k

Ram · Answer 2 · 2012-03-09

4

Entering edit mode

13.1 years ago

Aaronquinlan 12k

Also,

bedtools coverage -abam aln.bam -b regions.bed -d

ADD COMMENT • link updated 5.1 years ago by Ram 45k • written 13.1 years ago by Aaronquinlan 12k