I am thinking that the most important quantification parameter must be how many reads are mapped within annotated regions?
Hence it is not mapped reads per whole genome, but rather mapped reads per annotated region, like so:
coverage = n_mapped_reads * read_length / (n_annotated_regions * average_length_of_annotated_regions)
This way we get average number of times each annotated nucleotide was sequenced?
L
Hence it is not mapped reads per whole genome, but rather mapped reads per annotated region, like so:
coverage = n_mapped_reads * read_length / (n_annotated_regions * average_length_of_annotated_regions)
This way we get average number of times each annotated nucleotide was sequenced?
L
Comment