import itertools


def get_genotype_calls_combinations(allele_frequencies):
  """
  given a list of allele frequencies, 
  return all possible genotype call combinations
  for example, if len(allele_frequencies) = 6, one combination may be [0,5]
  """
  genotypes = [i for i in range(len(allele_frequencies))]
  genotype_combinations = list(itertools.product(genotypes, genotypes)) 
  genotype_calls = [[i[0], i[1]] for i in genotype_combinations]
  return genotype_calls