<div dir="ltr">Since the only use of a string is $hgvs_notation->{ref}, a guess is to change line 1220 to: <br><br> unless( defined($hgvs_notation->{ref}) && $hgvs_notation->{ref} eq "-" )<br><br><div>This masks the problem of why this didn't get defined in the first place which might be a more serious matter.<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 29, 2015 at 4:47 PM, Homann, Oliver <span dir="ltr"><<a href="mailto:ohomann@amgen.com" target="_blank">ohomann@amgen.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have been using the VEP to annotate a GATK-generated VCF containing indels from a WGS genotyping project, and have encountered a small number of entries that result in the following error (only when –hgvs is enabled):<u></u><u></u></p>
<p class="MsoNormal"> “Use of uninitialized value in string eq at /<mypath>/ensembl-tools-release-78/scripts/variant_effect_predictor/Bio/EnsEMBL/Variation/TranscriptVariationAllele.pm line 1220”<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Here are the offending lines in the perl module:<u></u><u></u></p>
<p class="MsoNormal"> 1219 ### Convert peptide to 3 letter code as used in HGVS<u></u><u></u></p>
<p class="MsoNormal"> 1220 unless( $hgvs_notation->{ref} eq "-"){<u></u><u></u></p>
<p class="MsoNormal"> 1221 $hgvs_notation->{ref} = Bio::SeqUtils->seq3(Bio::PrimarySeq->new(-seq => $hgvs_notation->{ref}, -id => 'ref', -alphabet => 'protein')) || "";<u></u><u></u></p>
<p class="MsoNormal"> 1222 }<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I tracked down one of the offending entries, but it isn’t immediately obvious to me why it is causing an error:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT MySample<u></u><u></u></p>
<p class="MsoNormal">13 21729289 . T TCAGTTTTCTTTGTTGCTGACATCTCGGATGTTCTGTCCATGTTTAAGGAACCTTTTA 118.19 . AC=1;AF=5.556e-03;AN=180;BaseQRankSum=0.744;DP=2652;FS=3.538;GQ_MEAN=66.68;GQ_STDDEV=20.39;InbreedingCoeff=-0.0085;MLEAC=1;MLEAF=5.556e-03;MQ=60.00;MQ0=0;MQRankSum=0.609;NCC=0;QD=0.13;ReadPosRankSum=-1.220e-01;SOR=0.252
GT:AD:DP:GQ:PL 0/1:41,8:.:99:160,0,14558<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Any ideas what the problem might be?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks,<u></u><u></u></p>
<p class="MsoNormal">Oliver Homann<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">P.S. Apart from this small glitch, VEP is teriffic!<u></u><u></u></p>
</div>
</div>
<br>_______________________________________________<br>
Dev mailing list <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
<br></blockquote></div><br></div>