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.