[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