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))