Quickfix for polyT_adpator_finder

shimlab · Apr 29, 2024 · 7cd4fc2 · 7cd4fc2
1 parent cdc518e
commit 7cd4fc2
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/blaze/polyT_adaptor_finder.py b/blaze/polyT_adaptor_finder.py
@@ -112,7 +112,11 @@ def find_poly_T(seq, poly_T_len=PLY_T_LEN,
             T_prop = helper.sliding_window_mean(read_code, poly_T_len)
             return np.where(T_prop >= min_match_prop)[0]  
 
+        strand = self._strand if not strand else strand
+        read = self.seq if not read else read
+
         if strand == '-':
+            seq = read[:num_nt]
             # find polyT
             ply_T_idx = find_poly_T(seq)
             d1, d2 = SEQ_SUFFIX_AFT_ADPT