index = index + 1
continue
else:
new_sequence = new_sequence + sequence[index]
index = index + 1
natural_out_sequences.write(new_sequence + "\n")
natural_out_sequences.close()
#list of files that contain L vs RSA data
data = []
#list of sequences with missing RSA values
bad_list = []
#Get all the RSA values using DSSP
seq_data = af.get_values(pdb_id, chain_id)
RSAValues = []
RSA = seq_data[1]
#print "RSA Length from dssp script: " + str(len(RSA))
index = 0
#Open the files you just created with the gapless sequences and read in the sequences.
fp_natural = open(natural_pdb_file_title)
natural_data_arr = fp_natural.readlines()
fp_natural.close()
fp_designed = open(duncan_designed_sequence_path + designed_pdb_file_title)
designed_data_arr = fp_designed.readlines()
fp_designed.close()
#filepointer for the partial result
index = index + 1
continue
else:
new_sequence = new_sequence + sequence[index]
index = index + 1
natural_out_sequences.write(new_sequence + "\n")
natural_out_sequences.close()
#list of files that contain L vs RSA data
data = []
#list of sequences with missing RSA values
bad_list = []
#Get all the RSA values using DSSP
seq_data = af.get_values(pdb_id, chain_id)
RSAValues = []
RSA = seq_data[1]
print "RSA Length from dssp script: " + str(len(RSA))
index = 0
#Open the files you just created with the gapless sequences
fp_natural = open(natural_pdb_file_title)
natural_data_arr = fp_natural.readlines()
fp_natural.close()
#filepointer for the partial result
fpW_natural = open("results_array_natural_" + pdb_id + "_" + chain_id + ".csv","w")
#write the RSA values:
fpW_natural.write(af.dump_csv_line(RSA))