funannotate
funannotate copied to clipboard
Funannotate test error
Are you using the latest release? 1.8.9
Describe the bug When running funannotate test -t all --cpus 90, I endup having the error below, despite of the "funannotate check --show-versions" show everything installed and with all path defined.
[May 16 01:02 PM]: Pulling out high quality Augustus predictions
[May 16 01:02 PM]: Found 314 high quality predictions from Augustus (>90% exon evidence)
[May 16 01:02 PM]: Running SNAP gene prediction, using training data: annotate/predict_misc/busco.final.gff3
[May 16 01:02 PM]: 2 predictions from SNAP
[May 16 01:02 PM]: Running GlimmerHMM gene prediction, using training data: annotate/predict_misc/busco.final.gff3
[May 16 01:04 PM]: 1,773 predictions from GlimmerHMM
[May 16 01:04 PM]: Summary of gene models passed to EVM (weights):
[May 16 01:04 PM]: EVM: partitioning input to ~ 35 genes per partition using min 1500 bp interval
Traceback (most recent call last):
File "/home/dar21/miniconda3/envs/funannotate/lib/python3.8/site-packages/funannotate/aux_scripts/funannotate-runEVM.py", line 479, in funannotate predict
test complete.
#########################################################
#########################################################
Running funannotate predict
BUSCO-mediated training unit testing
CMD: funannotate predict -i test.softmasked.fa --protein_evidence protein.evidence.fasta -o annotate --cpus 90 --species Awesome busco
#########################################################
#########################################################
Traceback (most recent call last):
File "/mnt/mountdoom/mnt/work/bioinfo/home/dar21/miniconda3/envs/funannotate/bin/funannotate", line 10, in
What command did you issue? funannotate test -t all --cpus 90
Logfiles Done it
OS/Install Information Included in the logfile
Sorry for late response, so in your logfile I see that SNAP seems to be failing, I think on debian systems the conda version of snap is broken. You can remove it and install with apt-get
which should fix SNAP and it will hopefully predict genes. I'm not sure, but I think perhaps EVM choked on some incomplete gene model that the broken SNAP executable generated.
[May 16 01:02 PM]: 2 predictions from SNAP
In a functioning version of SNAP you should see something like 1500 gene models predicted.