def _dict_seq_locus(list_c, loci_obj, seq_obj):
seqs = defaultdict(set)<EOL>for c in list_c.values():<EOL><INDENT>for l in c.loci2seq:<EOL><INDENT>[seqs[s].add(c.id) for s in c.loci2seq[l]]<EOL><DEDENT><DEDENT>common = [s for s in seqs if len(seqs[s]) > <NUM_LIT:1>]<EOL>seqs_in_c = defaultdict(float)<EOL>for c in list_c.values():<EOL><INDENT>for l in c.loci2seq:<EOL><INDENT>for s in c.loci2seq[l]:<EOL><INDENT>if s in common:<EOL><INDENT>pos = seq_obj[s].pos[l]<EOL>cov = <NUM_LIT:1.0> * loci_obj[l].coverage[pos]<EOL>if seqs_in_c[(s, c.id)] < cov:<EOL><INDENT>seqs_in_c[(s, c.id)] = cov<EOL><DEDENT><DEDENT><DEDENT><DEDENT><DEDENT>seqs_in_c = _transform(seqs_in_c)<EOL>return seqs_in_c<EOL>
return dict with sequences = [ cluster1, cluster2 ...]