[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