# 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_1394.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1394 [L=181] Description: # 1502099 # 1502641 # 1 # ID=1_1394;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.247 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-37 128.7 0.1 2.6e-36 124.9 0.1 1.9 1 Formyl_trans_N Formyl transferase ------ inclusion threshold ------ 0.053 13.9 0.0 12 6.3 0.0 2.1 2 CoA_binding_2 CoA binding domain 0.063 13.4 0.0 0.11 12.7 0.0 1.3 1 DUF2483 Hypothetical protein of unknown function (DUF Domain annotation for each model (and alignments): >> Formyl_trans_N Formyl transferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 124.9 0.1 4.7e-40 2.6e-36 2 181 .] 3 175 .. 2 175 .. 0.84 Alignments for each domain: == domain 1 score: 124.9 bits; conditional E-value: 4.7e-40 EEEEEESESSHHHHHHHHHHHHTTSECEEEEEEESCTTEHCHHHHHHCCHCCCECHHCCSCCCCCCCHHHHHHHHCTT-SEEE CS Formyl_trans_N 2 kivvlisgtgsnlqalldalkkgsleadvvlvvsnkdkaaglekaeaasiptkvfeekkleprsqfeeelaeelkelkadliv 84 ki+vl+sg+g+n+ +l+ k++++ + +++ +++ +a + a + +i + l+++++ ++l + +e k+dliv 23726_FusoPortal_Gene_1394 3 KIIVLVSGSGTNMLQLI----KNNIK--IDCIIADRECKA-KNIADEYKI-----DFVLLNRNKEISKNLLKIFEERKPDLIV 73 67788888887765544....45554..445555544333.222222222.....33336777899999************** PP EES--S---HHHHHSSTTSEEEEESSSTTTT.....BSS-HHHHHHHTTSSEEEEEEEE--SSSSTS-EEEEEEEE--TT--H CS Formyl_trans_N 85 lagymkllpkelvqavrlkilniHpSLLPrf.....kGaaaiqraleaGdketGvtihfvdeelDtGailaqkkveiiaedtk 162 lag++ +l +e+++++++ki+niHpSLLP++ +G++++q+++e+Gdke+G+t+h+v+ ++D G+i+ q+kv+i ++ 23726_FusoPortal_Gene_1394 74 LAGFLSILDGEILEKYKNKIINIHPSLLPKYggkgmYGLKVHQAVFENGDKESGCTVHYVTSNVDAGEIIGQEKVDISMAKSP 156 *********************************99999********************************************* PP HHHHHHHHHHHHHHHHHHH CS Formyl_trans_N 163 eeleqrvadaehkalveal 181 ee+++ v + e k+l++++ 23726_FusoPortal_Gene_1394 157 EEIQKIVLEREWKLLPRVV 175 ****************996 PP >> CoA_binding_2 CoA binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 6.3 0.0 0.0022 12 47 103 .. 37 94 .. 24 100 .. 0.78 2 ? 5.8 0.0 0.0031 18 25 76 .. 125 179 .. 107 180 .. 0.79 Alignments for each domain: == domain 1 score: 6.3 bits; conditional E-value: 0.0022 ---T....T--EEEE-S-HHHHGGGHHHHHHH--SEEEE..-TT---HHHHHHHHHTT- CS CoA_binding_2 47 aslaedpepvdlvdvftppeatpeivkealalgakavwl..qpgsedeeaaelaeeagi 103 +++a d+ ++d v + ++ e+ +++k ++++++ ++l ++++ d e++e+ +++ i 23726_FusoPortal_Gene_1394 37 KNIA-DEYKIDFVLLNRNKEISKNLLKIFEERKPDLIVLagFLSILDGEILEKYKNKII 94 4556.5699***********************999888722345777777777776655 PP == domain 2 score: 5.8 bits; conditional E-value: 0.0031 HHT--EEEESSS-..SEETT.EES----T....T--EEEE-S-HHHHGGGHHHHH CS CoA_binding_2 25 ehglpvvpVnpke..eeilg.ekvyaslaedpepvdlvdvftppeatpeivkeal 76 e+g++v+ V+ + +ei g ekv s+a++pe+++ +++ + + +p +vk+++ 23726_FusoPortal_Gene_1394 125 ESGCTVHYVTSNVdaGEIIGqEKVDISMAKSPEEIQKIVLEREWKLLPRVVKKLI 179 38899999955555599999566888999779*******************9986 PP >> DUF2483 Hypothetical protein of unknown function (DUF2483) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 2e-05 0.11 22 64 .. 122 165 .. 112 168 .. 0.85 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 2e-05 DUF2483 22 rPtevndtvkystdlrdarefdGledavvdmskh.kaikktvte 64 e tv y t+ da e G e+ + m+k + i+k v e 23726_FusoPortal_Gene_1394 122 GDKESGCTVHYVTSNVDAGEIIGQEKVDISMAKSpEEIQKIVLE 165 5668889************************9863678888876 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (181 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 542 (0.0324318); expected 334.2 (0.02) Passed bias filter: 363 (0.0217209); expected 334.2 (0.02) Passed Vit filter: 28 (0.00167544); 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.24u 0.23s 00:00:00.47 Elapsed: 00:00:00.32 # Mc/sec: 1644.29 // [ok]