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

Offline use of countSplitReads #50

Closed
AEstebanMar opened this issue Feb 22, 2023 · 2 comments
Closed

Offline use of countSplitReads #50

AEstebanMar opened this issue Feb 22, 2023 · 2 comments

Comments

@AEstebanMar
Copy link

AEstebanMar commented Feb 22, 2023

Hello, I am having some problems running your package due to (lack of) internet connectivity.

Our group works in a network which blocks internet connection for security reasons. We can download any necessary resources locally, but when it comes to execution we have no connection available.

The issue comes with the countSplitReads function, which ends up calling get_chrom_info_for_registered_UCSC_genome from genomeInfoDb. This function attempts to download chromInfo.txt.gz, which is not possible in our setting.

An issue was submitted to GenomeInfoDb describing a related problem (Bioconductor/GenomeInfoDb#26), and a solution was implemented to the package (quote reply). Would it be possible for countSplitReads to have a parameter passed down through function calls to use this "offline mode"? If so, we could first download chromosome information locally and access that cache in offline execution.

Thank you,
Álvaro.

@drewjbeh
Copy link
Collaborator

drewjbeh commented Jan 22, 2024

Hi @AEstebanMar,

I am trying to identify the source of the error but am having a hard time. countSplitReads doesn't actually call get_chrom_info_for_registered_UCSC_genome as you mention. I am looking for other function calls that try to access this info online. As part of that search, I have run FRASER with the test data on my local machine and disabled internet connectivity when doing so. This completes successfully.
A few things that may be the issue; which version of FRASER are you running? Which species/genome does the data originate from?
It would also be very useful to send some other things here to help debugging. For example; the list of commands you use up to the point you get an error, and the traceback and error output at the problematic step. Thanks!

@AEstebanMar
Copy link
Author

Hello, @drewjbeh,

I have updated FRASER since I originally posted the issue, and did not think to undo my temporary fix. I've tested it again after reading your answer, this time removing my fix, and it works properly. It seems the issue has been fixed at some point since my original question.

Thank you, marking as solved.

This issue was closed.
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

2 participants