-
Notifications
You must be signed in to change notification settings - Fork 32
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
WorkflowError - missing file #73
Comments
What CADD-scripts version you are using? If 1.7 then please use only Snakemake 7.x . For CADD-scripts1.7.1. your Snakemake version should be fine. (From your command I am pretty sure you use 1.7 and not 1.7.1. so please use the latest cadd scripts version. A simple upgrade of the repo should be fine. No other data needed). Do you run Snakemake locally or in a cluster environment? Some environments have difficulties to run it in /tmp. Maybe you don't use the CADD.sh script avoiding the /tmp directory. Then you have to modify this command: snakemake /tmp/tmp.Zmoud0EEpB/input.tsv.gz --use-conda --conda-prefix /data/analysis/src/CADD-scripts-1.7/envs/conda --cores 1 |
@visze Thank you for the comments. I am running CADD 1.7 and following the README.md directions per https://github.com/kircherlab/CADD-scripts, which specify using Snakemake version 8. I am running the CADD.sh script. I tried setting TMPDIR=~/caddtmp to avoid using /tmp, but getting similar missing file error. I tried CADD 1.7.1 with different/worse results. README.md mentions using apptainer/singularity, but no specifics. |
I am very sure you use CADD-scripts v1.7. For CADD-scripts v1.7.1 your command should look like: Lines 148 to 151 in 77df69a
But It looks like v1.7: Lines 126 to 127 in 203ee3b
Snakemake v1.7 requires snakemake 7.X which is mentioned in it's readme: Lines 56 to 59 in 203ee3b
You are referring to the latest (CADD-scripts v1.7.1 release) Readme which is correct: there snakemake 8.X sould be used. Apptainer will only work with CADD-scripts v1.7.1 and it is the default in Can you show me the Please |
Installed CADD 1.7.1 in /data/workspace/bchesnut/CADD-scripts-1.7.1 Linked /data/workspace/bchesnut/CADD-scripts-1.7.1/data to data location:
Set some environment variables:
Ran ./install.sh Ran ./CADD.sh -p -a -g GRCh38 -o ./output_inclAnno_GRCh38.tsv.gz ./test/input.vcf
|
Ok. Two things I see. First can you bgzip your input file to input.vcf.gz. but not sure if it change anything. Second which can be a trouble maker too: Paths have to be correctly set for apptainer images (apptainer command --bind). You have to bind a lot of them extra. E.g. tmp folder,.... Otherwise tmp of singularity image is used and then wiped and next time loaded not there anymore. I tested it on my end..it worked but you never know on other systems... So maybe first recommendation is to use only mamba first (-m) flag in the CADD.sh script? |
I am getting the following error while running the test script for CADD 1.7:
Verbose output using -p is attached. cadd-output.txt
I'm running Red Hat EL9 and Miniforge3 conda with snakemake 8.20.3
Thank you in advance for suggestions.
The text was updated successfully, but these errors were encountered: