Ensembl / ensembl-vep

The Ensembl Variant Effect Predictor predicts the functional effects of genomic variants
https://www.ensembl.org/vep
Apache License 2.0
456 stars 152 forks source link

`0` does not work as a variant identifier #1714

Closed dvg-p4 closed 5 months ago

dvg-p4 commented 5 months ago

Describe the issue

If a provided variant identifier is 0, it will be ignored and the default identifier will be used instead.

Additional information

System

Full VEP command line

~/ensembl-vep/vep --input_file repro_input.tsv --format ensembl --cache --offline --dir ~/.vep --assembly GRCh38 --output_file repro_output.vcf --vcf

Data files

repro_input.tsv:

1   1220794 1220794 G/T +   -2
1   1220795 1220795 C/T +   -1
1   1220796 1220796 G/A +   0
1   1220796 1220796 G/C +   1
1   1220797 1220797 G/A +   2
1   1220801 1220801 C/G +   3

$ grep -v "^##" repro_output.vcf | cut -f 1-5 | column -t:

#CHROM  POS      ID             REF  ALT
1       1220794  -2             G    T
1       1220795  -1             C    T
1       1220796  1_1220796_G/A  G    A
1       1220796  1              G    C
1       1220797  2              G    A
1       1220801  3              C    G
olaaustine commented 5 months ago

Hi @dvg-p4, Hope this meets you well? Thank you for your query. To better understand your issue, do you want to use 0 as a variant identifier ? Thank you Ola.

dvg-p4 commented 5 months ago

Yes, I'd like for "0" to be usable as a variant identifier.

olaaustine commented 5 months ago

Hi @dvg-p4. Thank you for your response. At this time, we do not support this feature. However, I will discuss the possibility of adding this functionality with our team for future updates. I will close this ticket now, please feel free to open another ticket or reopen this ticket if you have any issues. Thank you Ola