<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hello,
<div class=""><br class="">
</div>
<div class="">For the following variant</div>
<div class="">17<span class="Apple-tab-span" style="white-space:pre"> </span>7576915<span class="Apple-tab-span" style="white-space:pre">
</span>7577019<span class="Apple-tab-span" style="white-space:pre"> </span>TGTTGGGCAGTGCTAGGAAAGAGGCAAGGAAAGGTGATAAAAGTGAATCTGAGGCATAACTGCACCCTTGGTCTCCTCCACCGCTTCTTGTCCTGCTTGCTTACC/-<span class="Apple-tab-span" style="white-space:pre">
</span>+</div>
<div class=""><br class="">
</div>
<div class="">the correct protein consequence for the transcript ENST00000269305 should be </div>
<div class="">ENSP00000269305.4:p.Ala307ThrfsTer34</div>
<div class=""><br class="">
</div>
<div class="">As is confirmed by the VEP web interface, REST API and the VEP perl script.</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">However, when I try to get this using the Variation API, I get</div>
<div class="">ENSP00000269305.4:p.Ala307Thr</div>
<div class=""><br class="">
</div>
<div class="">This is how I’m using the API:</div>
<div class=""><br class="">
</div>
<div class=""><font face="Menlo" class="">$reg->load_registry_from_db(-host => '<a href="http://ensembldb.ensembl.org" class="">ensembldb.ensembl.org</a>', -species => 'homo sapiens', -user => 'anonymous', -port => 3337);</font></div>
<div class=""><font face="Menlo" class="">my $sa = $reg->get_adaptor(‘human’, 'core', 'slice');</font></div>
<div class=""><font face="Menlo" class="">my $vfa = $reg->get_adaptor('human', 'variation', 'VariationFeature');</font></div>
<div class=""><font face="Menlo" class="">my $slice = $sa->fetch_by_region(‘chromosome', 17);</font></div>
<div class=""><font face="Menlo" class="">my $vf = Bio::EnsEMBL::Variation::VariationFeature->new(-start => 7576915, -end => 7577019, -slice => $slice, </font></div>
<div class=""><font face="Menlo" class=""><span class="Apple-tab-span" style="white-space:pre"></span>-allele_string => ‘TGTTGGGCAGTGCTAGGAAAGAGGCAAGGAAAGGTGATAAAAGTGAATCTGAGGCATAACTGCACCCTTGGTCTCCTCCACCGCTTCTTGTCCTGCTTGCTTACC/-‘, -strand => 1, -adaptor =>
$vfa);<span class="Apple-tab-span" style="white-space:pre"> </span></font></div>
<div class=""><font face="Menlo" class="">my $cons = $vf->get_all_TranscriptVariations;</font></div>
<div class=""><font face="Menlo" class="">my @hgvsp = grep { defined } map { $_->hgvs_protein } map { @{$_->get_all_alternate_TranscriptVariationAlleles} } @$cons;</font></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I’m using v89 for the API and the VEP script.</div>
<div class=""><br class="">
</div>
<div class="">Is there something I’m doing wrong in my script?</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Bhavana</div>
<br class="">
<br class="">
<div class=""><br class="webkit-block-placeholder">
</div>
<div class="">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">
<br class="Apple-interchange-newline">
---------------------------------------------<br class="">
Bhavana Harsha<br class="">
Bioinformatician<br class="">
COSMIC<br class="">
Wellcome Trust Sanger Institute<br class="">
Hinxton, UK<br class="">
CB10 1SA<br class="">
</div>
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">
<br class="">
</div>
<br class="Apple-interchange-newline">
</div>
<br class="">
</body>
</html>