Get ccs_seqname from the a aligned reference rather than constructing…

… it from the subread. PiperOrigin-RevId: 512104765
google · Feb 24, 2023 · d3c94f0 · d3c94f0
1 parent 10114e6
commit d3c94f0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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)