Reverse quality string in force_merge_pairs

pinellolab · Mar 22, 2023 · b074699 · b074699
1 parent 587a56e
commit b074699
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/CRISPResso2/CRISPRessoShared.py b/CRISPResso2/CRISPRessoShared.py
@@ -1191,7 +1191,7 @@ def force_merge_pairs(r1_filename, r2_filename, output_filename):
         id2 = f2.readline()
         seq2 = reverse_complement(f2.readline().strip()) + "\n"
         plus2 = f2.readline()
-        qual2 = f2.readline()
+        qual2 = f2.readline().strip()[::-1]
 
         f_out.write(id1 + seq1 + seq2 + plus1 + qual1 + qual2)