# 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_599.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_599 [L=342] Description: # 664474 # 665499 # 1 # ID=1_599;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.292 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-98 329.4 29.4 2e-98 329.1 29.4 1.0 1 FecCD FecCD transport family ------ inclusion threshold ------ 0.15 11.8 1.0 0.81 9.5 1.0 2.0 1 BAAT_C BAAT / Acyl-CoA thioester hydrolase C terminal 1.7 8.4 3.1 11 5.8 0.3 2.4 2 DUF3767 Protein of unknown function (DUF3767) Domain annotation for each model (and alignments): >> FecCD FecCD transport family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 329.1 29.4 3.7e-102 2e-98 1 312 [] 14 333 .. 14 333 .. 0.98 Alignments for each domain: == domain 1 score: 329.1 bits; conditional E-value: 3.7e-102 FecCD 1 llllllllslllgavkispsevfeallgka........dsleaeivlelrlprvlaallvGaaLavaGallQtvtrNpLadPsi 76 ll+ ++++sl++gav +++s+v++++l+k+ +++ iv+ +r+pr+++a++vG+aLa++G+ +Q++++Np++d++i 23726_FusoPortal_Gene_599 14 LLIGIIIFSLFYGAVRVPISDVIKIILNKTglfnfeisKQSYIPIVFFVRFPRIMVAVIVGGALALCGCTMQSLLKNPIVDSGI 97 578899**********************9999999988888888**************************************** PP FecCD 77 lGissgaalgavlalllllgslslfalplaaflgallalllvlllarkrkgasplkliLaGiavsallsaltsllllllneeal 160 +Gissga+lgav+a+ l+l+++++fa+p+++ +al+++++++ +++ r+++++l liL+Gia++++++a+ts++l+ l e+++ 23726_FusoPortal_Gene_599 98 IGISSGASLGAVIAVSLGLTATNIFAMPIFSGTFALIISAIIYKISTLRGKTDNLLLILSGIAIGSFVGAITSVILTSLAETEM 181 ************************************************************************************ PP FecCD 161 esilfwllGslsgaswskvllllvllllllllllllarkldllalgdetAkslGvnvervrllllllvalltavavalvGpigF 244 ++++fw +Gsl+g++w+++++ l+ +++l ++l+++ ++l++l lg+e+AkslG+n++++r +l+++allta++v+++G+i+F 23726_FusoPortal_Gene_599 182 KEYIFWAMGSLNGRRWEHFFFGLIPITILSPILFYYGKELNILLLGEEEAKSLGINIKKIRGKILIIIALLTAISVCISGNITF 265 ************************************************************************************ PP FecCD 245 vGLivPhiarllvgadhrkllplsallGallllladllarlllsptelpvgivtaliGapffiylllr 312 vGLivPhi+r+++g+d+rkll +s+l+Ga++l++ dll+r++l+p e++vgivtaliGap+fiyl+++ 23726_FusoPortal_Gene_599 266 VGLIVPHILRKIIGSDNRKLLKSSFLAGACFLTFGDLLSRVVLAPKEISVGIVTALIGAPYFIYLIVK 333 ******************************************************************97 PP >> BAAT_C BAAT / Acyl-CoA thioester hydrolase C terminal # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 9.5 1.0 0.00015 0.81 12 68 .. 86 143 .. 83 251 .. 0.88 Alignments for each domain: == domain 1 score: 9.5 bits; conditional E-value: 0.00015 HHHCSTTB--SSEEEE.EETHHHHHHHHHHHHSSSEEEEEEES--SB--SSEEEETTE CS BAAT_C 12 yllrhpkvkkdkigll.GiSkGgelalllasflkkisavvlingsavvvegklvykde 68 ll++p v ig+ G S G+ +a++l+ + ++i a+ + +g+++ + +yk + 23726_FusoPortal_Gene_599 86 SLLKNPIVDSGIIGISsGASLGAVIAVSLGLTATNIFAMPIFSGTFALIISAIIYKIS 143 588999999999999548999***********************99999988888855 PP >> DUF3767 Protein of unknown function (DUF3767) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 5.8 0.3 0.002 11 21 46 .. 89 113 .. 82 139 .. 0.81 2 ? 2.6 0.2 0.019 1.1e+02 30 60 .. 156 191 .. 154 198 .. 0.76 Alignments for each domain: == domain 1 score: 5.8 bits; conditional E-value: 0.002 DUF3767 21 siPCaRdslltGigagavvGglrfvl 46 + P + ds+++Gi++ga++G++ +v 23726_FusoPortal_Gene_599 89 KNPIV-DSGIIGISSGASLGAVIAVS 113 55543.899************99885 PP == domain 2 score: 2.6 bits; conditional E-value: 0.019 DUF3767 30 ltGigagavvGglrfvlkk.....nvakaaNWavga 60 l Gi+ g vG+++ v+ + ++++ + Wa+g+ 23726_FusoPortal_Gene_599 156 LSGIAIGSFVGAITSVILTslaetEMKEYIFWAMGS 191 89***********99998855555666666666665 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (342 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 722 (0.0432025); expected 334.2 (0.02) Passed bias filter: 260 (0.0155577); expected 334.2 (0.02) Passed Vit filter: 17 (0.00101723); 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): 3 [number of targets reported over threshold] # CPU time: 0.26u 0.23s 00:00:00.49 Elapsed: 00:00:00.31 # Mc/sec: 3207.11 // [ok]