From 39a653d1b49c0e4fa5de1783e7aba1743fd66e32 Mon Sep 17 00:00:00 2001 From: Calvin Morrison Date: Thu, 13 Mar 2014 17:07:27 -0400 Subject: add coutning points for reverses --- src/score_mers.py | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/score_mers.py b/src/score_mers.py index 9ea9d4a..c10b04a 100755 --- a/src/score_mers.py +++ b/src/score_mers.py @@ -105,6 +105,10 @@ def populate_locations(input_fn, mers, mer): for line in strstream.stdout: mers[mer].pts.append(int(line)) + cmd = 'tac ' + input_fn + " | tr '[ACGT]' '[TGCA]' | strstreamone " + mer + " " + input_fn + strstream = Popen(cmd, stdout=PIPE, shell=True) + for line in strstream.stdout: + mers[mer].pts.append(int(line)) def apply_filters(combination): for mer in combination: -- cgit v1.2.3