write regions into fasta headers

kriese · Jun 3, 2022 · 322143d · 322143d
1 parent db303d0
commit 322143d
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/workflow/scripts/ref_cut.py b/workflow/scripts/ref_cut.py
@@ -45,7 +45,7 @@ def collect_liftover(args):
 
     liftover_chain_files = {}
     seq_coll = SequenceCollection()
-    seq_coll.references.append(NamedSequence(main_ref_name, main_ref_region.get_reference_seq(main_reference)))
+    seq_coll.references.append(NamedSequence(f"{main_ref_name} {main_ref_region}", main_ref_region.get_reference_seq(main_reference)))
 
     # If we have no information about how to handle liftovers, we dont to them
     if not Path(args.liftover_paths).exists():
@@ -77,7 +77,7 @@ def collect_liftover(args):
                 print("Alternative reference region is absurdely large. Liftover probably failed.")
             else:
                 alt_ref_str = alt_ref_region.get_reference_seq(alt_ref)
-                seq_coll.references.append(NamedSequence(alt_ref_name, alt_ref_str))
+                seq_coll.references.append(NamedSequence(f"{alt_ref_name} {alt_ref_region}", alt_ref_str))
         else:
             print(f"Region not found in {alt_ref_names[ref_idx]} reference!")