Commit 5d3374aa authored by Mikaël Salson's avatar Mikaël Salson Committed by Vidjil

sample_set.py: genescan, 100-600bp, more robust

parent 75bfde63
......@@ -703,13 +703,16 @@ def getFusedStats(file_name, res, dest):
tmp = {}
for c in data['clones']:
try:
arl = int(math.ceil(c['_average_read_length'][result_index]))
except:
continue
if arl > 0:
if arl not in tmp:
tmp[arl] = 0.0
tmp[arl] += float(c['reads'][result_index])
min_len = 0 #int(min(tmp.keys()))
max_len = 500 #int(max(tmp.keys()))
min_len = 100 #int(min(tmp.keys()))
max_len = 600 #int(max(tmp.keys()))
tmp_list = []
if mapped_reads == 0:
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment