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 -------------------------- src/below_melting_temperature.py | 26 ++++++++++++++++++++++++++ 2 files changed, 26 insertions(+), 26 deletions(-) delete mode 100644 below_melting_temperature.py create mode 100644 src/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 diff --git a/src/below_melting_temperature.py b/src/below_melting_temperature.py new file mode 100644 index 0000000..e2521db --- /dev/null +++ b/src/below_melting_temperature.py @@ -0,0 +1,26 @@ +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