Skip to content
Snippets Groups Projects
Commit 0c413f87 authored by juguang's avatar juguang
Browse files

enable to convert Bio::SeqFeature::FeaturePair to Bio::EnsEMBL::ProteinFeature

parent e3bb16f0
No related branches found
No related tags found
No related merge requests found
......@@ -48,6 +48,7 @@ use vars qw(@ISA);
use Bio::EnsEMBL::FeaturePair;
use Bio::EnsEMBL::RepeatConsensus;
use Bio::EnsEMBL::ProteinFeature;
use Bio::EnsEMBL::Utils::Converter;
use Bio::EnsEMBL::Utils::Converter::bio_ens;
@ISA = qw(Bio::EnsEMBL::Utils::Converter::bio_ens);
......@@ -73,6 +74,8 @@ sub _convert_single {
return $self->_convert_single_to_repeatFeature($pair);
}elsif($self->out eq 'Bio::EnsEMBL::FeaturePair'){
return $self->_convert_single_to_featurePair($pair);
}elsif($self->out eq 'Bio::EnsEMBL::ProteinFeature'){
return $self->_convert_single_to_proteinFeature($pair);
}else{
my $output_module = $self->out;
$self->throw("Cannot covert to [$output_module]");
......@@ -94,6 +97,16 @@ sub _convert_single_to_featurePair {
return $ens_fp;
}
sub _convert_single_to_proteinFeature {
my ($self, $pair) = @_;
my $featurePair = $self->_convert_single_to_featurePair($pair);
my $proteinFeature = Bio::EnsEMBL::ProteinFeature->new(
-feature1 => $featurePair->feature1,
-feature2 => $featurePair->feature2
);
return $proteinFeature;
}
sub _convert_single_to_repeatFeature {
my ($self, $pair) = @_;
my $feature1 = $pair->feature1;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment