@@ -83,8 +83,13 @@ getAzimuthScores <- function(sequences, fork=FALSE){

             sequences.valid <- rep(sequences.valid, 2)

         }

         .dumpToFile(sequences.valid, inputfile)

-        system2("python",

-                c(programFile, inputfile, outputfile))

+

+        pyBinary <- basilisk.utils:::getPythonBinary(env)

+

+        system2(c(pyBinary,

+                  programFile,

+                  inputfile,

+                  outputfile))

         scores <- read.table(outputfile)[,1]

         if (sum(good)==1){

             scores <- scores[1]

@@ -179,8 +179,14 @@ getCrispraiScores <- function(tss_df,

                            "getWeissmanScores.py",

                            package="crisprScore",

                            mustWork=TRUE)

-    system2("python",

-            c(program, rosterFile, modality, outputFile, verbose))

+

+    pyBinary <- basilisk.utils:::getPythonBinary(env)

+    system2(c(pyBinary,

+              program,

+              rosterFile,

+              modality,

+              outputFile,

+              verbose))

     scores <- read.table(outputFile)

     scores <- as.data.frame(scores)

@@ -89,8 +89,13 @@ getDeepSpCas9Scores <- function(sequences,

     #              outputfile)

     if (sum(good)>0){

         .dumpToFile(sequences.valid, inputfile)

-        system2("python",

-                c(programFile, inputfile, modelDir, outputfile))

+        pyBinary <- basilisk.utils:::getPythonBinary(env)

+

+        system2(c(pyBinary,

+                  programFile,

+                  inputfile,

+                  modelDir,

+                  outputfile))

         scores <- read.table(outputfile)[,1]

         #scores <- read.table(outputfile)

         scores <- scores/100