#! /usr/bin/python2.7

import os,subprocess,shlex

if __name__ == "__main__":
	inFile=open("PDBsequences.csv",'r')
	outFile=open("entropy_300.csv",'w')
	outFile.write("pdb\tEntropy\n")
	for line in inFile.readlines():
		rna = line.strip().split('\t')[1]
		pdbId = line.strip().split('\t')[0]
		#print rna
		cmd = "/home/amir/entropy/RNAentropy/RNAentropy_1.8.5/RNAentropy %s" %rna
		ent=subprocess.Popen(shlex.split(cmd),stdout=subprocess.PIPE)
		outEnt,errEnt = ent.communicate()
		outFile.write("%s\t%s\n"%(pdbId,outEnt.strip()))
	inFile.close()
	outFile.close()