Question: Ensembl VEP 89 Output Format
0
gravatar for bruce.moran
2.5 years ago by
bruce.moran670
Ireland
bruce.moran670 wrote:

I have just updated to VEP v89 from v81. Previously the output was an annotated VCF, but now it is a strange, new format. Can someone advise on how to get the former?

Thanks,

Bruce.

command:

VCF="input.vcf"
OUT=$(echo $VCF | sed 's/vcf/VEP89.vcf/')
STATS=$(echo $VCF | sed 's/vcf/VEP89.stats/')

${HOME}/ensembl-vep/vep \
    --cache \
    --offline \
    --merged \
    --fork 12 \
    --gencode_basic \
    --everything \
    --total_length \
    --allele_number \
    --no_escape \
    --check_existing \
    --xref_refseq \
    --variant_class \
    --dir ${HOME}/.vep \
    --fasta ${HOME}/human_g1k_v37.fasta \
    --assembly GRCh37 \
    --force_overwrite \
    --input_file $VCF \
    --output_file $OUT \
    --stats_file $STATS \
    --stats_text
vep ensembl • 809 views
ADD COMMENTlink modified 2.5 years ago by Emily_Ensembl20k • written 2.5 years ago by bruce.moran670
4
gravatar for Emily_Ensembl
2.5 years ago by
Emily_Ensembl20k
EMBL-EBI
Emily_Ensembl20k wrote:

--vcf will give you output in VCF

ADD COMMENTlink written 2.5 years ago by Emily_Ensembl20k

Thanks Emily. Out of curiosity, was that the default previously? I had just updated old scripts which never had that flag.

ADD REPLYlink written 2.5 years ago by bruce.moran670

I don't think it was, but I could be wrong.

ADD REPLYlink written 2.5 years ago by Emily_Ensembl20k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2264 users visited in the last hour