comparison find_str.py @ 16:ee773bbde3c3 draft

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
comparison
equal deleted inserted replaced
15:4af0486b8513 16:ee773bbde3c3
25 for name, seq in fa: 25 for name, seq in fa:
26 if args.specific: 26 if args.specific:
27 ssrs = pytrf.STRFinder( 27 ssrs = pytrf.STRFinder(
28 name, 28 name,
29 seq, 29 seq,
30 args.minreps, 30 2,
31 args.minreps, 31 1,
32 args.minreps, 32 2,
33 args.minreps, 33 2,
34 args.minreps, 34 2,
35 args.minreps, 35 2,
36 ) 36 )
37 else: 37 else:
38 ssrs = pytrf.STRFinder( 38 ssrs = pytrf.STRFinder(
39 name, 39 name,
40 seq, 40 seq,
53 ssr.motif, 53 ssr.motif,
54 ssr.repeat, 54 ssr.repeat,
55 ssr.length, 55 ssr.length,
56 ) 56 )
57 # pytrf reports a 1 based start position so start-1 fixes the bed interval lengths 57 # pytrf reports a 1 based start position so start-1 fixes the bed interval lengths
58 if args.specific and ssr.motif in specific: 58 if args.specific and ssr.motif in specific and ssr.repeat >= args.minreps:
59 bed.append(row) 59 bed.append(row)
60 elif args.mono and len(ssr.motif) == 1: 60 elif args.mono and len(ssr.motif) == 1:
61 bed.append(row) 61 bed.append(row)
62 elif args.di and len(ssr.motif) == 2: 62 elif args.di and len(ssr.motif) == 2:
63 bed.append(row) 63 bed.append(row)