[Biopython-dev] [Bug 2949] _parse_pdb_header_list: REVDAT is for oldest entry.

[bugzilla-daemon at portal.open-bio.org]

http://bugzilla.open-bio.org/show_bug.cgi?id=2949


crosvera at gmail.com changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
                 CC|                            |crosvera at gmail.com




------- Comment #1 from crosvera at gmail.com  2010-03-20 02:56 EST -------
(In reply to comment #0)
> [...]
>         elif key=="REVDAT":
>             #Modified by Paul T. Bathen to get most recent date instead of
> oldest date.
>             #Also added additional dict entries
>             if dict['release_date'] == "1909-01-08": #set in init
>                 rr=re.search("\d\d-\w\w\w-\d\d",tail)
>                 if rr!=None:
>                     dict['release_date']=_format_date(_nice_case(rr.group()))
> 
>                 dict['mod_number'] = hh[7:10].strip()
>                 dict['mod_id'] = hh[23:28].strip()
>                 dict['mod_type'] = hh[31:32].strip()

The Protein Data Bank Contents Guide (Version 3.20,
http://www.wwpdb.org/documentation/format32/sect2.html#REVDAT)  says that
modNum use the colums: 8-10. modId use the colums: 24-27. And modType use the
colum 32.
So the last part of your code should change to:

                 dict['mod_number'] = hh[7:9]
                 dict['mod_id'] = hh[23:26]
                 dict['mod_type'] = hh[31]


Regards.

--
Carlos Rios V.


-- 
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.