From 96fdab29cf1cf500246d6ac84f1334bf7fa93bd4 Mon Sep 17 00:00:00 2001 From: Calvin Morrison Date: Fri, 17 Jan 2014 13:42:09 -0500 Subject: add melting_range c code, update python, update makefile --- below_melting_temperature.py | 23 ++++++++--------------- 1 file changed, 8 insertions(+), 15 deletions(-) (limited to 'below_melting_temperature.py') diff --git a/below_melting_temperature.py b/below_melting_temperature.py index bc4157e..e2521db 100644 --- a/below_melting_temperature.py +++ b/below_melting_temperature.py @@ -1,8 +1,6 @@ import sys -from multiprocessing import Pool -def above_melting_temperature(kmer_with_count): - kmer = kmer_with_count.split("\t")[0] +def in_temp_range(kmer): A = kmer.count('A') C = kmer.count('C') @@ -16,18 +14,13 @@ def above_melting_temperature(kmer_with_count): else: melt_temp = 64.9 + 41*(G+C-16.4)/(A+T+G+C) - if melt_temp < max_melting_temp: - return kmer_with_count - else: - return 0 - + return min_melting_temp < melt_temp < max_melting_temp -lines = sys.stdin.readlines() -max_melting_temp = float(sys.argv[1]) +min_melting_temp = float(sys.argv[1]) +max_melting_temp = float(sys.argv[2]) -p = Pool() -output = p.map(above_melting_temperature, lines) -for line in output: - if line != 0: - sys.stdout.write(line) +output = [] +for line in sys.stdin: + if in_temp_range(line.split("\t")[0]): + sys.stdout.write(line) \ No newline at end of file -- cgit v1.2.3