I have with me now a group of about 20 SNPs that I would like to acquire all possible genotype combinations for. As an example, let's start off with three SNPs and their alleles.
SNP A1 A2 SNP1 A T SNP2 C G SNP3 T A
I want to start off by generating a list of all possible genotype permutations/combinations of these three SNPs, for example:
SNP1 SNP2 SNP3 AA CC TT AA CC TA AA CC AA AA CG TT AA CG TA AA CG AA AA GG TT AA GG TA AA GG AA ...
And so on, for what I expect to be 3^3 = 27 possible combinations.
From here, I hope to scale this up to my full group of ~20 SNPs. What is a good way of doing this, in Python or even in R?