[Bioperl-l] possible to skip parsing features when calling Bio::SeqIO::new?

Jason Stajich jason.stajich at duke.edu
Mon Jan 17 21:32:06 EST 2005


See the docs for Bio::Seq::SeqBuilder

I think this will work:
my $seqIterator = Bio::SeqIO->new("-file"   => "$file",
                                     "-format" => "genbank");
$seqIterator->sequence_builder->add_unwanted_slot('features');

If you additionally don't want the Annotations (references,etc)
$seqIterator->sequence_builder->add_unwanted_slot('features', 
'annotation');

[don't ask why one is plural and other singular... =)]

-jason
On Jan 17, 2005, at 9:00 PM, Scott Markel wrote:

> I'm using BioPerl 1.4 to read a GenBank sequence file.  Is there
> an option or parameter I can set in the
>
> my $::seqIterator = Bio::SeqIO->new("-file"   => "$file",
>                                     "-format" => "genbank");
>
> call that will cause the parser to skip the features?
>
> I checked BioPerl 1.5RC2 and didn't see any changes there
> that would address my question.
>
> Scott
>
> -- 
> Scott Markel, Ph.D.
> Principal Bioinformatics Architect  email:  smarkel at scitegic.com
> SciTegic Inc.                       mobile: +1 858 205 3653
> 9665 Chesapeake Drive, Suite 401    voice:  +1 858 279 8800, ext. 253
> San Diego, CA 92123                 fax:    +1 858 279 8804
> USA                                 web:    http://www.scitegic.com
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/



More information about the Bioperl-l mailing list