# 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_852.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_852 [L=335] Description: # 948809 # 949813 # -1 # ID=1_852;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.297 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-71 240.5 0.3 2.3e-71 240.1 0.3 1.2 1 NMT1 NMT1/THI5 like 9.6e-05 22.1 0.3 0.0011 18.6 0.0 2.4 2 NMT1_2 NMT1-like family ------ inclusion threshold ------ 0.018 14.5 0.2 0.028 13.8 0.2 1.2 1 VitK2_biosynth Menaquinone biosynthesis Domain annotation for each model (and alignments): >> NMT1 NMT1/THI5 like # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 240.1 0.3 4.1e-75 2.3e-71 1 216 [] 47 261 .. 47 261 .. 0.99 Alignments for each domain: == domain 1 score: 240.1 bits; conditional E-value: 4.1e-75 -GGGHHHHHHHHTTHHHHTTEEEEEE--SSTTHHHHHHHCTS-SEEEE-HHHHHHHHHTT--EEEEEEEE-S--EEEEEECCC- CS NMT1 1 npnhaglyvAkekGYFkeeGldVeivepadpsdatelvasgkadfgvsyqesllaarakglpvvsvaallqeplsglivlkdsg 84 n+nh+glyvAkekGYF eeG+d++i++pa+ +++ +lv ++ka +gv+ q+ + + akg+p++++aa++++++sg+i++k+ + 23726_FusoPortal_Gene_852 47 NTNHTGLYVAKEKGYFSEEGIDLDIKQPAN-ESTSDLVINNKAPMGVYFQDYMSSKLAKGAPITAIAAIIENNTSGIITNKKLN 129 89**************************99.99*************************************************** PP -SSGGGGTTSEEEE-S-TTHHHHHHHHHHHTT--GGG-EEEE-.TT-HHHHHHCTSSSEEECEETTTHHHHHHHTT--EEEE-G CS NMT1 85 ikspkDLkGkrigyvgspfekallkallakdggdpsdvtlvnvggtslakalltgkvDaaiggyenwegvelesegkelrifal 168 i+spk+L+G+++g++ p+e +l+ +++kd gd+s+v+lv+++ ++ +++l +g +Daa +y++w++++ +s g+e+++f++ 23726_FusoPortal_Gene_852 130 INSPKELAGHKYGTWDIPIELGMLQFIIEKDSGDFSKVELVPNTDDNSITPLSNGAFDAAP-VYYAWDKIMGDSLGIETNFFYY 212 ************************************************************9.********************** PP GGTT--.---SEEEEEEHHHHHH-HHHHHHHHHHHHHHHHHHHHSHHHH CS NMT1 169 adlgvp.dyyslvlitnekllkenpelvkkFlrAlkkGyeyalahPeea 216 +d+ + ++ys+v+i+n+++lk+n+e + k+lrA+kkGy+ya++hPeea 23726_FusoPortal_Gene_852 213 KDYAPElNFYSPVIIANNDYLKDNKEEATKILRAIKKGYQYAMEHPEEA 261 ************************************************8 PP >> NMT1_2 NMT1-like family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 18.6 0.0 1.9e-07 0.0011 10 83 .. 42 110 .. 34 146 .. 0.75 2 ? 0.5 0.1 0.067 3.7e+02 142 142 .. 269 269 .. 225 325 .. 0.54 Alignments for each domain: == domain 1 score: 18.6 bits; conditional E-value: 1.9e-07 NMT1_2 10 lGfialtdaaplvvakekglfakeGltvelekqaswaelrdkllagelDaahv.ltpapllidlGiggakadmkv 83 l ++++t+++ l vakekg+f++eG+++++ kq + +++d +++++ +++ + + + + +G ++ + 23726_FusoPortal_Gene_852 42 LDWVPNTNHTGLYVAKEKGYFSEEGIDLDI-KQPANESTSDLVINNKAPMGVYfQDYMSSKLAKGA-----PITA 110 5689***********************885.67777899******999999886666667777776.....3333 PP == domain 2 score: 0.5 bits; conditional E-value: 0.067 NMT1_2 142 w 142 + 23726_FusoPortal_Gene_852 269 A 269 2 PP >> VitK2_biosynth Menaquinone biosynthesis # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.8 0.2 5e-06 0.028 188 241 .. 228 285 .. 218 294 .. 0.77 Alignments for each domain: == domain 1 score: 13.8 bits; conditional E-value: 5e-06 EECCCHH...HHHHHHHHHHHHHHHHHHHSHHHHHHHHHCCSTTS.H.HHHHHHHHHH CS VitK2_biosynth 188 vrrdlaa...eaakeleealreSveyalahldevaeyaaeraqemdd.evlekhlelY 241 +++d+ + e+a+++ +a+++ ++ya++h++e+ae + ++a e ++ + + ++Y 23726_FusoPortal_Gene_852 228 ANNDYLKdnkEEATKILRAIKKGYQYAMEHPEEAAEILIKYAPELENkKAMIIESQKY 285 4444333345889*****************************9997755555555555 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (335 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 776 (0.0464337); expected 334.2 (0.02) Passed bias filter: 581 (0.0347654); expected 334.2 (0.02) Passed Vit filter: 51 (0.0030517); 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.28u 0.23s 00:00:00.51 Elapsed: 00:00:00.32 # Mc/sec: 3043.30 // [ok]