# 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)