[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