Commit e2a7af24 authored by Jan Hoeckesfeld's avatar Jan Hoeckesfeld
Browse files


parent 1bab1ebf
...@@ -81,14 +81,14 @@ generate_six_venndiagrams(type_kmer_counts, random_spatypes, snakemake.output['v ...@@ -81,14 +81,14 @@ generate_six_venndiagrams(type_kmer_counts, random_spatypes, snakemake.output['v
top = [] top = []
with open(snakemake.input['scores'], 'r') as f: with open(snakemake.input['scores'], 'r') as f:
for i, l in enumerate(f.readlines()): for i, l in enumerate(f.readlines()):
if i > 5: if i > 2:
break break
t = re.split(r'\t+', l.rstrip('\t\n')) t = re.split(r'\t+', l.rstrip('\t\n'))
if t[1] == "None": if t[1] == "None":
break break
top.append(t[0]) top.append(t[0])
U = set.union(*[set(expected_counts[spaTypeID].keys()) for spaTypeID in top]) U = set.intersection(*[set(expected_counts[spaTypeID].keys()) for spaTypeID in top])
outF = open(snakemake.output['topUkmers'], "w") outF = open(snakemake.output['topUkmers'], "w")
for line in U: for line in U:
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