Coverage for exon

[nvteja]

Hi everyone,

I want to calculate coverage for exon sequences using 454AlignmentInfo.tsv file. I have a bed file of exons. So to calculate coverage is it enough if i average the unqiue depth of every base in the exon to get total coverage for that exon ? or should I use average of total depth for each base ?


Thanks,
Teja.

[flxlex]

Unique depth is more correct, as Total depth includes the duplicate reads. These are (usually) the results of two beads ending up in the same (micro)reaction during emPCR, resulting in them obtaining the exact same ssDNA strand, resulting in basically the same flowpattern.