INNER CODE UNIT · Python
expanded_alleles
BirolLab/straglr · straglr_compare.py:175
expanded_alleles = set.intersection(*map(set, [vc[locus][0] for vc in vs_controls]))
if not expanded_alleles:
continue
supports = []
for allele in expanded_alleles:
for vc in vs_controls:
if allele in vc[locus][3]:
supports.append(vc[locus][3][allele])
break
control_alleles = []
pvals = []
control_alleles_list = []
for vc in vs_controls:
if vc[locus][1]:
control_alleles.append(vc[locus][1][0])
for allele in vc[locus][1][0].split(','):