# HG changeset patch # User fubar # Date 1721258352 0 # Node ID 45f690db0eaf7ca4ac4973638496ca71d97d2458 # Parent b5acf60b903d794fe1713c405d8d5bcdda681dfb planemo upload for repository https://github.com/fubar2/microsatbed commit d952bc313f408735456747c3d33e09a3170c8f59-dirty diff -r b5acf60b903d -r 45f690db0eaf dibed.bw Binary file dibed.bw has changed diff -r b5acf60b903d -r 45f690db0eaf find_str.py --- a/find_str.py Wed Jul 17 22:41:41 2024 +0000 +++ b/find_str.py Wed Jul 17 23:19:12 2024 +0000 @@ -11,14 +11,14 @@ Designed to build some of the microsatellite tracks from https://github.com/arangrhie/T2T-Polish/tree/master/pattern for the VGP. """ -def getDensity(name, bed, len, winwidth): - nwin = int(len / winwidth) +def getDensity(name, bed, chrlen, winwidth): + nwin = int(chrlen/winwidth) d = [0.0 for x in range(nwin+1)] for b in bed: nt = b[5] bin = int(b[1]/winwidth) d[bin] += nt - dw = [(name,x*winwidth,(x+1)*winwidth,float(d[x])) for x in range(nwin+1) if (x+1)*winwidth <= len] + dw = [(name,x*winwidth,(x+1)*winwidth,float(d[x])) for x in range(nwin+1) if (x+1)*winwidth <= chrlen] return dw def write_ssrs(args): diff -r b5acf60b903d -r 45f690db0eaf microsatbed.xml --- a/microsatbed.xml Wed Jul 17 22:41:41 2024 +0000 +++ b/microsatbed.xml Wed Jul 17 23:19:12 2024 +0000 @@ -92,6 +92,7 @@ help="Bed features will be output for every motif of the selected length(s) with the minimum required repeats or more"> + @@ -102,6 +103,7 @@ help="Bed features will be output for every motif of the selected length(s) with the minimum required repeats or more"> +