diff --git a/RNA_Seq_Cecere/RNA-seq.snakefile b/RNA_Seq_Cecere/RNA-seq.snakefile
index 03798aa83efc0bcf60306cf33fa7cf1af9f57d50..fd6421f59e7cd3fd2c727cbd3cd8da3876da0d7a 100644
--- a/RNA_Seq_Cecere/RNA-seq.snakefile
+++ b/RNA_Seq_Cecere/RNA-seq.snakefile
@@ -1367,10 +1367,10 @@ def do_linear_regression(data, x_col, y_col, y_min=0, fit_intercept=True, transf
     """
     if transform is not None:
         data = getattr(np, transform)(
-            data[[x_col, y_col]].dropna().query(f"{y_col} >= {y_min}")).replace(
+            data[[x_col, y_col]].dropna().query(f"`{y_col}` >= {y_min}")).replace(
                 [np.inf, -np.inf], np.nan).dropna()
     else:
-        data = data[[x_col, y_col]].dropna().query(f"{y_col} >= {y_min}")
+        data = data[[x_col, y_col]].dropna().query(f"`{y_col}` >= {y_min}")
     lm = LinearRegression(fit_intercept=fit_intercept)
     X = data.drop(y_col, axis=1)
     lm.fit(X, data[y_col])