View Single Post
Old 10-29-2013, 01:48 PM   #3
David Eccles (gringer)
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838

I would use bedtools (which can use BAM files as input). Use genomecov to get the coverage all across the genome:

Or just coverage if you have a specific gff/bed file with the region that you want to determine the coverage for:

edit: sorry, that's just for the base-pair coverage (independent of the base). I suspect that per-nucleotide counts in a particular region is straying into "needs a custom script to work it out from the samtools mpileup output" territory.

Last edited by gringer; 10-29-2013 at 01:53 PM.
gringer is offline   Reply With Quote