+ resolved an issue with the memory mapping in filterReads
[qpalma.git] / scripts / compile_dataset.py
index c4f3b88..aec5fc5 100644 (file)
@@ -191,12 +191,15 @@ def process_read(reReads,fRead,currentGene,dna_flat_files,test):
 
    # if we perform testing we cut a much wider region as we want to detect
    # how good we perform on a region
+   import random
+   ADD_POS = random.randint(1,20)
+
    if test:
-      up_cut = up_cut-500
+      up_cut = up_cut-ADD_POS
       if up_cut < 0:
          up_cut = 0
 
-      down_cut = down_cut+500
+      down_cut = down_cut+ADD_POS
 
       if down_cut > len(currentSeq):
          down_cut = len(currentSeq)