Entering edit mode
5.8 years ago
dec986
▴
370
I have some aligned BAM files and I want to calculate the sequencing depth, but numerous google searches have not produced any results. I could generate a tool by myself, but I want to use something which has been used by others so that there is less room for error.
What was wrong with, say, samtools depth? Or bedtools coverage? if you can explain why those were not suitable people might be better able to direct you to what is suitable.
In addition, are you interested in read-depth over a particular genomic position?; are you interested in the target depth of coverage to which your samples were sequenced?; the number of bases that reached a particular depth of coverage (e.g. 10x, 20x, 30x, etc?).
You must be more specific.