Skip to content
Snippets Groups Projects
Commit f0656809 authored by Daniel Rios's avatar Daniel Rios
Browse files

modified to use get_all_AlleleFeatures_Slice instead of get_all_differences_Slice

parent 21318ce8
No related branches found
No related tags found
No related merge requests found
......@@ -29,7 +29,7 @@ Bio::EnsEMBL::AlignStrainSlice - Represents the slice of the genome aligned with
-STRAINS => \@strainSlices);
#get coordinates of variation in alignSlice
my $alleleFeatures = $strainSlice1->get_all_differences_Slice();
my $alleleFeatures = $strainSlice1->get_all_AlleleFeature_Slice();
foreach my $af (@{$alleleFeatures}){
my $new_feature = $alignSlice->alignFeature($af, $strainSlice1);
print "Coordinates of the feature in AlignSlice are: ", $new_feature->start, "-", $new_feature->end, "\n";
......@@ -210,7 +210,7 @@ sub _get_indels{
#go throuh all the strains getting ONLY the indels (length_diff <> 0)
my @indels;
foreach my $strainSlice (@{$self->strains}){
my $differences = $strainSlice->get_all_differences_Slice(); #need to check there are differences....
my $differences = $strainSlice->get_all_AlleleFeatures_Slice(); #need to check there are differences....
if (defined $differences){
my @results = grep {$_->length_diff != 0} @{$differences};
push @indels, @results;
......
......@@ -24,10 +24,10 @@ Bio::EnsEMBL::StrainSlice - SubClass of the Slice. Represents the slice of the g
my $seq = $strainSlice->seq();
print $seq;
#get differences between this StrainSlice and the reference
my $differences = $strainSlice->get_all_differences_Slice();
foreach my $difference (@{$differences}){
print "Difference in position ", $difference->start,"-",$difference->end," in strain with allele ", $difference->allele_string,"\n";
#get allele features between this StrainSlice and the reference
my $afs = $strainSlice->get_all_AlleleFeatures_Slice();
foreach my $af (@{$afs}){
print "AlleleFeature in position ", $af->start,"-",$af->end," in strain with allele ", $af->allele_string,"\n";
}
#compare a strain against another strain
my $strainSlice_2 = $slice->get_by_strain($strain_name_2);
......
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