[Biopython] Sequence annotation (Features)

[Peter]

On Sat, Jan 16, 2010 at 8:38 AM, Anirban Bhattachariya <anbhat at utu.fi> wrote:
> Hi,
>
> I'm trying to download a protein sequence object (using ID or
> accession number) and then trying to print its variants (all
> variant sequences) from its features and annotations.I'm using
> pseudocholinesterase (http://www.uniprot.org/uniprot/P06276
> as an example since it has lot of natural variants.
>
> The problem is when I'm trying to access the features its
> saying "0 features" ; how can I access the features in
> Swiss-Prot file  like in genbank file format ( as in section
> 4.6  of the tutorial).

It's a know missing feature, although there is a patch here:
http://bugzilla.open-bio.org/show_bug.cgi?id=2235

You could help with testing/improving the patch in order to
get Bio.SeqIO to do this in future, or in the short term use
the underlying parser in Bio.SwissProt.

Regards,

Peter