[Bioperl-l] Removing SeqFeatures
Marc.Logghe at devgen.com
Thu Jun 2 09:23:34 EDT 2005
I think the only way is to run the remove_SeqFeatures() and catch the
removed features in an array.
Filter the array keeping the features you want to keep and add them to
your sequence object using
> -----Original Message-----
> From: bioperl-l-bounces at portal.open-bio.org
> [mailto:bioperl-l-bounces at portal.open-bio.org] On Behalf Of
> michael watson (IAH-C)
> Sent: Thursday, June 02, 2005 2:54 PM
> To: Bioperl-l at portal.open-bio.org
> Subject: [Bioperl-l] Removing SeqFeatures
> I want to take an EMBL sequence and remove some of the
> features, but keep all of the other information (accession,
> id, sequence, comments etc).
> So I have got some code which gives me the features that I
> want and the ones I don't want, but I can't figure out how to
> "remove" a feature.
> The only way I can think of doing this is to create a new
> sequence and transfer the data across from the old sequence
> bit by bit - surely there's a better way?
> Many thanks
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
More information about the Bioperl-l