aboutsummaryrefslogtreecommitdiff
path: root/below_melting_temperature.py
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2014-01-17 14:32:06 -0500
committerCalvin Morrison <mutantturkey@gmail.com>2014-01-17 14:32:06 -0500
commit5df2b9950b1748148b90c33b7a004056491f9e51 (patch)
tree74315ead07de8c49844943d1e1f04b9fe2720534 /below_melting_temperature.py
parent62d8086e777a7dd306eda23577003a113b32a8f0 (diff)
move melting temp py script to src
Diffstat (limited to 'below_melting_temperature.py')
-rw-r--r--below_melting_temperature.py26
1 files changed, 0 insertions, 26 deletions
diff --git a/below_melting_temperature.py b/below_melting_temperature.py
deleted file mode 100644
index e2521db..0000000
--- a/below_melting_temperature.py
+++ /dev/null
@@ -1,26 +0,0 @@
-import sys
-
-def in_temp_range(kmer):
-
- A = kmer.count('A')
- C = kmer.count('C')
- G = kmer.count('G')
- T = kmer.count('T')
-
- melt_temp = 0.0;
-
- if len(kmer) < 13:
- melt_temp = ((A+T) * 2) + ((C+G) * 4)
- else:
- melt_temp = 64.9 + 41*(G+C-16.4)/(A+T+G+C)
-
- return min_melting_temp < melt_temp < max_melting_temp
-
-min_melting_temp = float(sys.argv[1])
-max_melting_temp = float(sys.argv[2])
-
-
-output = []
-for line in sys.stdin:
- if in_temp_range(line.split("\t")[0]):
- sys.stdout.write(line) \ No newline at end of file