view reactivity_cal/separate_rna.py @ 5:7a8ddf1819b1 draft

Uploaded

#!/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()