From d3c94f083b153e1c7816a3ca848a780f44531680 Mon Sep 17 00:00:00 2001 From: danielecook Date: Fri, 24 Feb 2023 10:29:51 -0800 Subject: [PATCH] Get ccs_seqname from the a aligned reference rather than constructing it from the subread. PiperOrigin-RevId: 512104765 --- deepconsensus/preprocess/pre_lib.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/deepconsensus/preprocess/pre_lib.py b/deepconsensus/preprocess/pre_lib.py index ea9829b..a0d3403 100644 --- a/deepconsensus/preprocess/pre_lib.py +++ b/deepconsensus/preprocess/pre_lib.py @@ -1315,7 +1315,7 @@ def proc_feeder(): counter=main_counter, ) subreads = list(map(expand_clip_indent_count, read_set)) - ccs_seqname = '/'.join(subreads[0].name.split('/')[:2] + ['ccs']) + ccs_seqname = read_set[0].reference_name # Fetch ccs sequence and append to subread set. while True: ccs_bam_read = next(ccs_bam_h)