def test_duplicate_subreads(self): """Test that report doesn't crash on duplicate subreads""" output_dir = tempfile.mkdtemp(suffix="_mapping_stats") with warnings.catch_warnings(record=True) as w: report = to_report(self.bam_file, output_dir) shutil.rmtree(output_dir) self.assertTrue(w >= 4)
def setUpClass(cls): cls.output_dir = tempfile.mkdtemp(suffix="_mapping_stats") t = tempfile.NamedTemporaryFile( delete=False, suffix="mapping_report.json") t.close() cls.report_json = t.name cls.report = to_report(cls._get_input_file(), cls.output_dir) cls.report.write_json(cls.report_json) assert isinstance(cls.report, Report) log.info(pprint.pformat(cls.report.to_dict())) for table in cls.report.tables: log.info(str(table))
def setUpClass(cls): cls.output_dir = tempfile.mkdtemp(suffix="_mapping_stats") t = tempfile.NamedTemporaryFile(delete=False, suffix="mapping_report.json") t.close() cls.report_json = t.name cls.report = to_report(cls._get_input_file(), cls.output_dir) cls.report.write_json(cls.report_json) assert isinstance(cls.report, Report) log.info(pprint.pformat(cls.report.to_dict())) for table in cls.report.tables: log.info(str(table))
def setUpClass(cls): _to_p = lambda x: os.path.join(_GMAP_DATA_DIR, x) cls.output_dir = tempfile.mkdtemp(suffix="_mapping_stats") cls.input_fofn = _to_p('input.fofn') cls.aligned_reads_cmph5 = _to_p('aligned_reads.cmp.h5') t = tempfile.NamedTemporaryFile( delete=False, suffix="mapping_report.json") t.close() cls.report_json = t.name cls.report = to_report(cls.aligned_reads_cmph5, cls.output_dir) if isinstance(cls.report, Report): log.info(pprint.pformat(cls.report.to_dict())) for table in cls.report.tables: log.info(str(table))
def setUpClass(cls): _to_p = lambda x: op.join(_GMAP_DATA_DIR, x) cls.output_dir = tempfile.mkdtemp(suffix="_mapping_stats") cls.input_fofn = _to_p('input.fofn') cls.aligned_reads_cmph5 = _to_p('aligned_reads.cmp.h5') t = tempfile.NamedTemporaryFile(delete=False, suffix="mapping_report.json") t.close() cls.report_json = t.name cls.report = to_report(cls.aligned_reads_cmph5, cls.output_dir) if isinstance(cls.report, Report): log.info(pprint.pformat(cls.report.to_dict())) for table in cls.report.tables: log.info(str(table))