I am trying to recode alleles for my .bim file in PLINK , for example from A,B allele coding to A,C,G,T.
My file
1 rs101 0 4566795 A B
1 rs102 0 4640902 B A
PLINK requires that you use the --update-alleles command .
/plink --bfile mydata --update-alleles mylist.txt --make-bed --out newfile
with which the file mylist.txt contains five columns per row listing For example,
rs101 A B G T
rs102 A B A C
My question is when creating the mylist.txt file, for the last two columns, How do you determine the new allele code for first allele and the new allele code for other allele ( A is no G, B is now T)?
•
written
by
wondereye •
You're the one who wanted to recode them. Put any letters in there you please! Maybe the first A/B can be A/T, and the second can become that novel nucleotide X/Y.
Plink uses AB for a reason, that is linkage screening uses binary markers and doesn't care about the sequence. Since most of the high frequency SNPs will be intergenic, it's not going to have amino acid changes anyway.
If you for some reason want to have the real sequence, you need to find out where the file came from in the first place, because each microarray technology will use a semi-random strand (AT vs TA at a het site). If dbSNP said that rs101 is A/G, we can't tell if your data file is referring to G/A or A/G, or even reverse strand T/C and C/T.
That's something that will be annotated by the microarray manufacturer, or maybe available through bioconductor and Arrayexpress annotations. Maybe NCBI GEO will have platform information definitions too. I did this a few years ago with Affymetrix, and was forced to make an account with them to get to the reference data about my particular chip type.
•
written
by
karl.stamm
Thank you so much @karl.stamm. I'm in just learning how to use PLINK. I was hoping to find out if there were any base changes for the SNP but, I'm assuming I can look this up on dbSNP?
Log in to answer this question.
