在biopython中,我怎样才能获得100%的身份对齐?

时间:2014-11-05 21:46:18

标签: bioinformatics biopython

我找到了获得hsp身份的方法,但是没有对齐。一个对齐可以有多个hsp,这意味着hsp的得分(身份)不等于对齐。 当我在网页上爆炸时,我总是得到对齐的分数和身份。 这是我提出的任何逻辑问题,谁可以解释这个?

1 个答案:

答案 0 :(得分:0)

网页为您提供了某些HSP的身份,而不是完整的对齐方式。例如。这两个HSP:

1203 bits(651)  0.0     730/769(95%)    1/769(0%)   Plus/Minus
  106 bits(57)  6e-19     70/76(92%)    2/76(2%)    Plus/Minus

计算:

Max     Total   Cover   eVal ident
1203    1309    100%    0.0  95%

因此网络只是添加得分(1203 + 106 = 1309)并且可能平均eval和ident(730 + 30)/(769 + 76)= 0.9467~0.95。所以这就是你模仿网络BLAST“加入”统计数据的方式。根据记录,这两个HSP之间存在差距,web-BLAST未包含在统计数据中。