When forging a BSgenome I'm getting this error :
> "Error in .readSeedFile(x, verbose = verbose) : seed file > 'SARSCOV2-seed' must have exactly 1 record"
This is the code I'm running:
> library(BSgenome) myfile <- read.dcf("seed_file.txt", fields = NULL, > all = FALSE, keep.white = NULL) write.dcf(myfile , file = "seed.dcf", > append = FALSE, useBytes = FALSE, indent = 0.1 * getOption("width"), > width = 0.9 * getOption("width"),keep.white = NULL) > forgeBSgenomeDataPkg("seed.dcf")
Please can you tell me what I'm doing wrong?
Thanks a lot.