# 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/AA6/fasta_for_each_cluster/cluster_16.txt # target HMM database: merged_hmms/AA6.hmm # per-dom hits tabular output: hmmscan_out/domain_result/AA6/AA6_cluster_16.out.dm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 343409|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 [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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.5e-74 234.0 0.0 6.9e-74 233.7 0.0 1.1 1 AA6.hmm Domain annotation for each model (and alignments): >> AA6.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 233.7 0.0 6.9e-74 6.9e-74 1 193 [. 8 200 .. 8 202 .. 0.97 Alignments for each domain: == domain 1 score: 233.7 bits; conditional E-value: 6.9e-74 AA6.hmm 1 kiaviyYSsyghikqlAeavaeGveea.gaevdlrqvpetlpeevlakmea.a 51 k+++++YS+y h+++lA av+eG+e++ ++ev ++q++et+pee+l++ +a + 343409|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 8 KLYIVFYSMYLHVYKLALAVKEGAERNgNVEVSIFQIKETVPEEILTNIKApS 60 789************************999*********************88 PP AA6.hmm 52 kedvpvatpedLaeadaiifGsPtrfGnvaaqlkaFlDatgglwakgalkgKv 104 k+dvpv+tp+ L+ead+++fG+PtrfG+v+aqlk F+D+ g+lwa+gal+gK+ 343409|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 61 KPDVPVITPALLKEADGVLFGFPTRFGTVPAQLKDFIDSCGQLWATGALNGKL 113 8**************************************************** PP AA6.hmm 105 agvFtstgtlhGGqetTllsllttlahhGmiivglgykeafkelanleevngg 157 ag F st+t+h GqetT++++lt+lah+Gmi+v+ gy+++ +++ +e+ gg 343409|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 114 AGTFNSTQTQHSGQETTAFTFLTVLAHFGMIYVPNGYSSP--LIHEDKEIMGG 164 **************************************87..6799******* PP AA6.hmm 158 spyGastlagadgsrqpseeelelaraqGkrvaeia 193 sp+Gas++a + g ++ps++ele+ + qG+++a+ 343409|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 165 SPWGASVVAVNTGVQSPSDKELEIGKVQGENFAKYL 200 ********************************9975 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (256 residues searched) Target model(s): 1 (195 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: 46.01 // Query: 325478|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 [L=252] 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.6e-78 248.1 0.0 3.3e-78 247.8 0.0 1.1 1 AA6.hmm Domain annotation for each model (and alignments): >> AA6.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 247.8 0.0 3.3e-78 3.3e-78 1 195 [] 8 201 .. 8 201 .. 0.97 Alignments for each domain: == domain 1 score: 247.8 bits; conditional E-value: 3.3e-78 AA6.hmm 1 kiaviyYSsyghikqlAeavaeGveea.gaevdlrqvpetlpeevlakmea.a 51 k+ +++YS+y h+++lA+a++eG+e+ ++ev ++q+petl+eevl+km+a a 325478|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 8 KVFIVFYSTYSHVYKLAQAIKEGAEKEgKVEVSIYQIPETLSEEVLSKMGApA 60 899************************999*********************77 PP AA6.hmm 52 kedvpvatpedLaeadaiifGsPtrfGnvaaqlkaFlDatgglwakgalkgKv 104 k+dvpv+tp+ L+ead+i+fG+PtrfG+v+aqlk F+D+ g+lwa+gal+gK+ 325478|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 61 KPDVPVITPALLKEADGILFGFPTRFGTVPAQLKDFIDSCGQLWATGALNGKL 113 7**************************************************** PP AA6.hmm 105 agvFtstgtlhGGqetTllsllttlahhGmiivglgykeafkelanleevngg 157 +g+Ftst+t+hGGqe+T++++lt+l+h+Gm++v+ gy+++ +++ +ev gg 325478|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 114 GGIFTSTQTQHGGQESTAFTFLTVLTHFGMLYVPNGYSTP--LIQEDKEVIGG 164 **************************************87..57899****** PP AA6.hmm 158 spyGastlagadgsrqpseeelelaraqGkrvaeiakk 195 sp+Gas++a++d ++p+e+e+e+ r+qG+++a+ ++k 325478|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 165 SPWGASVVAKND-VPSPNEKEIEIGRIQGQNFANYVAK 201 **********98.6889****************98876 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (252 residues searched) Target model(s): 1 (195 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: 51.78 // Query: 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 [L=252] 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.9e-76 242.0 0.1 3.3e-76 241.3 0.0 1.4 2 AA6.hmm Domain annotation for each model (and alignments): >> AA6.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 241.3 0.0 3.3e-76 3.3e-76 3 194 .. 10 200 .. 8 201 .. 0.96 2 ? -2.9 0.1 0.22 0.22 33 50 .. 211 228 .. 204 239 .. 0.59 Alignments for each domain: == domain 1 score: 241.3 bits; conditional E-value: 3.3e-76 AA6.hmm 3 aviyYSsyghikqlAeavaeGveea.gaevdlrqvpetlpeevlakmea.ake 53 ++yYS+y h+++lA+a++eG+e+ +++v l+q+petl+eevl km+a ak+ 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 10 FIVYYSTYSHVYKLAKAIQEGAEKEgKVKVSLYQIPETLSEEVLGKMGApAKP 62 79***********************999*********************777* PP AA6.hmm 54 dvpvatpedLaeadaiifGsPtrfGnvaaqlkaFlDatgglwakgalkgKvag 106 dvpv+tp+ L+ead+i+fG+PtrfG+v+aq+k F+D+ g+lw++gal+gK++g 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 63 DVPVITPSLLKEADGILFGIPTRFGTVPAQVKDFIDSCGQLWMTGALNGKLGG 115 ***************************************************** PP AA6.hmm 107 vFtstgtlhGGqetTllsllttlahhGmiivglgykeafkelanleevnggsp 159 +Ftst+t+hGGqe+T++++lt+l+h+Gm++v+ gy+++ +++ +ev ggsp 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 116 IFTSTQTQHGGQESTAFTFLTILTHFGMLYVPNGYSTP--LIQEDKEVIGGSP 166 ************************************87..57899******** PP AA6.hmm 160 yGastlagadgsrqpseeelelaraqGkrvaeiak 194 +Gas++a +d ++p+++e+e+ r+qG+++a+ + 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 167 WGASVVARND-IPSPNDKEIEIGRIQGQNFANYLA 200 ********99.6789***************99766 PP == domain 2 score: -2.9 bits; conditional E-value: 0.22 AA6.hmm 33 lrqvpetlpeevlakmea 50 +v+e +++ev +++ + 361826|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 211 ATKVQEVTKAEVTKNNVQ 228 345666555666655544 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (252 residues searched) Target model(s): 1 (195 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: 55.74 // Query: 297203|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 [L=237] 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-76 241.0 0.0 4.7e-76 240.7 0.0 1.1 1 AA6.hmm Domain annotation for each model (and alignments): >> AA6.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 240.7 0.0 4.7e-76 4.7e-76 1 194 [. 42 235 .. 42 236 .. 0.98 Alignments for each domain: == domain 1 score: 240.7 bits; conditional E-value: 4.7e-76 AA6.hmm 1 kiaviyYSsyghikqlAeavaeGveea.gaevdlrqvpetlpeevlakmea.a 51 k+ +++Y++y hi+++A ++++G+e++ +++v l+q+petl++e+l+k+ a + 297203|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 42 KLFILFYTTYLHIYKMALEIKKGAEKQgKCDVYLYQIPETLSKEALKKVRApP 94 688************************999*********************99 PP AA6.hmm 52 kedvpvatpedLaeadaiifGsPtrfGnvaaqlkaFlDatgglwakgalkgKv 104 k++vp++t++ L++ad+i+fG+Ptr+G ++aqlk F+D+ g+lwakg+l+gK+ 297203|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 95 KPNVPYITTDLLRQADGILFGFPTRYGILPAQLKDFIDSCGQLWAKGELAGKM 147 9**************************************************** PP AA6.hmm 105 agvFtstgtlhGGqetTllsllttlahhGmiivglgykeafkelanleevngg 157 ag F+st+t+hGGqet +l++lttlah+G+i+v+ gy+++ +l+n +e+ gg 297203|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 148 AGTFVSTATQHGGQETSILTFLTTLAHFGIIYVPFGYSSP--HLNNDTEILGG 198 **************************************76..799******** PP AA6.hmm 158 spyGastlagadgsrqpseeelelaraqGkrvaeiak 194 sp+G+s++ag+dg+rqpse+ele+a+ qG++++++ + 297203|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 199 SPWGSSVIAGKDGTRQPSELELEIANLQGENFSKTLA 235 ********************************99875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (237 residues searched) Target model(s): 1 (195 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: 52.43 // Query: 325479|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 [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.1e-66 210.1 0.0 1.4e-66 209.8 0.0 1.1 1 AA6.hmm Domain annotation for each model (and alignments): >> AA6.hmm # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 209.8 0.0 1.4e-66 1.4e-66 1 194 [. 8 200 .. 8 201 .. 0.97 Alignments for each domain: == domain 1 score: 209.8 bits; conditional E-value: 1.4e-66 AA6.hmm 1 kiaviyYSsyghikqlAeavaeGveea.gaevdlrqvpetlpeevlakmea.a 51 k+++++YS+ ++++lA a++eG+e++ +++v+++qv+etl++++lak +a + 325479|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 8 KLLIVFYSTFMNVYKLALAIKEGAEKVnEIDVEIYQVKETLSDSILAKKGApE 60 789************************9***********************77 PP AA6.hmm 52 kedvpvatpedLaeadaiifGsPtrfGnvaaqlkaFlDatgglwakgalkgKv 104 k+d+p++tpe L+ead+i+fG+P fG ++aqlk F+D+ g+lwa+ga++gK+ 325479|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 61 KPDIPIITPELLKEADGILFGFPACFGLLPAQLKDFIDSCGSLWASGAMQGKL 113 7**************************************************** PP AA6.hmm 105 agvFtstgtlhGGqetTllsllttlahhGmiivglgykeafkelanleevngg 157 ag+F st +h GqetT+++++t+lah+Gmi+v+ gy+ + e+++ +e+ gg 325479|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 114 AGIFNSTELQHSGQETTAFTFITILAHFGMIFVPNGYS--CLEIHEDKEILGG 164 **************************************..56899******** PP AA6.hmm 158 spyGastlagadgsrqpseeelelaraqGkrvaeiak 194 sp+Gas++a++ + r+ps++e+++ar qG+++a+ + 325479|Conth1_GeneCatalog_proteins_20141003.aa.fasta|AA6 165 SPWGASVVASK-KIRSPSNKEKAIARCQGENFATYLA 200 *********86.5899****************98765 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (253 residues searched) Target model(s): 1 (195 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: 51.79 // [ok]