Mercurial > repos > tyty > structurefold
view get_reads/reactivity_cal/separate_rna.py @ 20:13024d65eee7 draft
Uploaded
author | tyty |
---|---|
date | Mon, 20 Oct 2014 14:31:27 -0400 |
parents | |
children |
line wrap: on
line source
#!/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()