Bio::Tools
Coil
Toolbar
Summary
Bio::Tools::Coil
Package variables
No package variables defined.
Included modules
Inherit
Synopsis
use Bio::Tools::Coil
my $parser = new Bio::Tools::Coil();
while( my $sp_feat = $parser->next_result($file) ) {
#do something
#eg
push @sp_feat, $sp_feat;
}
Description
Methods
Methods description
Title : _add_prediction() Usage : $obj->_add_prediction($feat) Function: internal Example : Returns : |
Title : _predictions_parsed Usage : $obj->_predictions_parsed Function: internal Example : Returns : TRUE or FALSE |
Title : _read_fasta Usage : obj->_read_fasta($file) Function: Internal(not to be used directly) Returns : Args : |
Title : _result Usage : $feat = $obj->_result() Function: internal Example : Returns : |
Title : create_feature Usage : obj->create_feature(\%feature) Function: Internal(not to be used directly) Returns : Args : |
Title : next_result Usage : while($feat = $coil->next_result($file)) { # do something } Function: Returns the next protein feature of the coil output file Returns : Args : |
Title : parse_results Usage : obj->parse_results Function: Parses the coil output. Automatically called by next_result() if not yet done. Example : Returns : |
Methods code
sub _add_prediction
{ my ($self, $gene) = @_;
if(! exists($self->{'_feats'})) {
$self->{'_feats'} = [];
}
push(@{$self->{'_feats'}}, $gene); } |
sub _predictions_parsed
{ my ($self, $val) = @_;
$self->{'_preds_parsed'} = $val if $val;
if(! exists($self->{'_preds_parsed'})) {
$self->{'_preds_parsed'} = 0;
}
return $self->{'_preds_parsed'}; } |
sub _read_fasta
{ local (*FILE) = @_;
my( $id , $seq , %name2seq); while (<FILE>) {
chomp; if (/^>(\S+)/) {
my $new_id = $1;
if ($id) {
$name2seq{$id} = $seq;
}
$id = $new_id ; $seq = "" ;
}
elsif (eof) {
if ($id) {
$seq .= $_ ; $name2seq{$id} = $seq;
}
}
else {
$seq .= $_
}
}
return %name2seq;
}
1; } |
sub _result
{ my ($self) = @_;
return undef unless(exists($self->{'_feats'}) && @{$self->{'_feats'}});
return shift(@{$self->{'_feats'}}); } |
sub create_feature
{ my ($self, $feat) = @_;
my $feature = Bio::SeqFeature::Generic->new(-seq_id => $feat->{name},
-start => $feat->{start},
-end => $feat->{end},
-score => $feat->{score},
-source => $feat->{source},
-primary => $feat->{primary},
-logic_name => $feat->{logic_name},
);
$feature->add_tag_value('evalue',0);
$feature->add_tag_value('percent_id','NULL');
$feature->add_tag_value("hid",$feat->{primary});
return $feature; } |
sub new
{ my($class,@args) = @_;
my $self = $class->SUPER::new(@args);
$self->_initialize_io(@args);
return $self; } |
sub next_result
{
my ($self,$resfile) = @_;
my $gene;
$self->parse_results($resfile) unless $self->_predictions_parsed();
$gene = $self->_result();
return $gene; } |
sub parse_results
{ my ($self,$resfile) = @_;
my $filehandle = $resfile;
my %result_hash =_read_fasta($filehandle); my @ids = keys %result_hash;
my @feats;
foreach my $id (keys %result_hash){
my $pep = reverse ($result_hash{$id});
my $count = my $switch = 0;
my ($start, $end);
while (my $aa = chop $pep) {
$count++;
if (!$switch && $aa eq "x") {
$start = $count;
$switch = 1;
}
elsif ($switch && $aa ne "x") {
$end = $count-1;
my (%feature);
$feature{name} = $id;
$feature{start} = $start;
$feature{end} = $end;
$feature{source} = "Coils";
$feature{primary} = 'ncoils';
($feature{program}) = 'ncoils';
$feature{logic_name} = 'Coils';
my $new_feat = $self->create_feature (\%feature);
$self->_add_prediction($new_feat);
$switch = 0;
}
}
}
$self->_predictions_parsed(1); } |
General documentation
User feedback is an integral part of the evolution of this and other
Bioperl modules. Send your comments and suggestions preferably to
the Bioperl mailing list. Your participation is much appreciated.
bioperl-l@bioperl.org - General discussion
http://bioperl.org/MailList.shtml - About the mailing lists
Report bugs to the Bioperl bug tracking system to help us keep track
of the bugs and their resolution. Bug reports can be submitted via
email or the web:
bioperl-bugs@bio.perl.org
http://bugzilla.bioperl.org/
The rest of the documentation details each of the object methods.
Internal methods are usually preceded with a _