So I've recently downloaded a dataset of genomes associated with a study. It's in the Bed/Bim/Fam format, which is what PLINK uses. However, how do I convert the dataset in the Bed/Bim/Fam files to a 23andme format for each individual in the dataset?
What I've done so far is: --bfile test --extract 23andMe.txt --make-bed --out test, where 23andMe.txt contains only markers for 23andme. Which is what I want.
Now, what I would like, is obviously to output a file with the same format as 23andme for each individual.
23andme format looks like this:
rsid chromosome position genotype
rs4477212 1 82154 AA
rs3094315 1 752566 AA
rs3131972 1 752721 GG
rs12124819 1 776546 AA
rs11240777 1 798959 AG
Any help would be appreciated. Thank you.