[Bioperl-l] SignalP bug in predict_protein_features?
david.vilanova at urbanet.ch
david.vilanova at urbanet.ch
Wed Jul 2 16:14:01 EDT 2003
Hello,
I was checking the Bio::Tools::Run::Signalp module and saw that in the code (in
the predict_protein_features) there are 2 different sizes for the subseq
depending on weather a sequence is passed or a file.
I guess the size should be the same right (meaning 40 or 50) ?
if ($seq->length>50){
my $sub_seq = $seq->subseq(1, 50);
$seq->seq($sub_seq);
}
while ( my $tmpseq = $in->next_seq() ) {
my $sub_seq = $tmpseq->subseq(1,40);
$tmpseq->seq($sub_seq);
$infile1 = $self->_writeSeqFile($tmpseq);
}
More information about the Bioperl-l
mailing list