Fix bug that resulted in the TSL getting incorrectly evaluated when d…

…etermining the best peptide
griffithlab · Oct 31, 2023 · a9b7f38 · a9b7f38
1 parent 60dcec2
commit a9b7f38
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pvactools/lib/aggregate_all_epitopes.py b/pvactools/lib/aggregate_all_epitopes.py
@@ -334,7 +334,7 @@ def get_best_binder(self, df):
         #subset protein_coding dataframe to only include entries with a TSL < maximum_transcript_support_level
         tsl_df = biotype_df[biotype_df['Transcript Support Level'] != 'NA']
         tsl_df = tsl_df[tsl_df['Transcript Support Level'] != 'Not Supported']
-        tsl_df = tsl_df[tsl_df['Transcript Support Level'] < self.maximum_transcript_support_level]
+        tsl_df = tsl_df[tsl_df['Transcript Support Level'] <= self.maximum_transcript_support_level]
         #if this results in an empty dataframe, reset to previous dataframe
         if tsl_df.shape[0] == 0:
             tsl_df = biotype_df