# HG changeset patch # User drosofff # Date 1403514875 14400 # Node ID 59dfb33ca70e6f6c8e6fe6f7f0c930be25193ebf # Parent ac8fd41b8c98e44f13cd3338faeeedccf84d7748 Uploaded diff -r ac8fd41b8c98 -r 59dfb33ca70e smRtools.py --- a/smRtools.py Sun Jun 22 18:34:29 2014 -0400 +++ b/smRtools.py Mon Jun 23 05:14:35 2014 -0400 @@ -304,6 +304,8 @@ for offset in self.readDict: for size in self.readDict[offset]: dicsize[size] = dicsize.get(size, 0) + 1 + for offset in range (min(dicsize.keys()), max(dicsize.keys())+1): + dicsize[size] = dicsize.get(size, 0) # to fill offsets with null values return dicsize def statsizes (self, upstream_coord=None, downstream_coord=None): @@ -393,7 +395,7 @@ return mylist def readcoverage (self, upstream_coord=None, downstream_coord=None, windowName=None): - '''This method has not been tested yet 15-11-2013''' + '''Use by MirParser tool''' upstream_coord = upstream_coord or 1 downstream_coord = downstream_coord or self.size windowName = windowName or "%s_%s_%s" % (self.gene, upstream_coord, downstream_coord)