# hmmscan :: search sequence(s) against a profile database # HMMER 3.3.2 (Nov 2020); http://hmmer.org/ # Copyright (C) 2020 Howard Hughes Medical Institute. # Freely distributed under the BSD open source license. # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query sequence file: /array2/siva/dbcansub_update/cazy2025/match_results/dbsub_data_2025/merged_clusters/CBM13/fasta_for_each_cluster/cluster_583.txt # target HMM database: merged_hmms/CBM13.hmm # per-dom hits tabular output: hmmscan_out/domain_result/CBM13/CBM13_cluster_583.out.dm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 176332|Micbo1_GeneCatalog_proteins_20151118.aa.fasta|CBM13 [L=163] Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.7e-15 42.0 10.1 5.8e-14 38.9 10.1 2.0 1 CBM13.hmm Domain annotation for each model (and alignments): >> CBM13.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 38.9 10.1 5.8e-14 5.8e-14 6 134 .. 9 159 .. 4 162 .. 0.79 Alignments for each domain: == domain 1 score: 38.9 bits; conditional E-value: 5.8e-14 CBM13.hmm 6 tyrlvnvnsgkcl..ggst.aaganvqlwdccn.ganQqWtltsdg..... 47 +++l+n + g cl +s+ ++g++v+ +++ n ++nQ W+l+++ 176332|Micbo1_GeneCatalog_proteins_20151118.aa.fasta|CBM13 9 VIALINYGAGTCLdlTASDpKSGNPVIGYNY-NaTSNQHWKLVKVPysask 58 67899********884444488888**9*95.54788******66566778 PP CBM13.hmm 48 ...tirlenhqsgkcLdvaasstaagaevqqytcngg.....anQqWrles 90 t+ l n+ +g +Ld ++ ++g++++ +++ g+ + Q Wr+++ 176332|Micbo1_GeneCatalog_proteins_20151118.aa.fasta|CBM13 59 vwpTWFLINQATGTALDLLNGGEDNGTQIVGWAN-GAasktnPHQLWRIIT 108 99988888888889***6566668*********6.66455445899***98 PP CBM13.hmm 91 vgsg..yykirnrasgkcLdvkgg..adgtpvqqwdc...ngganQqWsvl 134 +++ + i+n+ ++ +d+ +g ++ t+v+ w+ n+++ Q W+v+ 176332|Micbo1_GeneCatalog_proteins_20151118.aa.fasta|CBM13 109 ADNNfdVVAIQNVGTNTFVDLYNGspQNSTKVTGWAGdisNQNPHQLWKVK 159 8644435999**999*******8898889*******766466889999875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (163 residues searched) Target model(s): 1 (188 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.00u 0.00s 00:00:00.00 Elapsed: 00:00:00.00 # Mc/sec: 11.05 // Query: 432238|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 [L=185] Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.7e-19 57.0 6.0 2.3e-18 53.3 5.6 1.9 1 CBM13.hmm Domain annotation for each model (and alignments): >> CBM13.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 53.3 5.6 2.3e-18 2.3e-18 7 145 .. 11 170 .. 9 182 .. 0.80 Alignments for each domain: == domain 1 score: 53.3 bits; conditional E-value: 2.3e-18 CBM13.hmm 7 yrlvnvnsgkcl..ggst.aaganvqlwdccnganQqWtltsdg...... 47 +++n ++g cl +s+ ++g++v+ +++++++nQ W+l++++ 432238|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 11 LAFINYGTGTCLdlTASDpKSGTPVIGYRYNGTSNQHWKLVKVDyagakv 60 578999******884444599999*****666888******88778899* PP CBM13.hmm 48 ..tirlenhqsgkcLdvaasstaagaevqqytcngg...anQqWrlesvg 92 + l n+ +g +Ld + s +g++++ +t +g + Q Wrl +++ 432238|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 61 fpVYWLINQATGTALDLLNASKENGTQINGWTTGDGkpnPTQLWRLVTTN 110 *988888888889***5555558*********866677668*****9998 PP CBM13.hmm 93 .sgy.ykirnrasgkcLdvk.gg.adgtpvqqwdc...ngganQqWsvlk 135 +g+ + i+n+ sg +d+ g+ ad t+++ w+ + ++ Q W+v+ 432238|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 111 aNGAvVMIQNVGSGTFVDLYlGNaADRTKITGWAGniaEENPHQLWRVKV 160 89986888************777999********8454457789998877 PP CBM13.hmm 136 sdlkikliae 145 +++ ++ ++ 432238|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 161 ISQDDTNDTR 170 7666665554 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (185 residues searched) Target model(s): 1 (188 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.00u 0.00s 00:00:00.00 Elapsed: 00:00:00.00 # Mc/sec: 13.16 // Query: 2966|Pesfi1_GeneCatalog_proteins_20210226.aa.fasta|CBM13 [L=167] Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7e-20 58.2 10.2 1.6e-18 53.8 10.0 2.0 1 CBM13.hmm Domain annotation for each model (and alignments): >> CBM13.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 53.8 10.0 1.6e-18 1.6e-18 8 135 .. 12 163 .. 10 166 .. 0.84 Alignments for each domain: == domain 1 score: 53.8 bits; conditional E-value: 1.6e-18 CBM13.hmm 8 rlvnvnsgkcl...ggstaaganvqlwdccnganQqWtltsdg......tirl 51 +++n +s c+ gg+ +g ++ w c++g+nQ+W+l+++g +++ 2966|Pesfi1_GeneCatalog_proteins_20210226.aa.fasta|CBM13 12 TFTNLGSRTCIdlaGGNSGNGVAISGWKCHGGSNQKWQLQQVGktgpwpVYMI 64 678999*****987666688889*****998899******99999****99** PP CBM13.hmm 52 enhqsgkcLdva.asstaagaevqqytcngg.....anQqWrlesvg..sgyy 96 +n+ sg +d +++ a+g+++ ++ +++ + Q+Wr+ + + +g++ 2966|Pesfi1_GeneCatalog_proteins_20210226.aa.fasta|CBM13 65 KNQYSGTYMDLYlGGT-ADGTPIYGWQGSNTagsgnPHQRWRFVTGDpqNGNV 116 **99*9****986765.******999976443334469*****7665789996 PP CBM13.hmm 97 ki.rnrasgkcLdvk.gg.adgtpvqqwdc.....ngganQqWsvlk 135 ++ +n g +d+ gg a+gt + w+ +g++ QqW+v+ 2966|Pesfi1_GeneCatalog_proteins_20210226.aa.fasta|CBM13 117 VLiENIGAGTYVDLYlGGvANGTGINGWSLagdlsSGNPHQQWTVTV 163 66599888*******777***********8899978899****9875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (167 residues searched) Target model(s): 1 (188 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.00u 0.00s 00:00:00.00 Elapsed: 00:00:00.00 # Mc/sec: 12.83 // Query: 398464|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 [L=164] Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-15 43.7 2.4 2e-15 43.7 2.4 1.9 1 CBM13.hmm Domain annotation for each model (and alignments): >> CBM13.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 43.7 2.4 2e-15 2e-15 52 153 .. 13 126 .. 3 129 .. 0.85 Alignments for each domain: == domain 1 score: 43.7 bits; conditional E-value: 2e-15 CBM13.hmm 52 enhqsgkcLdvaasstaagaevqqytcngganQqWrlesvg.sg.....y 95 n sg cLd +as +g +v+ y++n ++nQ+W+l +v+ sg 398464|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 13 INYGSGTCLDLTASDPKSGNPVIGYNYNATSNQHWKLVKVPySGakvwpV 62 5667899***7766656**********999*******9998888899887 PP CBM13.hmm 96 ykirnrasgkcLdvkgg..adgtpvqqwdc....ngganQqWsvlksdlk 139 y ++n+a+g +Ld+ +g ++gt+++ w++ +++ Q W++ ++d++ 398464|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 63 YFLINQATGTALDLYNGgqENGTQITGWANgpasKTNPHQLWRLVTADAN 112 889999*********888999********9888555889*********99 PP CBM13.hmm 140 ikliaekvipdgnn 153 + ++ + i g++ 398464|Mictri1_GeneCatalog_proteins_20170510.aa.fasta|CBM13 113 SDVVMIQNIGTGTY 126 99999998887776 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (164 residues searched) Target model(s): 1 (188 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.00u 0.00s 00:00:00.00 Elapsed: 00:00:00.00 # Mc/sec: 12.49 // Query: 6479|Neodi1_GeneCatalog_proteins_20170711.aa.fasta|CBM13 [L=167] Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.9e-15 42.7 2.6 9.9e-15 41.4 2.5 1.6 1 CBM13.hmm Domain annotation for each model (and alignments): >> CBM13.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 41.4 2.5 9.9e-15 9.9e-15 5 131 .. 10 161 .. 6 165 .. 0.78 Alignments for each domain: == domain 1 score: 41.4 bits; conditional E-value: 9.9e-15 CBM13.hmm 5 gtyrlvnvnsgkcl...ggst...aaganvqlwdccnganQqWtltsdg.... 47 +t+++vn +++ + gg + +g++ + wd + nQqW+l++ 6479|Neodi1_GeneCatalog_proteins_20170711.aa.fasta|CBM13 10 KTVAFVNFGTNTAIdlsGG-HcdpPDGTPCIGWDAVLSENQQWKLVKLRdgpg 61 5678888888888866422.24568888988899554699******6555588 PP CBM13.hmm 48 ...tirlenhq.sgkcLdvaasstaagaevqqytcngg...anQqWrlesvg. 92 +r++n++ gk+Ld ++s++++++++ ++ g+ + Q W + +vg 6479|Neodi1_GeneCatalog_proteins_20170711.aa.fasta|CBM13 62 glpLFRIQNVKaPGKALDLYNGSSSDNTKITGWSH-GSttnDHQLWYIHPVGe 113 9**9*******77*****8777779*********7.555667899******77 PP CBM13.hmm 93 ...sgy.ykirnrasgkcLdvk.gg.adgtpvqqwdc..ng.ganQqW 131 g+ + ++n+ ++ +d+ gg a+gt v+ w+ g+ Q W 6479|Neodi1_GeneCatalog_proteins_20170711.aa.fasta|CBM13 114 fqnRGNvVMLENYGTNTFVDLYlGGaANGTSVVGWQGslKAlGSHQLW 161 7776664778************767999*****999633333456666 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (167 residues searched) Target model(s): 1 (188 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 1 (1); expected 0.0 (0.02) Passed Vit filter: 1 (1); expected 0.0 (0.001) Passed Fwd filter: 1 (1); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.01u 0.00s 00:00:00.01 Elapsed: 00:00:00.00 # Mc/sec: 13.00 // [ok]