cpg.py 文件源码

def __call__( self, block ):
        if not block: return block
        if len(block.components) < 2:
            return
        cpglist = list_cpg_restricted( \
            string.upper(block.components[0].text), \
            string.upper(block.components[1].text) )

        # now we have a fast list of CpG columns, iterate/mask
        self.masked += len(cpglist)
        self.total += len(block.components[0].text)
        for component in block.components:
            component.text = mask_columns( cpglist, component.text, self.mask )

        return block

# Inclusive. Mask out all sites that are not non-CpG sites.