--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/reactivity_cal/separate_rna.py	Mon Oct 20 14:32:13 2014 -0400
@@ -0,0 +1,42 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+import sys
+from parse_dis_pac import *
+
+
+dist_file = sys.argv[1]
+cdna_file = sys.argv[2]
+rrna_file = sys.argv[3]
+
+dist = parse_dist(dist_file)
+dist = dist[1]
+hc = file(cdna_file, 'w')
+hr = file(rrna_file, 'w')
+
+for t in dist:
+    if t.find('AT') != -1:
+        hc.write(t)
+        hc.write('\n')
+        for i in range(len(dist[t])-1):
+            hc.write(dist[t][i])
+            hc.write('\t')
+        i = i+1
+        hc.write(dist[t][i])
+        hc.write('\n')
+    else:
+        hr.write(t)
+        hr.write('\n')
+        for i in range(len(dist[t])-1):
+            hr.write(dist[t][i])
+            hr.write('\t')
+        i = i+1
+        hr.write(dist[t][i])
+        hr.write('\n')
+
+hc.close()
+hr.close()
+
+    
+        
+