<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello All, <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">As part of a pipeline, I’ve been using the following code snippet to grab COSMIC variant IDs from specific regions of the human genome – this all worked wonderfully until the beginning of February this year and I’m wondering if the code
base has changed with the most recent ensembl release? Any advice on how to do this now, bearing in mind I am definitely a bit of a Bioinformatics amateur?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">my $slice_adaptor = $registry->get_adaptor('Human', 'Core', 'Slice');<o:p></o:p></p>
<p class="MsoNormal">my $slice = $slice_adaptor->fetch_by_region($coords);<o:p></o:p></p>
<p class="MsoNormal">my $vs = $slice->get_all_somatic_VariationFeatures();<o:p></o:p></p>
<p class="MsoNormal">my @variants = @{$vs};<o:p></o:p></p>
<p class="MsoNormal">foreach (@variants){<o:p></o:p></p>
<p class="MsoNormal"> my $name = $_->name();<o:p></o:p></p>
<p class="MsoNormal"> push (@variant_list, $name);<o:p></o:p></p>
<p class="MsoNormal">};<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">When this was working the end result was pushing of COSV IDs into the final @variant_list array which I was then able to print against regions in a final report. This array now ends up empty so either the “->name()” parameter or something
to do with COSMIC data storage amongst the “get_all_somatic_VariationFeatures” dataset may have been altered…<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Jon Williams <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>