[ensembl-dev] Variation set: fetch_all_by_Variation() outputting hash values

Gong, Henry Henry.Gong at ucsf.edu
Tue Jul 22 23:42:09 BST 2014


Hi Anja,

It looks like I had assumed that simply calling VariationSet would produce a name, similar to
my @vstates = @{$tva->variation_feature->get_all_validation_states()};
print
        "\t",join(",", @vstates); #evidence
Though I see now that an evidence value is qualitatively different from a variation set.
I don't intend on calling all the variations per set; my application is only looking at potentially hundreds (per transcript or per gene), but I'll consider using iterators to save a little memory.

I added your suggestion as
    foreach my $tva(@{$tvas}) {
        my @vsets = @{$vs_adaptor->fetch_all_by_Variation($tva->variation_feature->variation)};
        print $tva->variation_feature->display_id, "\t";
        foreach my $vs(@vsets) {
            my $set_name = $vs->name();
            print $set_name,",";
        }
        print "\n";
    }
and it does indeed work!

Thanks,
Henry Gong
Junior Specialist
UCSF Dept of Neuro Surgery
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20140722/677c68d3/attachment.html>


More information about the Dev mailing list