[Bioperl-l] parsering signal-p output
Diego Chaves Moreno
diegochavesm at gmail.com
Thu Apr 8 20:05:02 UTC 2010
Hi everyone,
im a new user of bioperl modules. I tried to make a parser to signal-p output using Bio::Tols::Signalp (code forward). The parser is working and i need use the information in $seq_id variable (in red) but when i tried to print this variable the value is empty, the same is with $name_seq variable (blue).
Thanks,
Diego Chaves Moreno
Los Andes University
Bogota-Colombia
code
========================
#!/usr/bin/perl
use strict;
use Bio::Tools::Signalp;
my $file = $ARGV[0];
my $parser = Bio::Tools::Signalp->new(-file =>$file );
while( my $sp_feat = $parser->next_result()) {
my $start = $sp_feat->start();
my $end = $sp_feat->end();
my $len = $sp_feat->length();
my $strand = $sp_feat->strand();
my $score = $sp_feat->score();
my $frame = $sp_feat->frame();
my $featname = $sp_feat->display_name();
my $seq_id = $sp_feat->seq_id();
my $whole_seq = $sp_feat->entire_seq();
my $name_seq = $sp_feat->seq_id();
while(my $location = $sp_feat->location()){
my $location_type = $location->location_type();
my $start_L = $location->start();
my $end_L = $location->end();
my $minstart = $location->min_start();
my $maxstart = $location->max_start();
my $start_pos_type = $location->start_pos_type();
my $minend = $location->min_end();
my $maxend = $location->max_end();
my $seqidL = $location->seq_id();
More information about the Bioperl-l
mailing list