aboutsummaryrefslogtreecommitdiff
path: root/below_melting_temperature.py
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2014-01-17 13:42:09 -0500
committerCalvin Morrison <mutantturkey@gmail.com>2014-01-17 13:42:09 -0500
commit96fdab29cf1cf500246d6ac84f1334bf7fa93bd4 (patch)
tree570ce4f6e06254e73be3657e16333f415d7e25c6 /below_melting_temperature.py
parentc5442da2baef615119e65f10e874886be452958c (diff)
add melting_range c code, update python, update makefile
Diffstat (limited to 'below_melting_temperature.py')
-rw-r--r--below_melting_temperature.py23
1 files changed, 8 insertions, 15 deletions
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