Open avportal opened 1 month ago
My condition is similar to yours.It shows: [1] "hmmalign --trim --dna -o RasperGade16S_align//seq.align --mapali /home/tom/R/x86_64-pc-linux-gnu-library/3.6/RasperGade16S/extdata/aligned_remove_dot.fasta /home/tom/R/x86_64-pc-linux-gnu-library/3.6/RasperGade16S/extdata/16S_core_gg.hmm /home/tom/R/x86_64-pc-linux-gnu-library/3.6/RasperGade16S/extdata/Demo/demo.SILVA.fasta" sh: 1: hmmalign: not found Error in system(command = cmd, intern = TRUE) : error in running command And I don't know why it doesn't work...
Hi there,
I am trying to run pred.GCN = predict_16SGCN_from_sequences(), following the information provided @ https://github.com/wu-lab-uva/RasperGade16S
However, I get the following:
No FASTA file supplied, running prediction using demo sequences in /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library/RasperGade16S/extdata/Demo/demo.SILVA.fasta
[1] "hmmalign --trim --dna -o RasperGade16S_align//seq.align --mapali /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library/RasperGade16S/extdata/aligned_remove_dot.fasta /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library/RasperGade16S/extdata/16S_core_gg.hmm /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library/RasperGade16S/extdata/Demo/demo.SILVA.fasta" [1] "esl-reformat -o RasperGade16S_align//seq.afa -u --gapsym=- afa RasperGade16S_align//seq.align" sh: esl-reformat: command not found Error in system(command = cmd, intern = TRUE) : error in running command
Any assistance will be highly appreciated.
Angel