# 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/CBM18/fasta_for_each_cluster/cluster_110.txt # target HMM database: merged_hmms/CBM18.hmm # per-dom hits tabular output: hmmscan_out/domain_result/CBM18/CBM18_cluster_110.out.dm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 582565|Pirfi3_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=544] 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.4e-06 13.2 17.5 8e-06 12.3 17.5 1.5 1 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 12.3 17.5 8e-06 8e-06 2 35 .. 499 535 .. 498 541 .. 0.89 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 8e-06 CBM18.hmm 2 cgkqaggntCplnaCCskwGyC...GttsdYCg..kgCq 35 cg +ag+++C +n+CCs GyC Gt+s C+ +gCq 582565|Pirfi3_GeneCatalog_proteins_20160330.aa.fasta|CBM18 499 CGIHAGESKCSDNKCCSILGYCqkyGTSS--CKidNGCQ 535 999*******************9998888..96548999 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (544 residues searched) Target model(s): 1 (38 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: 39.44 // Query: 326471|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=568] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 0.00017 8.1 17.7 0.00027 7.4 17.7 1.4 1 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 7.4 17.7 0.00027 0.00027 2 36 .. 523 559 .. 522 564 .. 0.91 Alignments for each domain: == domain 1 score: 7.4 bits; conditional E-value: 0.00027 CBM18.hmm 2 cgkqaggntCplnaCCskwGyCGttsdYCg..kgCqs 36 cg +ag n+C n+CC GyC + ++ C+ +gCq 326471|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 523 CGINAGYNRCGRNLCCTFNGYCVNDENSCKieNGCQF 559 9*9**************************76589*96 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (568 residues searched) Target model(s): 1 (38 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: 62.07 // Query: 15095|PirE2_1_GeneCatalog_proteins_20110421.aa.fasta|CBM18 [L=616] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- [No hits detected that satisfy reporting thresholds] Domain annotation for each model (and alignments): [No targets detected that satisfy reporting thresholds] Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (616 residues searched) Target model(s): 1 (38 nodes) Passed MSV filter: 1 (1); expected 0.0 (0.02) Passed bias filter: 0 (0); expected 0.0 (0.02) Passed Vit filter: 0 (0); expected 0.0 (0.001) Passed Fwd filter: 0 (0); expected 0.0 (1e-05) Initial search space (Z): 1 [actual number of targets] Domain search space (domZ): 0 [number of targets reported over threshold] # CPU time: 0.00u 0.00s 00:00:00.00 Elapsed: 00:00:00.00 # Mc/sec: 147.19 // Query: 267109|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=589] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8e-06 12.3 20.4 1.3e-05 11.6 20.4 1.4 1 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.6 20.4 1.3e-05 1.3e-05 2 37 .. 545 582 .. 544 586 .. 0.93 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 1.3e-05 CBM18.hmm 2 cgkqaggntCplnaCCskwGyCGttsdYCg..kgCqsn 37 cg +ag n+C +++CCs+ GyC ++ C+ +gCqsn 267109|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 545 CGINAGYNKCGDKSCCSSDGYCVKDENSCKieNGCQSN 582 9*9**************************76589**98 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (589 residues searched) Target model(s): 1 (38 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: 73.44 // Query: 267108|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=591] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 0.00019 7.9 17.7 0.00032 7.2 17.7 1.4 1 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 7.2 17.7 0.00032 0.00032 2 36 .. 546 582 .. 545 587 .. 0.88 Alignments for each domain: == domain 1 score: 7.2 bits; conditional E-value: 0.00032 CBM18.hmm 2 cgkqaggntCplnaCCskwGyCGttsdYCg..kgCqs 36 cg +ag n+C n+CC + GyC + C+ +gCq 267108|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 546 CGINAGYNRCGRNLCCASDGYCVKDEYSCKvqNGCQF 582 9*9*******************998888964479995 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (591 residues searched) Target model(s): 1 (38 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: 72.24 // Query: 292121|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=578] 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.4e-08 21.1 15.8 2.8e-08 20.2 15.8 1.5 1 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.2 15.8 2.8e-08 2.8e-08 2 35 .. 532 567 .. 531 570 .. 0.94 Alignments for each domain: == domain 1 score: 20.2 bits; conditional E-value: 2.8e-08 CBM18.hmm 2 cgkqaggntCplnaCCskwGyCGttsdYCgk..gCq 35 cg +ag n+C++n+CCs GyC + +d C++ gCq 292121|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 532 CGINAGFNKCQNNLCCSAEGYCSNDADSCQNdlGCQ 567 9*9**************************9877999 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (578 residues searched) Target model(s): 1 (38 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: 63.19 // Query: 411511|Pirfi3_GeneCatalog_proteins_20160330.aa.fasta|CBM18|CBM18 [L=724] 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-12 33.7 51.9 4.1e-11 29.2 22.9 2.7 2 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 9.6 21.0 5.9e-05 5.9e-05 2 36 .. 591 628 .. 590 633 .. 0.88 2 ! 29.2 22.9 4.1e-11 4.1e-11 1 37 [. 642 676 .. 642 677 .. 0.94 Alignments for each domain: == domain 1 score: 9.6 bits; conditional E-value: 5.9e-05 CBM18.hmm 2 cgkqaggntCplnaCCskwGyC.GttsdYCgk..gCqs 36 cg ++g + C +++CCs GyC ++s YC++ gCq 411511|Pirfi3_GeneCatalog_proteins_20160330.aa.fasta|CBM18|CBM18 591 CGINSGLKSCGSSKCCSDDGYCyKSSSSYCKTsqGCQI 628 99999*****************888899**85449*96 PP == domain 2 score: 29.2 bits; conditional E-value: 4.1e-11 CBM18.hmm 1 qcgkqaggntCplnaCCskwGyCGttsdYCgkgCqsn 37 +cgk g C++++CCs +GyCGt+s+YC++gCqsn 411511|Pirfi3_GeneCatalog_proteins_20160330.aa.fasta|CBM18|CBM18 642 RCGK--GYGSCKSGYCCSYYGYCGTSSQYCDAGCQSN 676 6997..889**************************97 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (724 residues searched) Target model(s): 1 (38 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: 35.84 // Query: 243538|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 [L=591] 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-06 15.0 18.6 1.2e-06 15.0 18.6 2.0 2 CBM18.hmm Domain annotation for each model (and alignments): >> CBM18.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.6 0.1 0.37 0.37 25 31 .. 273 279 .. 272 283 .. 0.77 2 ! 15.0 18.6 1.2e-06 1.2e-06 2 36 .. 545 581 .. 544 588 .. 0.93 Alignments for each domain: == domain 1 score: -2.6 bits; conditional E-value: 0.37 CBM18.hmm 25 ttsdYCg 31 +t++YC 243538|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 273 YTEKYCV 279 6899996 PP == domain 2 score: 15.0 bits; conditional E-value: 1.2e-06 CBM18.hmm 2 cgkqaggntCplnaCCskwGyCGttsdYCg..kgCqs 36 cg +ag+n+C +++CCs GyC ++ C+ +gCqs 243538|Anasp1_GeneCatalog_proteins_20160330.aa.fasta|CBM18 545 CGINAGNNKCGNDLCCSYNGYCVDDKNLCKieNGCQS 581 9****************************76589**9 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (591 residues searched) Target model(s): 1 (38 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: 75.05 // [ok]