fixed bed strand if none provided (#9)

jonas-fuchs · Nov 3, 2023 · ed3a24d · ed3a24d
1 parent 7247533
commit ed3a24d
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/bamdash/scripts/data.py b/bamdash/scripts/data.py
@@ -269,6 +269,8 @@ def bed_to_dict(bed_file, coverage_df, ref, min_cov):
         bed_dict["bed annotations"][f"{start} {stop}"] = {}
         bed_dict["bed annotations"][f"{start} {stop}"]["start"] = start
         bed_dict["bed annotations"][f"{start} {stop}"]["stop"] = stop
+        # always add strand dummy information in case its not in bed file
+        bed_dict["bed annotations"][f"{start} {stop}"]["strand"] = "NA"
         # check for additional info
         if len(line) > 3:
             for element, classifier in zip(line[3:], possible_classifiers):