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

Scott Markel smarkel at scitegic.com
Mon Jan 17 21:53:05 EST 2005


Jason,

Excellent!  Thank you.

Scott

Jason Stajich wrote:

> 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
>>
> -- 
> Jason Stajich
> jason.stajich at duke.edu
> http://www.duke.edu/~jes12/
> 
> 

-- 
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



More information about the Bioperl-l mailing list