motif.py 文件源码

def mapping(self):
        bases = ['A', 'T', 'C', 'G']
        motifs = {}
        for i in range(6):
            for motif in itertools.product(bases, repeat=i+1):
                motif = "".join(list(motif))
                if not is_motif(motif):
                    continue

                smotif = self.standard(motif)
                if smotif not in motifs:
                    motifs[smotif] = []

                if motif not in motifs[smotif]:
                    motifs[smotif].append(motif)

        return motifs