Mercurial > repos > abims-sbr > pairwise

diff scripts/S03_run_second_blast.py @ 10:cb5f2cc3dedc draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/abims-sbr/adaptsearch commit 7d45b8e07bf99420df9b195617515b8002b5eaf3
author abims-sbr
date Fri, 06 Jul 2018 02:51:09 -0400
parents 1e0c547d88fe
children
line wrap: on
line diff
--- a/scripts/S03_run_second_blast.py	Tue Jul 03 10:52:31 2018 -0400
+++ b/scripts/S03_run_second_blast.py	Fri Jul 06 02:51:09 2018 -0400
@@ -21,7 +21,7 @@
         translated_file = open(f_name, 'w')
         with open(args.query_file, 'r') as file:
             for name, seq in itertools.izip_longest(*[file]*2):
-                s = Seq(seq.strip('\n').upper(), IUPAC.unambiguous_dna)
+                s = Seq(seq.strip('\n').upper(), IUPAC.ambiguous_dna)
                 translated_file.write(name.strip('\n')+'_orf_1\n')
                 translated_file.write(s.translate()._data+'\n')
                 translated_file.write(name.strip('\n')+'_orf_2\n')
@@ -42,4 +42,4 @@
         exit()
 
 if __name__ == "__main__":
-    main()
\ No newline at end of file
+    main()