Mercurial > repos > tyty > structurefold
diff reactivity_cal/separate_rna.py @ 5:7a8ddf1819b1 draft
Uploaded
author | tyty |
---|---|
date | Mon, 15 Sep 2014 14:52:52 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/reactivity_cal/separate_rna.py Mon Sep 15 14:52:52 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() + + + +