# 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_1266.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1266 [L=156] Description: # 1383134 # 1383601 # -1 # ID=1_1266;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.231 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5e-10 39.3 0.0 6.9e-10 38.8 0.0 1.3 1 Acetyltransf_10 Acetyltransferase (GNAT) domain 9e-09 35.5 0.1 1.1e-08 35.2 0.1 1.1 1 Acetyltransf_1 Acetyltransferase (GNAT) family 1.2e-08 35.2 0.1 2.1e-08 34.4 0.1 1.4 1 Acetyltransf_7 Acetyltransferase (GNAT) domain 0.0045 16.5 0.2 0.0073 15.9 0.1 1.4 1 MOZ_SAS MOZ/SAS family ------ inclusion threshold ------ 0.017 15.2 0.5 0.022 14.9 0.5 1.1 1 Acetyltransf_4 Acetyltransferase (GNAT) domain 0.05 13.0 0.5 0.084 12.3 0.1 1.5 2 Dala_Dala_lig_C D-ala D-ala ligase C-terminus 0.057 14.0 0.0 0.064 13.8 0.0 1.2 1 Acetyltransf_3 Acetyltransferase (GNAT) domain Domain annotation for each model (and alignments): >> Acetyltransf_10 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 38.8 0.0 2.9e-13 6.9e-10 20 116 .. 44 145 .. 13 150 .. 0.81 Alignments for each domain: == domain 1 score: 38.8 bits; conditional E-value: 2.9e-13 Acetyltransf_10 20 ealrerieegealllvaeeggkivGvielrdgs......hisllfVdpeyqgqGigkaLleaveeeaekdglkleltvnaspy 96 +++ r++++ ++++ +++k vG+i++ d + +is +f+ pey+ +Gi+++ ++ +ee +++ + t+ + + 23726_FusoPortal_Gene_1266 44 SNIEMRLKQNFTFFYFICVDNKKVGAIRIIDYKekyknkRISPIFILPEYRNKGIAQSAIKICEEIHGNNNWE-LSTILQEKG 125 356666766777777777888********7777799999*******************999999997666555.66777777* PP Acetyltransf_10 97 avpfYeklGFkavgeeqeen 116 +YeklG++ +g++q+ n 23726_FusoPortal_Gene_1266 126 NCYLYEKLGYHSTGKTQVIN 145 **************999765 PP >> Acetyltransf_1 Acetyltransferase (GNAT) family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 35.2 0.1 4.6e-12 1.1e-08 8 117 .] 28 135 .. 15 135 .. 0.76 Alignments for each domain: == domain 1 score: 35.2 bits; conditional E-value: 4.6e-12 HHHHTHHHHHHHHHHHHHHHHHSTTS..EEEEEEETTEEEEEEEEEEEETTSTTEEEEEEEEESGGGTTSSHHHHHHHHHHHH CS Acetyltransf_1 8 lvesffppnedaseeefeesfsleae..gllvaeedgeliGfisllkipdekggvayieslaVspeyrgkGlGsaLlealeer 88 l +++ ++ + ase + + +l + ++++++ d++ +G i++ + + k+++ i +++ peyr+kG+++ ++ +ee+ 23726_FusoPortal_Gene_1266 28 LDKYQDFETNPASEPISNIEMRLKQNftFFYFICVDNKKVGAIRIIDYKE-KYKNKRISPIFILPEYRNKGIAQSAIKICEEI 109 44444444455555553333444444346666***********9999999.9**************************99998 PP HHHTTESEEEEEEETTCHHHHHHHHHTTE CS Acetyltransf_1 89 areegikkivlaveedndaairlYeklGF 117 +++ +l + ++ + lYeklG+ 23726_FusoPortal_Gene_1266 110 HGNNNW---ELSTILQEKGNCYLYEKLGY 135 888844...66666778888********8 PP >> Acetyltransf_7 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 34.4 0.1 8.8e-12 2.1e-08 6 75 .. 58 136 .. 52 137 .. 0.67 Alignments for each domain: == domain 1 score: 34.4 bits; conditional E-value: 8.8e-12 EEEEETTEEEEEEEEEEES....TCTEE....EEE-GGGTTSSHHHHHHHH.TT.SEEEEE..EEEETHHHHHHHHHHTTEE CS Acetyltransf_7 6 fVaeddgkivGfialsled....aeitr....laVdpeyrgkGiGraLlea.ae.kgiile..qlelnaekaikfYeklGFe 75 + + d+k vG i++ ++ + ++++peyr+kGi++ +++ +e +g + + l+++ + +YeklG++ 23726_FusoPortal_Gene_1266 58 YFICVDNKKVGAIRIIDYKekykN---KrispIFILPEYRNKGIAQSAIKIcEEiHGNNNWelSTILQEKGNCYLYEKLGYH 136 4456789999*9999999875431...25678******************642244422221223333334455******96 PP >> MOZ_SAS MOZ/SAS family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.9 0.1 3.1e-06 0.0073 48 105 .. 53 109 .. 45 127 .. 0.70 Alignments for each domain: == domain 1 score: 15.9 bits; conditional E-value: 3.1e-06 MOZ_SAS 48 pFlfYvlteadeeg...ahivGyFskekesaeeynlaCilvlPpyqrkgyGklliefsYeL 105 +F+f+ ++ +d+++ ++i+ y kek +++ ++ i +lP+y++kg + i++ e+ 23726_FusoPortal_Gene_1266 53 NFTFFYFICVDNKKvgaIRIIDY--KEK--YKNKRISPIFILPEYRNKGIAQSAIKICEEI 109 45554444455544223668888..555..55899***************99999885554 PP >> Acetyltransf_4 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.9 0.5 9.1e-06 0.022 3 137 .. 7 140 .. 5 146 .. 0.78 Alignments for each domain: == domain 1 score: 14.9 bits; conditional E-value: 9.1e-06 Acetyltransf_4 3 aatqeDleailrly..aeyrvnt.aitfeyevssveefekrlesl.dsekqylflvfesgrlvGyatlrqfkkakdwkaelsv 81 +t++D ++i+++ a+++ ++fe +++e + +e++ +++ +++ + + + vG+ ++ k +++k + + 23726_FusoPortal_Gene_1266 7 KITLKDTKEIWEMQvkAFKELLDkYQDFET-NPASEPIS-NIEMRlKQNFTFFYFICVDNKKVGAIRIIDYKEKYKNKRISPI 87 578999999999886555444441444554.55555554.566665778899999999*******999988888999999999 PP Acetyltransf_4 82 ylardnnkeglgrellnalvevlkemnienllaeiasnninaivflkalGFekvgi 137 ++ + +++g+++ +++ ++e+ ++n e +++ +++++ ++++lG+ +g+ 23726_FusoPortal_Gene_1266 88 FILPEYRNKGIAQSAIKICEEIHGNNNWE---LSTILQEKGNCYLYEKLGYHSTGK 140 99999****************98777764...566668999999999999998887 PP >> Dala_Dala_lig_C D-ala D-ala ligase C-terminus # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.3 0.0 1 2.5e+03 131 149 .. 9 27 .. 6 30 .. 0.76 2 ? 12.3 0.1 3.5e-05 0.084 108 171 .. 70 131 .. 44 133 .. 0.85 Alignments for each domain: == domain 1 score: -2.3 bits; conditional E-value: 1 -HHHHHHHHHHHHHHHHHT CS Dala_Dala_lig_C 131 ekeveekiqelalkaykal 149 + + +++i e+ +ka+k l 23726_FusoPortal_Gene_1266 9 TLKDTKEIWEMQVKAFKEL 27 5567889999999999876 PP == domain 2 score: 12.3 bits; conditional E-value: 3.5e-05 SSSS-HHHHHC.TT-EEESS----HHHHHHHHHHHHHHHHHTT--SEEEEEEEEETTTEEEEEE CS Dala_Dala_lig_C 108 kgfydyeakysddkaqlevpaklekeveekiqelalkaykaldlrGlaRvDfflteegeivlnE 171 ++dy++ky+ + ++ l++ +++ i + a+k+ + ++++ + +l+e+g+ +l E 23726_FusoPortal_Gene_1266 70 IRIIDYKEKYK--NKRISPIFILPEYRNKGIAQSAIKICEEIHGNNNWELSTILQEKGNCYLYE 131 4799******7..4455555678999999********************************988 PP >> Acetyltransf_3 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.8 0.0 2.7e-05 0.064 4 103 .. 5 105 .. 2 134 .. 0.66 Alignments for each domain: == domain 1 score: 13.8 bits; conditional E-value: 2.7e-05 Acetyltransf_3 4 LrplteeDlealfelrsdpevmry...lpe..pltleeaeewlerawaeeseeegyg..waivlkdtgviGsvglrdleeede 79 L ++t +D +++ e ++ +++ + ++ + + +e + +e ++ +++++ + i++ ++ +G +++ d++e+ + 23726_FusoPortal_Gene_1266 5 LMKITLKDTKEIWE-MQVKAFKELldkYQDfeTNPASEPISNIE--MRL---KQNFTffYFICVDNK-KVGAIRIIDYKEKYK 80 55566666666655.44444444322111122333333344444..333...244444477777666.6************99 PP Acetyltransf_3 80 eaelg.ymlapsfwgkGyateavra 103 ++ + + p++++kG+a+ a + 23726_FusoPortal_Gene_1266 81 NKRISpIFILPEYRNKGIAQSAIKI 105 9999967777**********99875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (156 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 742 (0.0443992); expected 334.2 (0.02) Passed bias filter: 452 (0.0270464); expected 334.2 (0.02) Passed Vit filter: 41 (0.00245333); expected 16.7 (0.001) Passed Fwd filter: 7 (0.000418861); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 7 [number of targets reported over threshold] # CPU time: 0.24u 0.23s 00:00:00.47 Elapsed: 00:00:00.31 # Mc/sec: 1462.89 // [ok]