Closed:Not running the the code for compare2Sequences(
4
0
Entering edit mode
@susobhanbaidya-13780
Last seen 5.3 years ago

I have installed the  required packages properly. Now I crate two similar sequence by name X1.fa and X2.fa .

And try to run the following code. But getting uncanny error.

Please help me 

 

> library(CRISPRseek)
> outputDir <- getwd()
> inputFile1Path <- "X1.fa"
> inputFile2Path <- "X2.fa" 
> REpatternFile <- system.file("extdata", "NEBenzymes.fa", package = "CRISPRseek")
> seqs <- compare2Sequences(inputFile1Path, inputFile2Path, outputDir = outputDir, 
+ REpatternFile = REpatternFile, overwrite = TRUE)
search for gRNAs for input file1...
Validating input ...
Searching for gRNAs ...
<simpleError in findgRNAs(inputFilePath, findPairedgRNAOnly = findPairedgRNAOnly,     annotatePaired = annotatePaired, pairOutputFile = pairOutputFile,     PAM = PAM, PAM.location = PAM.location, gRNA.pattern = gRNA.pattern,     PAM.size = PAM.size, gRNA.size = gRNA.size, min.gap = min.gap,     max.gap = max.gap, name.prefix = gRNA.name.prefix, format = format,     featureWeightMatrixFile = featureWeightMatrixFile, baseBeforegRNA = baseBeforegRNA,     baseAfterPAM = baseAfterPAM, calculategRNAEfficacy = TRUE,     efficacyFile = efficacyFile): inputfile specified as X1.fa does not exists!>
search for gRNAs for input file2...
Validating input ...
Searching for gRNAs ...
<simpleError in findgRNAs(inputFilePath, findPairedgRNAOnly = findPairedgRNAOnly,     annotatePaired = annotatePaired, pairOutputFile = pairOutputFile,     PAM = PAM, PAM.location = PAM.location, gRNA.pattern = gRNA.pattern,     PAM.size = PAM.size, gRNA.size = gRNA.size, min.gap = min.gap,     max.gap = max.gap, name.prefix = gRNA.name.prefix, format = format,     featureWeightMatrixFile = featureWeightMatrixFile, baseBeforegRNA = baseBeforegRNA,     baseAfterPAM = baseAfterPAM, calculategRNAEfficacy = TRUE,     efficacyFile = efficacyFile): inputfile specified as X2.fa does not exists!>
[1] "Scoring ..."
Error in .Call2("new_input_filexp", filepath, PACKAGE = "XVector") : 
  cannot open file 'X1.fa'

 

software error CRISPRseek • 133 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 445 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6