I would like to get the average coverage of all the captured bases in a bam file. What would be the best way to do this? What I am looking is a simple one number like 40X. Given that there may be millions of bases sequenced in a next gen study I would like to get the overall average coverage for all these bases as a single number.
i have used DepthOfcoverage module but i am not getting average coverage its giving coverage for each base
note: My bam file is only for one gene not for whole genome.
Please let me know