From ed3a24dbebf0091a2171c169a47962484b5e7d04 Mon Sep 17 00:00:00 2001 From: Jonas Fuchs <78491186+jonas-fuchs@users.noreply.github.com> Date: Fri, 3 Nov 2023 15:57:46 +0100 Subject: [PATCH] fixed bed strand if none provided (#9) --- bamdash/scripts/data.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/bamdash/scripts/data.py b/bamdash/scripts/data.py index 6d6d917..054f059 100644 --- 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):