Mercurial > repos > tyty > structurefold
comparison predict/predict_RNAs.py @ 69:bd02fcf26f7a draft
Uploaded wrapper that correctly handles structure prediction without constraints.
| author | devteam |
|---|---|
| date | Mon, 24 Nov 2014 18:21:30 -0500 |
| parents | fb80870002a3 |
| children |
comparison
equal
deleted
inserted
replaced
| 68:fb80870002a3 | 69:bd02fcf26f7a |
|---|---|
| 53 fh.write('>'+id_s) | 53 fh.write('>'+id_s) |
| 54 fh.write('\n') | 54 fh.write('\n') |
| 55 fh.write(seqs[id_s]) | 55 fh.write(seqs[id_s]) |
| 56 fh.close() | 56 fh.close() |
| 57 if not flag: | 57 if not flag: |
| 58 command = shlex.split('Fold %s %s' % (os.path.join(syspath, temp.txt), os.path.join(output_directory, '%s.ct' % id_s))) | 58 command = shlex.split('Fold %s %s' % (os.path.join(syspath, 'temp.txt'), os.path.join(output_directory, '%s.ct' % id_s))) |
| 59 subprocess.call(command) | 59 subprocess.call(command) |
| 60 else: | 60 else: |
| 61 if id_s in react: | 61 if id_s in react: |
| 62 fh = file(os.path.join(syspath, "constraint.txt"), 'w') | 62 fh = file(os.path.join(syspath, "constraint.txt"), 'w') |
| 63 make_plot(react[id_s], id_s, output_directory) #make a plot of the distribution of the reactivites of the input RNA | 63 make_plot(react[id_s], id_s, output_directory) #make a plot of the distribution of the reactivites of the input RNA |