diff options
Diffstat (limited to 'below_melting_temperature.py')
-rw-r--r-- | below_melting_temperature.py | 26 |
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 |