Bio::EnsEMBL::Variation
AlleleFeature
Toolbar
Summary
Bio::EnsEMBL::Variation::AlleleFeature - A genomic position for an allele in a sample.
Package variables
No package variables defined.
Included modules
Inherit
Synopsis
# Allele feature representing a single nucleotide polymorphism
$af = Bio::EnsEMBL::Variation::AlleleFeature->new
(-start => 100,
-end => 100,
-strand => 1,
-slice => $slice,
-allele_string => 'A',
-variation_name => 'rs635421',
-variation => $v);
...
# a allele feature is like any other ensembl feature, can be
# transformed etc.
$af = $af->transform('supercontig');
print $af->start(), "-", $af->end(), '(', $af->strand(), ')', "\n";
print $af->name(), ":", $af->allele_string();
# Get the Variation object which this feature represents the genomic
# position of. If not already retrieved from the DB, this will be
# transparently lazy-loaded
my $v = $af->variation();
Description
This is a class representing the genomic position of a allele in a sample
from the ensembl-variation database. The actual variation information is
represented by an associated Bio::EnsEMBL::Variation::Variation object. Some
of the information has been denormalized and is available on the feature for
speed purposes. A AlleleFeature behaves as any other Ensembl feature.
See
Bio::EnsEMBL::Feature and
Bio::EnsEMBL::Variation::Variation.
Methods
Methods description
Arg [1] : string $newval (optional)
The new value to set the allele attribute to
Example : $allele = $obj->allele_string()
Description: Getter/Setter for the allele attribute.
Returntype : string
Exceptions : none
Caller : general
Status : At Risk |
Arg [1] : reference to string $seqref
Arg [2] : int $start of the seq_ref
Example : $sequence = 'ACTGAATATTTAAGGCA';
$af->apply_edit(\$sequence,$start);
print $sequence, "\n";
Description: Applies this AlleleFeature directly to a sequence which is
passed by reference.
If either the start or end of this AlleleFeature are not defined
this function will not do anything to the passed sequence.
Returntype : reference to the same sequence that was passed in
Exceptions : none
Caller : Slice
Status : At Risk |
Args : none
Example : my @sources = @{$af->get_all_sources()};
Description : returns a list of all the sources for this
AlleleFeature
ReturnType : reference to list of strings
Exceptions : none
Caller : general
Status : At Risk
: Variation database is under development. |
Arg [1] : (optional) Bio::EnsEMBL::Variation::Individual $individual
Example : $p = $af->individual();
Description: Getter/Setter for the individual associated with this feature.
If not set, and this AlleleFeature has an associated adaptor
an attempt will be made to lazy-load the individual from the
database.
Returntype : Bio::EnsEMBL::Variation::Individual
Exceptions : throw on incorrect argument
Caller : general
Status : At Risk |
Arg [1] : none
Example : my $diff = $af->length_diff();
Description: Returns the difference in length caused by applying this
AlleleFeature to a sequence. This may be be negative (deletion),
positive (insertion) or 0 (replacement).
If either start or end are not defined 0 is returned.
Returntype : int
Exceptions : none
Caller : general
Status : At Risk |
Arg [-ADAPTOR] :
see superclass constructor
Arg [-START] :
see superclass constructor
Arg [-END] :
see superclass constructor
Arg [-STRAND] :
see superclass constructor
Arg [-SLICE] :
see superclass constructor
Arg [-VARIATION_NAME] :
string - the name of the variation this feature is for (denormalisation
from Variation object).
Arg [-SOURCE] :
string - the name of the source where the SNP comes from
Arg [-VARIATION] :
int - the variation object associated with this feature.
Arg [-VARIATION_ID] :
int - the internal id of the variation object associated with this
identifier. This may be provided instead of a variation object so that
the variation may be lazy-loaded from the database on demand.
Arg [-SAMPLE_ID] :
int - the internal id of the sample object associated with this
identifier. This may be provided instead of the object so that
the population/individual may be lazy-loaded from the database on demand.
Arg [-ALLELE_STRING] :
string - the allele for this AlleleFeature object.
Example :
$af = Bio::EnsEMBL::Variation::AlleleFeature->new
(-start => 100,
-end => 100,
-strand => 1,
-slice => $slice,
-allele_string => 'A',
-variation_name => 'rs635421',
-source => 'Celera',
-sample_id => $sample_id,
-variation => $v);
Description: Constructor. Instantiates a new AlleleFeature object.
Returntype : Bio::EnsEMBL::Variation::AlleleFeature
Exceptions : none
Caller : general
Status : At Risk |
Args : None
Example : $allele = $obj->ref_allele_string()
Description: Getter for the reference allele.
Returntype : string
Exceptions : none
Caller : general
Status : At Risk |
Arg [1] : string $source (optional)
The new value to set the source attribute to
Example : $source = $vf->source()
Description: Getter/Setter for the source attribute
Returntype : string
Exceptions : none
Caller : general
Status : At Risk |
Arg [1] : (optional) Bio::EnsEMBL::Variation::Variation $variation
Example : $v = $af->variation();
Description: Getter/Setter for the variation associated with this feature.
If not set, and this AlleleFeature has an associated adaptor
an attempt will be made to lazy-load the variation from the
database.
Returntype : Bio::EnsEMBL::Variation::Variation
Exceptions : throw on incorrect argument
Caller : general
Status : At Risk |
Arg [1] : string $newval (optional)
The new value to set the variation_name attribute to
Example : $variation_name = $obj->variation_name()
Description: Getter/Setter for the variation_name attribute. This is the
name of the variation associated with this feature.
Returntype : string
Exceptions : none
Caller : general
Status : At Risk |
Methods code
sub allele_string
{ my $self = shift;
return $self->{'allele_string'} = shift if(@_);
return $self->{'allele_string'} if ($self->{'_half_genotype'});
return join('|',split (//,unambiguity_code($self->{'allele_string'})));
} | sub apply_edit
{
my $self = shift;
my $seqref = shift;
if(ref($seqref) ne 'SCALAR') {
throw("Reference to scalar argument expected");
}
if(!defined($self->{'start'}) || !defined($self->{'end'})) {
return $seqref;
}
my $len = $self->length;
substr($$seqref, $self->{'start'}-1, $len) = $self->{'allele_string'} if ($self->{'allele_string'} ne '-');
substr($$seqref, $self->{'start'}-1, 0) = $self->{'allele_string'} if ($self->{'allele_string'} eq '-');
return $seqref;} | sub get_all_sources
{ my $self = shift;
my @sources;
my %sources;
if ($self->{'adaptor'}){
map {$sources{$_}++} @{$self->{'adaptor'}->get_all_synonym_sources($self)};
$sources{$self->source}++;
@sources = keys %sources;
return\@ sources;
}
return\@ sources;
}
1;} | sub individual
{ my $self = shift;
if(@_) {
if(!ref($_[0]) || !$_[0]->isa('Bio::EnsEMBL::Variation::Individual')) {
throw("Bio::EnsEMBL::Variation::Individual argument expected");
}
$self->{'individual'} = shift;
}
elsif(!defined($self->{'individual'}) && $self->{'adaptor'} &&
defined($self->{'_sample_id'})) {
my $ia = $self->{'adaptor'}->db()->get_IndividualAdaptor();
$self->{'individual'} = $ia->fetch_by_dbID($self->{'_sample_id'});
if (!defined $self->{'individual'}){
warning("AlleleFeature attached to Strain, not Individual");
}
}
return $self->{'individual'};} | sub length_diff
{
my $self = shift;
return 0 if(!defined($self->{'end'}) || !defined($self->{'start'}));
return length($self->{'allele_string'}) - ($self->{'end'} - $self->{'start'} + 1) if ($self->{'allele_string'} ne '-');
return 0 - ($self->{'end'} - $self->{'start'} +1) if ($self->{'allele_string'} eq '-');
} | sub new
{ my $caller = shift;
my $class = ref($caller) || $caller;
my $self = $class->SUPER::new(@_);
my ($allele, $var_name, $variation, $variation_id,$population, $sample_id, $source) =
rearrange([qw(ALLELE_STRING VARIATION_NAME
VARIATION VARIATION_ID SAMPLE_ID SOURCE)], @_);
$self->{'allele_string'} = $allele;
$self->{'variation_name'} = $var_name;
$self->{'variation'} = $variation;
$self->{'_variation_id'} = $variation_id;
$self->{'_sample_id'} = $sample_id;
$self->{'source'} = $source;
return $self;} | sub new_fast
{ my $class = shift;
my $hashref = shift;
return bless $hashref, $class; } | sub ref_allele_string
{ my $self = shift;
my $reference_allele;
if ( ref ($self->slice) eq 'Bio::EnsEMBL::Slice' ){
$reference_allele = $self->seq;
}
else{
$reference_allele = $self->slice->ref_subseq($self->start,$self->end,$self->strand);
}
return $reference_allele;} | sub source
{ my $self = shift;
return $self->{'source'} = shift if(@_);
return $self->{'source'};} | sub variation
{ my $self = shift;
if(@_) {
if(!ref($_[0]) || !$_[0]->isa('Bio::EnsEMBL::Variation::Variation')) {
throw("Bio::EnsEMBL::Variation::Variation argument expected");
}
$self->{'variation'} = shift;
}
elsif(!defined($self->{'variation'}) && $self->{'adaptor'} &&
defined($self->{'_variation_id'})) {
my $va = $self->{'adaptor'}->db()->get_VariationAdaptor();
$self->{'variation'} = $va->fetch_by_dbID($self->{'_variation_id'});
}
return $self->{'variation'};} | sub variation_name
{ my $self = shift;
return $self->{'variation_name'} = shift if(@_);
return $self->{'variation_name'};} | General documentation