Entering edit mode
4.7 years ago
geizetomazetto
•
0
Hi there,
I have several genomes fasta files for screening based on HMM. I wrote this very simple script. It seems that works - the output files are right. However, I also got this message
"Incorrect number of command line arguments. Usage: hmmscan [-options] <hmmdb> <seqfile>
where most common options are: -h : show brief help on version and usage
To see more help on available options, do /Users/alexjones/miniconda3/bin/hmmscan -h"
Someone knows why I get this message? The script is wrong?
I appreciate any tips. Thanks.
for file in ./*.fasta;
do
hmmscan --domtblout ${file}.out.dm Target_HMM.txt ${file} > ${file}.out;
done
