This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Convert list of rs IDs from one format to another

Hello, I downloaded a huge list of eQTLs from the GTEx portal:

https://www.gtexportal.org/home/datasets

There is a column with the rs IDs with the format: "1_115746_C_T_b37". I would like to convert them to rs IDs of the form "rs1234567"? (these are just examples). In the documentation it says these IDs are RS IDs from dbSNP 147

Thanks !

snp identifiers eqtl gtex

Can you give context on the IDs? How did you get them? What's the biological context?

Hi Hussain Ather, I just edited my post

1 answer

Hello Mr Locuace,

this "ID" looks to me like it describes:

  • the chromosome
  • position
  • reference allele
  • alternative allele
  • reference genome

So to get the corresponding rs id, one solution is to extract the informations of the id and convert it to a vcf file. You can then annotate the ID column. There was a similar thread some time ago.

Let's start to create a vcf file. I assume you have a file id.txt which contain your id's in every line:

awk -F_ -v OFS="\t" '{print $1,$2,".",$3,$4,".",".","."}' id.txt|sort -k1,1V -k2,2g > ids.vcf

Now you can take this vcf file and annotate it with e.g. SnpSift:

java -jar SnpSift.jar annotate -id dbSNP.vcf.gz  ids.vcf > ids_annotated.vcf

Of course first you have to download the dbSnp file.

fin swimmer

Vielen Dank finswimmer !! ;)

dbSNP release must be based on b37/GRCh37 build.

Log in to answer this question.