[ensembl-dev] VEP run with condel, sift and polyphen options

binit treesa binit.treesa at gmail.com
Fri Oct 12 05:45:21 BST 2012


>
> Now that I have subscribed in Ensemble dev..
>
> It would be very helpful for me if you check this vep issue with the
> condel module.
>
> Thanks
> Treesa Binit
>
>
> On Wed, Oct 10, 2012 at 4:29 PM, <dev-owner at ensembl.org> wrote:
>
>> The Ensembl dev mailing list only accepts postings from people who are
>> subscribed. You can subscribe or unsubscribe at
>> http://lists.ensembl.org/mailman/listinfo/dev
>>
>>
>>
>> ---------- Forwarded message ----------
>> From: binit treesa <binit.treesa at gmail.com>
>> To: dev at ensembl.org
>> Cc: Neeba Dijo <neebasebastian at gmail.com>
>> Date: Wed, 10 Oct 2012 16:29:35 +0530
>> Subject: VEP run with condel, sift and polyphen options
>> Hi,
>>
>> I have downloaded the VEP version 2.6 and tried to run this script with a
>> variant file. It executed succesfully and generated the output file .
>> But when I executed this with the options for "sift, polyphen and condel
>> plugin", it throws some erros. Eventhough it calculated all the score and
>> generated the output file.
>>
>> And I'm confused on the correctness of these calculated score.
>>
>> It would be better if you check and comment about the output generated
>> here.
>>
>> Details of the vep.pl run is given below
>>
>> *Input File :test.txt*
>>
>> #CHROM_ID       START   ID      REF     ALT
>> 1       94586536        .       C       A
>> 1       94586536        .       C       G
>>
>>
>> *Output file : testOutput.txt*
>>
>> ## ENSEMBL VARIANT EFFECT PREDICTOR v2.6
>> ## Output produced at 2012-10-10 06:37:22
>> ## Connected to homo_sapiens_core_68_37 on ensembldb.ensembl.org
>> ## Using API version 68, DB version 68
>> ## Extra column keys:
>> ## DISTANCE : Shortest distance from variant to transcript
>> ## PolyPhen : PolyPhen prediction
>> ## SIFT : SIFT prediction
>> ## Condel    : Consensus deleteriousness score for an amino acid
>> substitution based on SIFT and PolyPhen-2
>> #Uploaded_variation    Location    Allele    Gene    Feature
>> Feature_type    Consequence    cDNA_position    CDS_position
>> Protein_position    Amino_acids    Codons    Existing_variation    Extra
>> 1_94586536_C/A    1:94586536    A    ENSG00000198691
>> ENST00000370225    Transcript    missense_variant,splice_region_variant
>> 153    66    22    K/N    aaG/aaT    -
>> PolyPhen=possibly_damaging(0.859);Condel=deleterious(0.885);SIFT=tolerated(0.09)
>> 1_94586536_C/A    1:94586536    A    ENSG00000198691
>> ENST00000535735    Transcript    missense_variant,splice_region_variant
>> 146    66    22    K/N    aaG/aaT    -
>> PolyPhen=benign(0.235);Condel=deleterious(0.583);SIFT=tolerated(0.07)
>> 1_94586536_C/G    1:94586536    G    ENSG00000198691
>> ENST00000370225    Transcript    missense_variant,splice_region_variant
>> 153    66    22    K/N    aaG/aaC    -
>> PolyPhen=possibly_damaging(0.859);Condel=deleterious(0.885);SIFT=tolerated(0.09)
>> 1_94586536_C/G    1:94586536    G    ENSG00000198691
>> ENST00000535735    Transcript    missense_variant,splice_region_variant
>> 146    66    22    K/N    aaG/aaC    -
>> PolyPhen=benign(0.235);Condel=deleterious(0.583);SIFT=tolerated(0.07)
>>
>>
>> *Command used*
>>
>> perl variant_effect_predictor.pl -i test.txt -o testOutput.txt -sift b
>> -polyphen b  -plugin
>> Condel,/users/tbinit/.vep/Plugins/config/Condel/config,b --force
>>
>> *The status of the processing is shown *
>>
>> readline() on closed filehandle SIFT at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 136.
>> readline() on closed filehandle POLYPHEN at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 147.
>> 2012-10-10 05:50:00 - Loaded plugin: Condel
>> 2012-10-10 05:50:02 - Starting...
>> 2012-10-10 05:50:02 - Detected format of input file as vcf
>> Use of uninitialized value $data[4] in string eq at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 316, <GEN0> line 4.
>> Use of uninitialized value $ref in length at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 349, <GEN0> line 4.
>> Use of uninitialized value $alt in pattern match (m//) at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 352, <GEN0> line 4.
>> Use of uninitialized value $alt in split at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 425, <GEN0> line 4.
>> Use of uninitialized value $alt in pattern match (m//) at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 434, <GEN0> line 4.
>> Use of uninitialized value $ref in concatenation (.) or string at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 501, <GEN0> line 4.
>> Use of uninitialized value $alt in concatenation (.) or string at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 501, <GEN0> line 4.
>> Use of uninitialized value $data[2] in string eq at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 501, <GEN0> line 4.
>> Use of uninitialized value in pattern match (m//) at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 1944, <GEN0> line 4.
>> Use of uninitialized value in substitution (s///) at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 1944, <GEN0> line 4.
>> Use of uninitialized value in string eq at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 1945, <GEN0> line 4.
>> Use of uninitialized value in pattern match (m//) at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 1950, <GEN0> line 4.
>> Use of uninitialized value in concatenation (.) or string at
>> /depts/pd/refDataManagment/vep_2.6/variant_effect_predictor/Bio/EnsEMBL/Variation/Utils/VEP.pm
>> line 1951, <GEN0> line 4.
>> WARNING: Start  or end -1 coordinate invalid on line 4
>> 2012-10-10 05:50:02 - Read 2 variants into buffer
>> 2012-10-10 05:50:02 - Skipping 0 non-variant loci
>> 2012-10-10 05:50:02 - Reading transcript data from cache and/or database
>> [================================================================================================================================]
>> [ 100% ]
>> 2012-10-10 05:50:06 - Retrieved 9 transcripts (0 mem, 0 cached, 9 DB, 0
>> duplicates)
>> 2012-10-10 05:50:06 - Analyzing chromosome 1
>> 2012-10-10 05:50:06 - Analyzing variants
>> [================================================================================================================================]
>> [ 100% ]
>> 2012-10-10 05:50:08 - Calculating consequences
>> [>
>> ]    [ 0% ]Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 203, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 212, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 203, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 216, <GEN0> line 4.
>> [==============================================================>
>> ]   [ 50% ]Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 203, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 212, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 203, <GEN0> line 4.
>> Use of uninitialized value in subtraction (-) at
>> /users/tbinit/.vep/Plugins/config/Condel/Condel.pm line 216, <GEN0> line 4.
>> [================================================================================================================================]
>> [ 100% ]
>> 2012-10-10 05:50:11 - Processed 2 total variants (0 vars/sec, 0 vars/sec
>> total)
>> 2012-10-10 05:50:11 - Finished!
>>
>>
>> Thanks
>> Treesa Binit
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20121012/bfa6607f/attachment.html>


More information about the Dev mailing list