fix: forget to lowercase the record tx_id when restranding

IGDRion · Jan 13, 2024 · 925f1b7 · 925f1b7
1 parent 3a39e48
commit 925f1b7
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/bin/filter_gtf_ndr.py b/bin/filter_gtf_ndr.py
@@ -122,8 +122,9 @@ def filter_count_matrix(file, transcripts, wr):
 
     with open("unformat.novel.filter.gtf", "w") as wr:
         for record in GTF.parse_by_line(args.gtf):
-            if "transcript_id" in record and record["transcript_id"].lower() in filter:
-                record.strand = strand_dict[record["transcript_id"]]
+            tx_id = record["transcript_id"].lower()
+            if "transcript_id" in record and tx_id in filter:
+                record.strand = strand_dict[tx_id]
                 print(record, file=wr)
 
     with open("counts_transcript.filter.txt", "w") as wr: