Mercurial > repos > fubar > microsatbed

diff find_str.py @ 16:ee773bbde3c3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/fubar2/microsatbed commit d952bc313f408735456747c3d33e09a3170c8f59-dirty
author fubar
date Wed, 17 Jul 2024 05:08:00 +0000
parents 1085e094cf5f
children 264d79548d19
line wrap: on
line diff
--- a/find_str.py	Wed Jul 17 04:28:28 2024 +0000
+++ b/find_str.py	Wed Jul 17 05:08:00 2024 +0000
@@ -27,12 +27,12 @@
             ssrs = pytrf.STRFinder(
                 name,
                 seq,
-                args.minreps,
-                args.minreps,
-                args.minreps,
-                args.minreps,
-                args.minreps,
-                args.minreps,
+                2,
+                1,
+                2,
+                2,
+                2,
+                2,
             )
         else:
             ssrs = pytrf.STRFinder(
@@ -55,7 +55,7 @@
                 ssr.length,
             )
             # pytrf reports a 1 based start position so start-1 fixes the bed interval lengths
-            if args.specific and ssr.motif in specific:
+            if args.specific and ssr.motif in specific and ssr.repeat >= args.minreps:
                 bed.append(row)
             elif args.mono and len(ssr.motif) == 1:
                 bed.append(row)