[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