Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Problems encountered in using DNDSCV:[E::faidx_fetch_seq2] Failed to retrieve block. (Seeking in a compressed, .gzi unindexed, file?) Error in h(simpleError(msg, call)) : 在为'strsplit'函数选择方法时评估'x'参数出了错: record 1 (9:12285795-12285926) failed #88

Open
manchengyanyu opened this issue Mar 20, 2023 · 0 comments

Comments

@manchengyanyu
Copy link

path_cds_table = system.file("extdata", "5.txt", package = "dndscv", mustWork = TRUE)
path_genome_fasta = system.file("extdata", "ref.fa", package = "dndscv", mustWork = TRUE)
buildref(cdsfile=path_cds_table, genomefile=path_genome_fasta, outfile = "example_output_refcds.rda", excludechrs="MT")
[1/3] Preparing the environment...
[2/3] Building the RefCDS object...
[E::faidx_fetch_seq2] Failed to retrieve block. (Seeking in a compressed, .gzi unindexed, file?)
Error in h(simpleError(msg, call)) :
在为'strsplit'函数选择方法时评估'x'参数出了错: record 1 (9:12285795-12285926) failed
file: E:/R/R-4.2.1/library/dndscv/extdata/ref.fa
In addition: Warning message:
In buildref(cdsfile = path_cds_table, genomefile = path_genome_fasta, :
35670 unique gene IDs (column 1) found. 16245 unique gene names (column 2) found. Consider combining gene names and gene IDs or replacing gene names by gene IDs to avoid losing genes (see useids argument in ? buildref)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant