Hey everyone! I followed this tutorial (https://support.bioconductor.org/p/125514/) in order to create a package for Papio anubis genome from NCBI. I originally wanted to use Ensembl's assembly but BSgenome mainly takes UCSC and NCBI.
Here is my seed file:
Package: BSgenome.Panubis.NCBI.Panubis1.0
Title: Papio anubis full genome (NCBI Panubis1.0)
Description: Papio anubis genome assembly as provided by NCBI (accession GCA_008728515.1), submitted by University of California, San Francisco
Version: 1.0.0
organism: Papio anubis
common_name: Olive baboon
genome: Panubis1.0
provider: NCBI
provider_version: Panubis1.0
release_date: 2019/10/2
source_url: https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/008/728/515/GCF_008728515.1_Panubis1.0/GCF_008728515.1_Panubis1.0_genomic.fna.gz
organism_biocview: Papio_anubis
BSgenomeObjname: Panubis
SrcDataFiles: GCF_008728515.1_Panubis1.0_genomic.fna from https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/008/728/515/GCF_008728515.1_Panubis1.0/GCF_008728515.1_Panubis1.0_genomic.fna.gz
seqs_srcdir: /scratch/nguyenth/bcm/trial/Panubis1.0_og
seqnames: c(paste0("NC_0449",sprintf("%0d",76:97),".1"),("NC_020006.2"))
circ_seqs: "NC_020006.2"
This is what I ran and the result:
">forgeBSgenomeDataPkg("Panubis-seed")
Creating package in ./BSgenome.Panubis.NCBI.Panubis1.0 Loading 'NC_044976.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044976.1.fa' ... DONE Loading 'NC_044977.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044977.1.fa' ... DONE Loading 'NC_044978.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044978.1.fa' ... DONE Loading 'NC_044979.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044979.1.fa' ... DONE Loading 'NC_044980.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044980.1.fa' ... DONE Loading 'NC_044981.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044981.1.fa' ... DONE Loading 'NC_044982.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044982.1.fa' ... DONE Loading 'NC_044983.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044983.1.fa' ... DONE Loading 'NC_044984.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044984.1.fa' ... DONE Loading 'NC_044985.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044985.1.fa' ... DONE Loading 'NC_044986.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044986.1.fa' ... DONE Loading 'NC_044987.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044987.1.fa' ... DONE Loading 'NC_044988.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044988.1.fa' ... DONE Loading 'NC_044989.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044989.1.fa' ... DONE Loading 'NC_044990.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044990.1.fa' ... DONE Loading 'NC_044991.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044991.1.fa' ... DONE Loading 'NC_044992.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044992.1.fa' ... DONE Loading 'NC_044993.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044993.1.fa' ... DONE Loading 'NC_044994.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044994.1.fa' ... DONE Loading 'NC_044995.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044995.1.fa' ... DONE Loading 'NC_044996.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044996.1.fa' ... DONE Loading 'NC_044997.1' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_044997.1.fa' ... DONE Loading 'NC_020006.2' sequence from FASTA file '/scratch/nguyenth/bcm/trial/Panubis1.0_og/NC_020006.2.fa' ... DONE Writing all sequences to './BSgenome.Panubis.NCBI.Panubis1.0/inst/extdata/single_sequences.2bit' ... Error in .TwoBits_export(mapply(.DNAString_to_twoBit, object, seqnames), : UCSC library operation failed In addition: Warning messages: 1: In call_fun_in_BSgenomeForge("forgeBSgenomeDataPkg", ...) : forgeBSgenomeDataPkg() has moved to the BSgenomeForge package. Please call BSgenomeForge::forgeBSgenomeDataPkg() to get rid of this warning. 2: In .TwoBits_export(mapply(.DNAString_to_twoBit, object, seqnames), : mustOpen: Can't open ./BSgenome.Panubis.NCBI.Panubis1.0/inst/extdata/single_sequences.2bit to write: No such file or directory
1 answer
I have solved it! It was a corrupt installation of BSgenome. Here is the GitHub that helped me: https://github.com/Bioconductor/BSgenomeForge/pull/46 enter link description here
Log in to answer this question.