Hi, I am a student currently studying to obtain my master's degree and for my dissertation we are thinking about using Python in order to find all of the possible combinations of genotypes for a set of 24 SNPs. I am very unfamiliar with coding, so I have no idea where to start or how we should go about getting the outcome that we are looking for. Any help or advice would be greatly appreciated! :)

Our samples came from ancient DNA, so there was a low genome coverage to begin with. Therefore, a lot of the SNPs we are looking at (looking at the 24 SNPs required for use with HIrisPlex) were either missing or had low coverage. Here is what I have so far for one of our samples (possible number of alleles for each SNP);

```
SNP 1 - 0
SNP 2 - 0
SNP 3 - 1
SNP 4 - 0, 1
SNP 5 - 0, 1
SNP 6 - 0
SNP 7 - 0, 1
```

...

We need to find out all of the possible combinations for all 24 SNPs, so we can obtain a range of probabilities for each hair and eye color so that we can estimate the predicted phenotype for each individual. Any help or advice anyone has on how we would go about doing this would be greatly appreciated! Thank you in advance! :)

Amanda

I added markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below: