Build profile HMMs and scan sequences against databases
Build profile HMMs and scan against sequence databases
hmmbuildhmmsearchhmmscanhmmpresshmmfetchhmmstathmmemithmmconverthmmalign$ hmmbuild output.hmm input_alignment.sto$ hmmsearch --tblout results.txt myprofile.hmm protein_db.fasta$ hmmscan --domtblout domain_results.txt Pfam-A.hmm query_sequence.fasta