This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ensembl data

Hi, I used the Ensembl website for SNPs annotation. I uploaded to their website a VCF file with a list of SNPs. Example of SNPs:

#Chrom  Pos Ref Alt 
1   13116   T   G   
1   13118   A   G   
1   16682   G   A   
1   900161  C   G   
1   902288  G   A   
1   980460  G   A   
1   1362903 G   C   
1   1414714 A   G
1   1420704 C   T   
1   1560103 C   T
1   1600156 C   G   

And I got thos files as a result : enter image description here

The problem is that the result file from the website doesnt have the REF and ALT column so i cant know exactly which SNP it stands for example : 1:13116-13116 is chrom 1 pos 13116 but no information in REF and ALT allele to connect it to my SNP file as there are multiple entries for 1:13116-13116. Is there any way to "fix" this output in Ensembl? is it connected to :"Transcript database to use" on ensembl (https://www.ensembl.org/Homo_sapiens/Tools/VEP?db=core;expand_form=true;tl=YSu3e94WHDBwtREC-8879777) Thank you :)

ensembl vcf snp

Do not leave your previous questions unresolved.

1 answer

You live, you learn. Always add a unique identifier in the input so you can map the output back. I add an ID column with ${CHR}_${POS}_${REF}_${ALT} as the values for this purpose.

Add this column and re-run VEP.

the problem is that the ID column in my VCF from where I got the SNPs is NA for all SNPs so I identifty SNPs based on CHROM-POR-REF-ALT combination, so is there any other way ?

No, you're going to need to re-run it. Prepare the input before providing it to VEP.

Log in to answer this question.