[Biopython] hsp.identities (solved)

Mike Williams dmikewilliams at gmail.com
Mon Aug 10 18:23:31 UTC 2009


Well, one hour after posting my question, I found the answer in the
list archives:
http://portal.open-bio.org/pipermail/biopython-dev/2006-April/002347.html
<quote from list message>
What happens is that if the Blast output looks
like this:

 Identities = 28/87 (32%), Positives = 44/87 (50%), Gaps = 12/87 (13%)

then the text-based parser returns
hsp.identities = (28, 87)
hsp.positives = (44, 87)
hsp.gaps = (12, 87)

while the XML parser returns
hsp.identities = 28
hsp.positives = 44
hsp.gaps = 12
; we can get the 87 from len(hsp.query).
</quote>

Cheers,

Mike



More information about the Biopython mailing list