[BioPython] GenBank.FeatureParser and multiple accessionnumbers

[Jason Stajich]

On Thu, 24 Oct 2002, Jeffrey Chang wrote:

> I agree.  Perhaps the first accession number should be seq.id, and
> *all* the accessions (including the first one) should be an annotation
> like seq.annotations["accessions"].  They certainly should be
> accessible.
>
> How does bioperl handle this?
>
separate methods for declared 'RichSequences' (GenBank/EMBL/SwissProt
parser produce these which have a richer interface than basic sequences
with features).

In our case we have
get_secondary_accessions() which returns a list of accessions
to distinguish primary from secondary accessions.

> Jeff
>
>
> On Thu, Oct 24, 2002 at 04:18:30PM +0200, Andreas Kuntzagk wrote:
> > Hi,
> >
> > at the moment the GenBank.FeatureParser stores only the first accession
> > number of an entry. (right?)
> > Would it make sense to store all of them? (e.g. as an annotation)
> > I would like to search a dictionary of SeqRecords by an accession
> > number.
> >
> > ciao, Andreas "murple" Kuntzagk
> >
> > --
> > Andreas Kuntzagk <andreas.kuntzagk@mdc-berlin.de>
> >
> > _______________________________________________
> > BioPython mailing list  -  BioPython@biopython.org
> > http://biopython.org/mailman/listinfo/biopython
> _______________________________________________
> BioPython mailing list  -  BioPython@biopython.org
> http://biopython.org/mailman/listinfo/biopython
>

-- 
Jason Stajich
Duke University
jason at cgt.mc.duke.edu