What is the proper way to calculate coverage from the total amount of reads per gene/exon target? I used bedTools coverageBed to output the per base reads in for a target in a bed file. I then used awk to average all reads per target and filter those to specific genes. Now I need to calcuulate the x coverage using the reads. I am not sure what the best formula is?
Basically, in the file I have there are a list of genes with reads, say 300. I am trying to use the reads to calculate coverage. Thank you
.
Basically, in the file I have there are a list of genes with reads, say 300. I am trying to use the reads to calculate coverage. Thank you
.
Comment