From 5df2b9950b1748148b90c33b7a004056491f9e51 Mon Sep 17 00:00:00 2001 From: Calvin Morrison Date: Fri, 17 Jan 2014 14:32:06 -0500 Subject: move melting temp py script to src --- below_melting_temperature.py | 26 -------------------------- 1 file changed, 26 deletions(-) delete mode 100644 below_melting_temperature.py (limited to 'below_melting_temperature.py') 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 -- cgit v1.2.3