[Bioperl-l] poblems with object's locations
Hilmar Lapp
hlapp at gnf.org
Thu Jul 10 18:44:59 EDT 2003
On Thursday, July 10, 2003, at 05:43 AM, Marc Logghe wrote:
> Hi Sergiy,
> don't worry, it is only a typo:
> it is Bio::SeqIO, with a 'q' not with a 'g'
And for use lib only specify the root of the tree:
>>
>> use lib '/usr/lib/perl5/site_perl/5.8.0/Bio/SeqIO';
>>
>> use lib '/usr/lib/perl5/site_perl/5.8.0/Bio/Seq';
instead do
use lib '/usr/lib/perl5/site_perl/5.8.0';
which suspiciously looks like your standard site-perl path, so you
probably can leave that out completely as it's in perl's search path
already.
-hilmar
>>
>> $seqin = Bio::SegIO->new( -format =>'EMBL',-file=>'record.EMBL');
>>
>> $seqout = Bio::SeqIO->new( -format => 'Fasta', -file => '>output.fa');
>>
>> while ( my $seqobj->next_seq()){
>>
>> print "Seen sequence",$seqobj->display_id,"start of seq",
>>
>> substr($seqobj->seq,1,10),"\n";
>>
>> if ($seqobj->moltipe eq "dna") {
>>
>> $rev=$seqobj->revcom();
>>
>> $id = $seqobj->display_id();
>>
>> $id="$id.rev";
>>
>> $rev->display_id($id);
>>
>> $seqout->write_seq($rev);
>>
>> }
>>
>> foreach $feat ($seqobj->top_SeqFeature()){
>>
>> if ($feat->primary_tag eq 'exon')
>>
>> { print STDOUT "Location",$feat->start,":",
>>
>> $feat->end,"CFF[".feat->gff_string,"]\n";
>>
>> }}
>>
>> }
>>
>>
>>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
--
-------------------------------------------------------------
Hilmar Lapp email: lapp at gnf.org
GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
-------------------------------------------------------------
More information about the Bioperl-l
mailing list