From 65e9bc6772f4538078d0eb1d78ff7b08df6d9c4e Mon Sep 17 00:00:00 2001
From: Calvin Morrison <mutantturkey@gmail.com>
Date: Mon, 24 Mar 2014 15:36:53 -0400
Subject: use biopython's temperature calculator... slow as dirt

---
 src/below_melting_temperature.py | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

(limited to 'src')

diff --git a/src/below_melting_temperature.py b/src/below_melting_temperature.py
index e2521db..e91e777 100644
--- a/src/below_melting_temperature.py
+++ b/src/below_melting_temperature.py
@@ -1,5 +1,8 @@
+#!/usr/bin/env python
 import sys
+from Bio.SeqUtils.MeltingTemp import Tm_staluc
 
+# naiive
 def in_temp_range(kmer):
 
 	A = kmer.count('A')
@@ -22,5 +25,5 @@ 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
+	if min_melting_temp < Tm_staluc(line.split("\t")[0]) < max_melting_temp:
+		sys.stdout.write(line)
-- 
cgit v1.2.3