I've been running a blastp search to compare the Histone H1t amino acid sequences of Ursus Maritimus and Mus Musculus, and have run into a bit of a problem. When calculating the E() Value, the BLAST program gives a value of 1e-58, which would be fine, but when I try to calculate the E() Value myself using the equation:
Where E is the expected value, m is the total number of residues in the database/subject species (in this case 212), n is the number of residues in the query sequence (209) and S’ is the bit score (171).
Which lead to a problem for me. The bit score is calculated correctly from both BLAST and my calculations, so I'm wondering what the problem could be?
Any help at all would be greatly appreciated! Minas