I have several different DNA exome and RNA sequences from the Illumina GA2 platform. For verification purposes and quality control I'd like to know how many bases the aligned sequences cover (I've already got aligned SAM files produced by Bowtie for DNA and TopHat for RNA). From the base counts I'd then like to be able to figure out the sequences' % coverage of the entire genome. Does anyone have any advice on how to go about doing this in the most efficient way possible?
Thanks!
Thanks!