report.py 文件源码

def _set_mapping_metrics(self, read_type, read_type_set):
        for genome, region in itertools.product(self.genomes, self.regions):
            is_read_type = (genome, region) in read_type_set
            reads_frac = self._get_metric_attr('reads_frac', genome, region, read_type)
            reads_frac.add(1, filter=is_read_type)

        for region in self.regions:
            is_read_type = any([(genome, region) in read_type_set for genome in self.genomes])
            multi_reads_frac = self._get_metric_attr('reads_frac', cr_constants.MULTI_REFS_PREFIX, region,
                                                     read_type)
            multi_reads_frac.add(1, filter=is_read_type)