[ensembl-dev] VEP (ver78) TranscriptVariationAllele.pm unitialized value error with --hgvs enabled
Homann, Oliver
ohomann at amgen.com
Tue Feb 3 16:52:56 GMT 2015
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20150203/32aa6b4a/attachment.html>
More information about the Dev
mailing list