[Biopython] Fwd: Re: BioJava-like seqres alignment for Bio.PDB
Bryan Lunt
lunt at ctbp.ucsd.edu
Wed Jun 30 18:13:37 UTC 2010
It is indeed fiddly,
But I needed it immediately.
http://gist.github.com/459005
Basically you call
mapAll( open("pdb1234.ent"), structureObjectPreviouslyReadWith-Bio.PDB)
and get back a dictionary of chain IDs -> lists of tuples, where
position in the list corresponds to position in the seqres.
It _is_ ugly, and makes my eyes bleed.
Sorry.
-Bryan
More information about the Biopython
mailing list