I have ~ 400 text files of BLASTp results.
NP_387917.1 ADZ06570.1 44.29 289 153 4 2 289 1 282 7e-77 236
I have a file that reads in files ALL my P_1 files, sorts the e values, and then calculates a threshold. I would like to convert this into a function. I have made some attempts to define a function that can read in all my files, but as of now, the files only reads ALL P_1 files, or ALL P_2 files. I also feel like the function is not sorting the values correctly.
def threshold_proteinsVStarget(protein_Files): # List of Lactobacillus databases created from BLAST commands Lactobacillus_DB = [ 'L_GRL1112', 'L_214','L_CTV-05','L_JV-V01','L_ST1','L_MV-1A', 'L_202-4','L_224-1','L_JV-V03', 'L_MV-22','L_DSM_13335', 'LactinV_03V1b', 'SPIN_1401G', 'UPII_143_D', 'L_1153','L_269-3', 'L_JV-V16','L_49540'] list_e_values =  for prot in protein_Files: for db in Lactobacillus_DB: line = open (prot + db + '.out').readline() print line line2 = line.strip('\n') fields = line2.split('\t') e_val = float(fields) list_e_values.append(e_val) return list_e_values
file1 = ['P_3']