[Biopython-dev] [Bug 2481] New: bitscore not parsed.
bugzilla-daemon at portal.open-bio.org
bugzilla-daemon at portal.open-bio.org
Thu Apr 3 20:00:56 EDT 2008
http://bugzilla.open-bio.org/show_bug.cgi?id=2481
Summary: bitscore not parsed.
Product: Biopython
Version: 1.45
Platform: PC
OS/Version: Linux
Status: NEW
Severity: normal
Priority: P2
Component: Main Distribution
AssignedTo: biopython-dev at biopython.org
ReportedBy: sbassi at gmail.com
>>> from Bio.Blast import NCBIXML
>>> fr=NCBIXML.parse(open('/media/disk/GENES/INTA9/BLAST/seqspMOS.xml')).next()
>>> fr.descriptions[0].title
u'gnl|BL_ORD_ID|0 pMOS (vector para mt INTA)'
>>> fr.descriptions[0].bits
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
AttributeError: Description instance has no attribute 'bits'
To fix this, two files must be modified:
NCBIXML.py
and
Record.py
In NCBIXML.py, 2 changes:
Line 94 : "method = self._secure_name('_end_' + name.replace("-","_"))" #name
in the xml file is "bit-score" and the function should be named like this but
can only be named "_end_Hsp_bit_score" hence change from - to _ resolve the
issue and should not disturb the rest, this method could/should be also applied
line 63 for StartElement
Lines 409-410 uncommented
In Record.py:
Add "self.bits = None" # in line 68
This bug was reported to me by Yoan Jacquemin when testing my code to convert
BLAST XML output to HTML.
After applying this modifications, it works:
>>> from Bio.Blast import NCBIXML
>>> f_in='/mnt/hda2/bio/3vsT.xml'
>>> fr=NCBIXML.parse(open(f_in)).next()
>>> fr.descriptions[0].bits
32.210500000000003
--
Configure bugmail: http://bugzilla.open-bio.org/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are the assignee for the bug, or are watching the assignee.
More information about the Biopython-dev
mailing list