INNER CODE UNIT · Python
cluster
wheaton5/souporcell · consensus.py:303
cluster = cell_clusters[cell-1]
for (locus, counts) in loci_counts.items():
ref = 0
alt = 0
if counts[0] > 0:
ref = counts[0]
if counts[1] > 0:
alt = counts[1]
if locus-1 in loci_for_soup:
cluster_allele_counts_soup[loci_for_soup[locus-1]][cluster][0] += ref
cluster_allele_counts_soup[loci_for_soup[locus-1]][cluster][1] += alt
average_allele_expression_soup[loci_for_soup[locus-1]][0] += ref
average_allele_expression_soup[loci_for_soup[locus-1]][1] += alt
if not locus-1 in locus_index:
continue
cluster_allele_counts[locus_index[locus-1]][cluster][0] += ref
cluster_allele_counts[locus_index[locus-1]][cluster][1] += alt
average_allele_expression[locus_index[locus-1]][0] += ref