[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