# 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_680.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_680 [L=113] Description: # 772639 # 772977 # 1 # ID=1_680;partial=00;start_type=ATG;rbs_motif=AGxAGG/AGGxGG;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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.2e-28 99.6 0.0 1.3e-28 99.4 0.0 1.0 1 HIT HIT domain 3.8e-24 85.3 0.0 4.4e-24 85.0 0.0 1.0 1 DcpS_C Scavenger mRNA decapping enzyme C-term binding ------ inclusion threshold ------ 0.015 14.9 0.0 0.021 14.5 0.0 1.2 1 GalP_UDP_tr_C Galactose-1-phosphate uridyl transferase, C-te Domain annotation for each model (and alignments): >> HIT HIT domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 99.4 0.0 2.4e-32 1.3e-28 2 97 .. 11 107 .. 10 108 .. 0.94 Alignments for each domain: == domain 1 score: 99.4 bits; conditional E-value: 2.4e-32 GCC-CCEEEE-SSEEEEE-SS-SSSSEEEEEESSTGCSGGGCC...HHHHHHHHHHHHHHHHHHHTCTTCCCEEEEEEESGGGT CS HIT 2 keipstvvyeddlvlafldidpkapghlLvipkkhlselhdlt...eeelakllsavkkvakalqkkfknkglrvvinngahag 82 keip+++vyedd+v+af di p ap+h+Lv+pkk++++++d+t ++k+ +++ k+ak+++ ++g+rvv n+++h+g 23726_FusoPortal_Gene_680 11 KEIPANIVYEDDDVIAFKDIAPVAPVHVLVVPKKEIPTINDITdedTLLIGKVYRVIGKLAKEFGI--DKDGYRVVSNCNEHGG 92 89*****************************************66556789**************9..789************* PP -SSSS-EEEEEECXT CS HIT 83 qsvahlHiHviprre 97 q+v+h+H+H+i++++ 23726_FusoPortal_Gene_680 93 QTVFHIHFHLIGGEK 107 ************986 PP >> DcpS_C Scavenger mRNA decapping enzyme C-term binding # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 85.0 0.0 8e-28 4.4e-24 3 101 .. 4 104 .. 2 110 .. 0.96 Alignments for each domain: == domain 1 score: 85.0 bits; conditional E-value: 8e-28 HHHHHHTTSS-B-EEEE-SSEEEEE-SS-SSSEEEEEEESS--SSGGG--TTCCCHHHHHHHHHHHHHHHTT-T..-EEEE--- CS DcpS_C 3 vfnilakgkepervvyedekfvvlkDikpkaklhllvipkrdikslrdLtkehlpllkemrekakkvieekyne..vdaeelrv 84 +f +++++++p+++vyed++++++kDi p a++h+lv+pk++i++++d+t+e+ l+ ++++++ k+++e + +++ ++ 23726_FusoPortal_Gene_680 4 LFTKIINKEIPANIVYEDDDVIAFKDIAPVAPVHVLVVPKKEIPTINDITDEDTLLIGKVYRVIGKLAKEFGIDkdGYRVVSNC 87 799*******************************************************************99889999****** PP STTTT--S-S--EEEEE CS DcpS_C 85 gvhyqPsvyhlHlHvis 101 h ++v+h+H+H+i 23726_FusoPortal_Gene_680 88 NEHGGQTVFHIHFHLIG 104 **************997 PP >> GalP_UDP_tr_C Galactose-1-phosphate uridyl transferase, C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.5 0.0 3.7e-06 0.021 33 119 .. 18 101 .. 4 110 .. 0.76 Alignments for each domain: == domain 1 score: 14.5 bits; conditional E-value: 3.7e-06 EEE-SSEEEE--TT--STT-EEEEESS--SSGGG--HHHHHHHHHHHHHHHHHHHHHCTS---EEEEEE---SCCCTTTT--EE CS GalP_UDP_tr_C 33 velkdsflavvpylakWPfevlllrkaevkslaeladeekedlaailkklarkydnlfetsfeysmGlhqaplnaeednswlHa 116 v ++d+ +a + P++vl+++k+e++ +++++de++ + + + k f + + + + n + +++ +H+ 23726_FusoPortal_Gene_680 18 VYEDDDVIAFKDIAPVAPVHVLVVPKKEIPTINDITDEDTLLI-GKVYRVIGKLAKEFGIDKDGYRVV--SNCNEHGGQTVFHI 98 44455455555566677********************988654.567778888889999887644445..56677777888999 PP EEE CS GalP_UDP_tr_C 117 Hfe 119 Hf+ 23726_FusoPortal_Gene_680 99 HFH 101 985 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (113 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 349 (0.0208832); expected 334.2 (0.02) Passed bias filter: 322 (0.0192676); expected 334.2 (0.02) Passed Vit filter: 24 (0.00143609); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.21u 0.23s 00:00:00.44 Elapsed: 00:00:00.31 # Mc/sec: 1059.66 // [ok]