# hmmscan :: search sequence(s) against a profile database # HMMER 3.1b2 (February 2015); http://hmmer.org/ # Copyright (C) 2015 Howard Hughes Medical Institute. # Freely distributed under the GNU General Public License (GPLv3). # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query sequence file: AA/23726_FusoPortal_Gene_39.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_39 [L=324] Description: # 50976 # 51947 # -1 # ID=1_39;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.307 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 9.2e-36 121.6 0.9 1.8e-35 120.7 0.9 1.5 1 ETF_alpha Electron transfer flavoprotein FAD-binding dom 1.2e-32 113.2 0.3 1.9e-32 112.5 0.3 1.3 1 ETF Electron transfer flavoprotein domain ------ inclusion threshold ------ 0.02 15.2 4.0 0.29 11.5 0.0 2.5 3 Tex_YqgF Tex protein YqgF-like domain 0.074 13.5 1.5 3.5 8.0 0.3 2.7 3 Usp Universal stress protein family 0.097 12.2 0.0 0.15 11.6 0.0 1.4 1 FAA_hydrolase Fumarylacetoacetate (FAA) hydrolase family Domain annotation for each model (and alignments): >> ETF_alpha Electron transfer flavoprotein FAD-binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 120.7 0.9 5.5e-39 1.8e-35 1 84 [] 198 281 .. 198 281 .. 0.98 Alignments for each domain: == domain 1 score: 120.7 bits; conditional E-value: 5.5e-39 ETF_alpha 1 vdLeeakvvvsgGrGlkseenlelleelaealgaavgasRavvdagwipeerqvGqtGktvaPklyiavGisGaiqhlaGmkes 84 v+Leea+v+v+gGrG++s+en+el+++la+ +g+ vga+R++++++wip+++qvGq+Gk+vaPklyia+GisGa+qh+ G+ s 23726_FusoPortal_Gene_39 198 VNLEEAEVIVAGGRGMGSKENFELVKQLADVCGGVVGATRPATEDEWIPRSHQVGQSGKIVAPKLYIACGISGATQHVSGIMGS 281 68******************************************************************************9876 PP >> ETF Electron transfer flavoprotein domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 112.5 0.3 5.6e-36 1.9e-32 2 181 .. 6 180 .. 5 181 .. 0.93 Alignments for each domain: == domain 1 score: 112.5 bits; conditional E-value: 5.6e-36 EEEE-EE.TTEE-HHHHHHHHHHHHHHHHHT--EEEEEEES....--THHHHGGGT-SEEEEEE-GGGSS--HHHHHHHHHHHHH CS ETF 2 lvlaehgvdgklnpvdlealeaArrlkekaggevtavvvGppeaeealrealaamGadkvilvddpalagydaeatakalaalik 86 +v++e d++ +v+leal+ A+++ ++++++v av++G+ + +ea +++ +Gad+v++v+ + + ea+ +al ++ + 23726_FusoPortal_Gene_39 6 MVYIETT-DNSPINVSLEALTLAKKISKENNEQVMAVLIGE-NLDEAAKKCFD-YGADEVLYVEENKKE---LEAVGNALIDVKE 84 7899997.555556************99888**********.55555477777.************996...99*********** PP HH--SEEEEESSHHHHHHHHHHHHHTT---EEEESEEEEE.EEEEEEEECTTTTEEEEE--S.SSSEEEEE-TTSS----S-TT- CS ETF 87 kedpdlvlagataigkqvaprlAalLgvplitdvteleveggkltvtRpieggnglaevvea.dlpavltvrpgafepaaimaak 170 k++p+++++g+t gk++a +A+++++p d+++++ e+gk +t p+++gn+l+ev+ d++ v+ vr+ga++++ + ak 23726_FusoPortal_Gene_39 85 KYNPSVIFLGSTINGKDLANIVASKMKTPAFVDAVNVKYENGKYLMTLPMYSGNILKEVTFDgDKTLVIAVRSGACKKEVCENAK 169 ****************************************9*************************************99***** PP --EEEE--HHH CS ETF 171 kkeveevevad 181 ++ev++ ++ad 23726_FusoPortal_Gene_39 170 SGEVKKEKAAD 180 ******99987 PP >> Tex_YqgF Tex protein YqgF-like domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.3 0.1 3.3 1.1e+04 70 92 .. 23 45 .. 19 46 .. 0.67 2 ? 11.5 0.0 8.7e-05 0.29 27 82 .. 59 114 .. 47 130 .. 0.86 3 ? 4.7 0.3 0.011 37 22 78 .. 170 226 .. 151 234 .. 0.70 Alignments for each domain: == domain 1 score: -3.3 bits; conditional E-value: 3.3 HHHHHHHHHHHGGG--EEEEE-- CS Tex_YqgF 70 teklvaelikeeelkvkyvivsE 92 + +l++++ ke++++v +v++ E 23726_FusoPortal_Gene_39 23 ALTLAKKISKENNEQVMAVLIGE 45 55666777777777777777766 PP == domain 2 score: 11.5 bits; conditional E-value: 8.7e-05 EEEEE--STTT--HHHHHHHHHHHHHHTT--EEEEESSTTHHHHHHHHHHHHHHGG CS Tex_YqgF 27 etaviyphepkkkveeakkklkklvkkykveliaiGnGtasreteklvaelikeee 82 + +v+y++e+kk++e+ ++l ++ +ky+ ++i +G +++ ++va+ +k+ + 23726_FusoPortal_Gene_39 59 ADEVLYVEENKKELEAVGNALIDVKEKYNPSVIFLGSTINGKDLANIVASKMKTPA 114 5789*********************************************9999853 PP == domain 3 score: 4.7 bits; conditional E-value: 0.011 TS-EEEEEEE--STTT--HHHHHHHHHHHHHHTT--EE.EEESSTTHHHHHHHHHHHH CS Tex_YqgF 22 tgkvletaviyphepkkkveeakkklkklvkkykveli.aiGnGtasreteklvaeli 78 +g+v +++ + +k++e +++++ v+ + e+i a G G +s+e +lv++l 23726_FusoPortal_Gene_39 170 SGEVKKEKAADK-NLFTKIAEIVQEISETVNLEEAEVIvAGGRGMGSKENFELVKQLA 226 455544444433.345678888899999999999887625679**********99875 PP >> Usp Universal stress protein family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 0.9 0.2 0.16 5.4e+02 14 67 .. 19 72 .. 1 85 [. 0.67 2 ? 3.1 0.1 0.035 1.2e+02 76 112 .. 60 96 .. 19 102 .. 0.69 3 ? 8.0 0.3 0.001 3.5 71 116 .. 170 215 .. 141 228 .. 0.84 Alignments for each domain: == domain 1 score: 0.9 bits; conditional E-value: 0.16 HHHHHHHHHHHHHHHCTSEEEEEEEEESSHCCCCCCCHHHHHHHHHHHHHHHHH CS Usp 14 eskkAlewAlelakargaelillhvidprasgaasdadeeaeeeelelaeaele 67 s +Al A + k+++ +++++ + + + +a + d a e ++++ +++ 23726_FusoPortal_Gene_39 19 VSLEALTLAKKISKENNEQVMAVLIGENLDEAAKKCFDYGADEVLYVEENKKEL 72 466788888888888888888888877777777666666666666666665554 PP == domain 2 score: 3.1 bits; conditional E-value: 0.035 HCTTCEEEEEEEESSHHHHHHHHHHHTTSSEEEEESC CS Usp 76 eaggvkvevvveegdpaeeilevaeeedadliVvGsr 112 ++ ++e + e ++ ++++v e+++++ i +Gs+ 23726_FusoPortal_Gene_39 60 DEVLYVEENKKELEAVGNALIDVKEKYNPSVIFLGST 96 4444455556666677888899999999999999986 PP == domain 3 score: 8.0 bits; conditional E-value: 0.001 HHHHHHCTTCEEEEEEEESSHHHHHHHHHHHTTSSEEEEESCCSST CS Usp 71 akaaaeaggvkvevvveegdpaeeilevaeeedadliVvGsrgrsg 116 + ++++ ++ + ++ + +ei e+++ e+a+ iV G rg+++ 23726_FusoPortal_Gene_39 170 SGEVKKEKAADKNLFTKIAEIVQEISETVNLEEAEVIVAGGRGMGS 215 4567889999999999***************************876 PP >> FAA_hydrolase Fumarylacetoacetate (FAA) hydrolase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 4.6e-05 0.15 144 187 .. 170 213 .. 104 245 .. 0.82 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 4.6e-05 TTEEEEEEEG.GGTTTSHHHHHHHHHHHH-B-TT-EEEEEES... CS FAA_hydrolase 144 ngeevqdsnt.kdmifspaeliahlsagltLepGdliltGTpsGv 187 +ge++++ + k++ +++ae+++++s+++ Le ++i++G G 23726_FusoPortal_Gene_39 170 SGEVKKEKAAdKNLFTKIAEIVQEISETVNLEEAEVIVAGGR-GM 213 6788888888799**************************987.43 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (324 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 965 (0.0577429); expected 334.2 (0.02) Passed bias filter: 611 (0.0365606); expected 334.2 (0.02) Passed Vit filter: 43 (0.002573); expected 16.7 (0.001) Passed Fwd filter: 5 (0.000299186); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.29u 0.24s 00:00:00.53 Elapsed: 00:00:00.32 # Mc/sec: 2943.37 // [ok]