Entering edit mode
In my hands, the readGenBank function in the genbankr package works very well with smaller GenBank files (<500 KB), but gives me a C stack usage error with any larger files.
In the example below chrI is 442Kb and chrIII is 638Kb. Am I doing something wrong and/or do the larger files need some pre-treatment?
I've tried other files as well - larger files always throw this error. Small files work just fine. I am currently running genbankr version 1.20.0 on R version 4.1.2.
Thanks!
> library(genbankr)
> chrI <- readGenBank("/Users/USER/Desktop/chrI.gb")
Translation product seems to be missing for 1 of 101 CDS annotations. Setting to ''
No exons read from genbank file. Assuming sections of CDS are full exons
No transcript features (mRNA) found, using spans of CDSs
> chrII <- readGenBank("/Users/USER/Desktop/chrII.gb")
Error: C stack usage 7970208 is too close to the limit
> sessionInfo()
R version 4.1.2 (2021-11-01)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Big Sur 11.5.2
Matrix products: default
LAPACK: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats4 parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] BSgenome_1.60.0 rtracklayer_1.52.1 GenomicAlignments_1.28.0
[4] Rsamtools_2.8.0 Biostrings_2.60.2 XVector_0.32.0
[7] SummarizedExperiment_1.22.0 MatrixGenerics_1.4.3 matrixStats_0.61.0
[10] GenomeInfoDbData_1.2.6 GenomicFeatures_1.44.2 AnnotationDbi_1.54.1
[13] Biobase_2.52.0 GenomicRanges_1.44.0 GenomeInfoDb_1.28.4
[16] IRanges_2.26.0 S4Vectors_0.30.2 BiocGenerics_0.38.0
[19] genbankr_1.20.0