Fix previous extraneous gRNA conv call, to give targets in output file

H2muller · Aug 1, 2024 · 35a4c7a · 35a4c7a
1 parent 4c08769
commit 35a4c7a
Show file tree

Hide file tree

Showing 2 changed files with 17,317 additions and 17,315 deletions.
diff --git a/CROPSR.py b/CROPSR.py
@@ -458,7 +458,9 @@ def main():
 
                     lesser_list = Complete_dataset[index_range:index_range+count]
 
-                    sequences = [np.frombuffer(bytes(str(get_gRNA_sequence(item[6]).replace('U','T')).upper(),"ascii"), 'uint8') if len(item[6]) == 30
+                    seqs = [str(get_gRNA_sequence(get_reverse_complement(item[6])).replace('U','T')).upper() for item in lesser_list]
+
+                    sequences = [np.frombuffer(bytes(str(get_gRNA_sequence(get_reverse_complement(item[6])).replace('U','T')).upper(),"ascii"), 'uint8') if len(item[6]) == 30
                     else np.empty(30,) for item in lesser_list ]
 
                     score = rs1_score(np.array(sequences))