# 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_894.txt # target HMM database: merged_hmms/CBM13.hmm # per-dom hits tabular output: hmmscan_out/domain_result/CBM13/CBM13_cluster_894.out.dm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 199132|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 [L=266] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.2e-17 48.9 0.0 2.1e-10 27.3 0.0 2.2 2 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 ! 19.2 0.0 6.5e-08 6.5e-08 94 133 .. 86 130 .. 78 134 .. 0.84 2 ! 27.3 0.0 2.1e-10 2.1e-10 9 92 .. 135 235 .. 130 238 .. 0.72 Alignments for each domain: == domain 1 score: 19.2 bits; conditional E-value: 6.5e-08 CBM13.hmm 94 gyykirnrasgkcLdvkgg..adgtpvqqwdc...ngganQqWsv 133 g + i+++ +gk+Ldv gg t++++w++ n+++nQ W 199132|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 86 GIFYIKSALNGKVLDVRGGsrKPDTRIILWPQkfgNDRDNQLWFA 130 567799999********9875447*******999999******54 PP == domain 2 score: 27.3 bits; conditional E-value: 2.1e-10 CBM13.hmm 9 lvnvnsgkcl.........ggstaaganvqlwdcc.n..ganQqWtltsdg 47 +vnv+s +cl + + a + +++++ + + ++nQ++ l++dg 199132|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 135 IVNVGSRLCLdirgslvpdDLRQLAETRIIQYTRKsRveAHNQRFGLSPDG 185 799999999998765555322224445688888774474566999999999 PP CBM13.hmm 48 tirlenhqsgkcLdvaasstaagaevqqytc....ngg..anQqWrlesvg 92 +i h + ++ dv + s+ +ga++++y++ ngg + Q+W e+v+ 199132|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 186 YIYCLAHPN-LVFDVRGESDRDGATLIIYPRkhsdNGGavKHQRWLCEPVS 235 999888886.89999988878999999999944433335558999888774 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (266 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: 16.10 // Query: 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 [L=256] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.2e-16 45.9 0.0 3.6e-10 26.5 0.0 2.1 2 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 ! 24.5 0.0 1.5e-09 1.5e-09 49 123 .. 78 158 .. 72 162 .. 0.75 2 ! 26.5 0.0 3.6e-10 3.6e-10 9 92 .. 125 225 .. 121 246 .. 0.74 Alignments for each domain: == domain 1 score: 24.5 bits; conditional E-value: 1.5e-09 CBM13.hmm 49 irlenhqsgkcLdvaasstaagaevqqytcn.gg..anQqW.rl 88 + +++ +gk+Ldv ++s + +++++++ g+ +nQ W + 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 78 FYIKSGLNGKVLDVRGGSRKPDTRIILWPQKfGTdrDNQLWfAV 121 557788899****9977768*********85333459****444 PP CBM13.hmm 89 esvgsgyykirnrasgkcLdvkgg..adg......tpvqqwdc 123 + i+n+ s +cLd+ g+ d t+++q+++ 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 122 --D----SFIVNVGSRLCLDIRGSlvPDDlqllaeTRIIQYTR 158 ..3....369************989944323555678888876 PP == domain 2 score: 26.5 bits; conditional E-value: 3.6e-10 CBM13.hmm 9 lvnvnsgkcl........ggst.aaganvqlwdcc.n..ganQq 40 +vnv+s +cl + + a + +++++ + + ++nQ+ 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 125 IVNVGSRLCLdirgslvpDDLQlLAETRIIQYTRKsRveAHNQR 168 78999999998875554422223455568888866447456699 PP CBM13.hmm 41 WtltsdgtirlenhqsgkcLdvaasstaagaevqqytc....ng 80 + l++dg+i h + ++ dv + s+ +ga++++y++ ng 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 169 FGLSPDGYIYCLAHPN-LVFDVRGESDRDGATLIIYPRklaeNG 211 9999999888888875.899998888789999999999666656 PP CBM13.hmm 81 g..anQqWrlesvg 92 g + Q+W e+v+ 325274|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 212 GavKHQRWLCEPVS 225 66778999888875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (256 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: 16.46 // Query: 305687|Gonbut1_GeneCatalog_proteins_20160804.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.3e-18 52.8 0.6 5.9e-14 38.9 0.2 2.3 2 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 0.2 5.9e-14 5.9e-14 45 133 .. 4 95 .. 2 101 .. 0.84 2 ! 18.7 0.0 9.3e-08 9.3e-08 51 132 .. 55 140 .. 49 152 .. 0.70 Alignments for each domain: == domain 1 score: 38.9 bits; conditional E-value: 5.9e-14 CBM13.hmm 45 sdgtirlenhqsgkcLdvaasstaagaevqqytcngg...anQqWrlesv 91 +dg++ + ++ +g+++dv++ s+ +ga+v +++ ++ +nQ W ++++ 305687|Gonbut1_GeneCatalog_proteins_20160804.aa.fasta|CBM13 4 PDGYFYIISKANGRAVDVDGASQKDGATVLLWPQKHNtdrDNQLW-FSED 52 6898889999999******77779*********76555779****.7667 PP CBM13.hmm 92 gsgyykirnrasgkcLdvkgg..adgtpvqqwdc...ngganQqWsv 133 + ++n++s+k+Ld+ gg + ++++q + + + nQ+W 305687|Gonbut1_GeneCatalog_proteins_20160804.aa.fasta|CBM13 53 ----GFLVNKSSQKVLDIRGGnvEESAKIIQFRRkviTEALNQRWGY 95 ....78***999*******87777789******98877788****55 PP == domain 2 score: 18.7 bits; conditional E-value: 9.3e-08 CBM13.hmm 51 lenhqsgkcLdva.asstaagaevqqytc...ngganQqWrlesvgsgyy 96 l+n++s+k+Ld+ ++ a+++q + + + nQ+W ++ + + 305687|Gonbut1_GeneCatalog_proteins_20160804.aa.fasta|CBM13 55 LVNKSSQKVLDIRgGNV-EESAKIIQFRRkviTEALNQRWGYS-D----G 98 56777789999875554.6888888888844446689999993.3....3 PP CBM13.hmm 97 kirnras.gkcLdvkggadgtpvqqwdc....ng.ganQqWs 132 ++ a+ + +L v +++dgtp+++ ++ QqW 305687|Gonbut1_GeneCatalog_proteins_20160804.aa.fasta|CBM13 99 VVYCLADpKFVLSVRSNDDGTPLILTRKkfqfDQqHEHQQWL 140 455557767999999889999999988755553335556774 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.01u 0.00s 00:00:00.01 Elapsed: 00:00:00.00 # Mc/sec: 11.66 // Query: 221220|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 [L=253] 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.2e-17 50.9 0.0 2.8e-17 49.7 0.0 1.5 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 ! 49.7 0.0 2.8e-17 2.8e-17 5 139 .. 93 236 .. 89 243 .. 0.82 Alignments for each domain: == domain 1 score: 49.7 bits; conditional E-value: 2.8e-17 CBM13.hmm 5 gtyrlvnvnsgkcl...ggstaaganvqlwdcc.n..gan.QqWtltsdgt 48 g + +++v s+k+l ggs++a+++++lw+ + + +++ Q W ++ +g 221220|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 93 GLFYIRSVLSHKVLdvrGGSKKADTQIVLWSQKkDgiQTQnQLWLYH-EG- 141 668899********999888899999*****8863475555***999.44. PP CBM13.hmm 49 irlenhqsgkcLdvaassta...agaevqqytcngg...anQqWrlesvgs 93 ++n++s ++Ld+ +s + + a+++qy++ + nQ + + ++ 221220|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 142 -FVVNVNSRLVLDIRGSLSPndlDTAKIIQYNRKEKeeaHNQLFAF--DE- 188 .578898889***8755544578*********97666667******..53. PP CBM13.hmm 94 gyykirnras.gkcLdvkgg.ad.gtpvqqwdcng.ganQqWsvlksdlk 139 ++i +a ++Ldv+ + +d g v+++++++ + +Q Ws+++++++ 221220|Umbra1_GeneCatalog_proteins_20121109.aa.fasta|CBM13 189 --GRIYCQAKpHLVLDVADKaSDeGVEVIINERDDdRVSQIWSFEEVQAS 236 ..6899988867******8886679*******9999********998875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (253 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: 20.77 // Query: 361462|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 [L=253] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.4e-17 48.8 0.0 3.7e-16 46.1 0.0 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 ! 46.1 0.0 3.7e-16 3.7e-16 5 139 .. 93 236 .. 90 243 .. 0.81 Alignments for each domain: == domain 1 score: 46.1 bits; conditional E-value: 3.7e-16 CBM13.hmm 5 gtyrlvnvnsgkcl...ggstaaganvqlwdcc.n...ganQqW 41 g + +++v s+k+l ggs++a+++++lw+ + + + nQ W 361462|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 93 GLFYIRSVLSHKVLdvrGGSKKADTQIVLWSQKkDgmqTLNQLW 136 5688999*******999888899999*****8863463355999 PP CBM13.hmm 42 tltsdgtirlenhqsgkcLdvaassta...agaevqqytcngg. 81 ++ +g ++n++s ++Ld+ +s + + a+++qy++ + 361462|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 137 LYH-EG--FIVNVNSRLVLDIRGSLSPndlDTARIIQYNRKEKe 177 999.45..577999889***8755544578*********97555 PP CBM13.hmm 82 ..anQqWrlesvgsgyykirnras.gkcLdvkgg..adgtpvqq 120 nQ + + ++ ++i +a ++Ldv+ + + g +++ 361462|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 178 daHNQLFAY--DE---GRIYCQAKpHLVLDVADKssEEGVEIII 216 447******..53...6899988867******7785559***** PP CBM13.hmm 121 wdcng.ganQqWsvlksdlk 139 +++++ + +Q Ws+++++++ 361462|Umbsp_AD052_1_GeneCatalog_proteins_20160912.aa.fasta|CBM13 217 NERDDdRVSQIWSFEEVQAS 236 **9999********998875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (253 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: 18.62 // [ok]