Correct percent identity based on gene coverage.
1
0
Entering edit mode
2.8 years ago
genya35 ▴ 40

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

alignment • 674 views
ADD COMMENT
0
Entering edit mode
2.8 years ago
Joe 19k

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

Normalised ID = ID x (Coverage / 100)

ADD COMMENT

Login before adding your answer.

Traffic: 1594 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6