[ensembl-dev] VEP (ver78) TranscriptVariationAllele.pm unitialized value error with --hgvs enabled

Sarah Hunt seh at ebi.ac.uk
Thu Feb 5 12:21:41 GMT 2015


Hi Oliver,

Thanks for the detailed report. This was a similar issue with a stop 
inserted before the stop not being handled well. It is fixed in branch78 
so the consequence of the variant has changed from 'stop_gained' to 
'stop_retained'. The HGVS string reflects the synonymous status.

Best wishes,

Sarah


On 03/02/2015 16:52, Homann, Oliver wrote:
>
> Apologies, forgot to clarify that the error occurs in line 1221 of the 
> new version of TranscriptVariationAllele.pm:
>
>    1221 unless( $hgvs_notation->{ref} eq "-"){
>
>    1222 $hgvs_notation->{ref}  = 
> Bio::SeqUtils->seq3(Bio::PrimarySeq->new(-seq => 
> $hgvs_notation->{ref}, -id => 'ref',  -alphabet => 'protein')) || "";
>
>    1223   }
>
> Regards,
>
> Oliver
>
> *From:*dev-bounces at ensembl.org [mailto:dev-bounces at ensembl.org] *On 
> Behalf Of *Homann, Oliver
> *Sent:* Tuesday, February 03, 2015 8:50 AM
> *To:* Ensembl developers list
> *Subject:* Re: [ensembl-dev] VEP (ver78) TranscriptVariationAllele.pm 
> unitialized value error with --hgvs enabled
>
> Hello Sarah,
>
> Your fix worked for all but one of my ‘uninitialized value’ errors.  
> Would you mind trying to track down the issue with this remaining 
> problematic entry?
>
> #CHROM POS        ID            REF         ALT         QUAL FILTER   
> INFO      FORMAT              SAMPLENAME
>
> 19 42355761             rs374945441       C CTAGAAA            
> 885.10   . 
> AC=1;AF=5.556e-03;AN=180;BaseQRankSum=-1.922e+00;DB;DP=2670;FS=2.561;GQ_MEAN=75.28;GQ_STDDEV=91.16;InbreedingCoeff=-0.0056;MLEAC=1;MLEAF=5.556e-03;MQ=60.00;MQ0=0;MQRankSum=-2.160e-01;NCC=0;QD=9.83;ReadPosRankSum=0.353;SOR=0.669 
> GT:AD:DP:GQ:PL              0/0:17,0:17:34:0,34,505
>
> Thanks,
>
> Oliver
>
>
>
> _______________________________________________
> Dev mailing list    Dev at ensembl.org
> Posting guidelines and subscribe/unsubscribe info: http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20150205/8df60d42/attachment.html>


More information about the Dev mailing list