Minimum coverage info on Miseq

[ramujana]

Hi,
does anyone can suggest how to extract MINIMUM coverage information of particular ROI after MiSeq run? (after standard fastq aligment to reference genome+manifest file and generating bam files)

[westerman]

It depends on how longer your ROI (region of interest) is. If not too long you should be able to do multiple '-c' (count) via samtools one base at a time and then choose the smallest count. I haven't tried this and it would take a bit of shell scripting but it should be feasible as a one-off solution.

[WhatsOEver]

"samtools depth" gives you the per base coverage and has a "-r" parameter with which you can define your ROI.
Alternatively, genomecov from bedtools works quite nicely for such questions (http://bedtools.readthedocs.org/en/l...genomecov.html)

[westerman]

'WhatsOEver' has the much better solution. I knew that samtools could do the work -- it is a rather powerful program -- but how to do it was escaping me in my work-induced haste of reading SeqAnswers.