Correct percent identity based on gene coverage.

Question

[genya35]

Hello,

When I run blastn on a sequence of a specific gene it returns percent identity. However, I would like to be able to correct for the coverage of the gene (How well the read covers the identified gene targeted by the primers). For example, if the coverage of the gene is 98% I would like to correct the % identity for the missing part of sequence. Could someone please suggest what formula to use?

Thanks

Answer 1

[Joe]

It should be sufficient to simply scale ID by coverage as a percent, e.g:

Normalised ID = ID x (Coverage / 100)