# 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_900.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_900 [L=265] Description: # 1002898 # 1003692 # -1 # ID=1_900;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.244 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.9e-56 190.4 3.2 3.3e-56 190.2 3.2 1.0 1 RrnaAD Ribosomal RNA adenine dimethylase 2.3e-05 24.9 0.0 3.7e-05 24.3 0.0 1.4 1 Methyltransf_25 Methyltransferase domain 0.00022 21.7 0.0 0.00042 20.8 0.0 1.4 1 Methyltransf_11 Methyltransferase domain 0.00025 20.6 0.0 0.00037 20.1 0.0 1.3 1 NodS Nodulation protein S (NodS) 0.0019 18.0 0.0 0.0028 17.4 0.0 1.5 1 Methyltransf_23 Methyltransferase domain 0.0032 17.1 1.0 0.01 15.5 1.0 1.8 1 Microcin Colicin E1 (microcin) immunity protein ------ inclusion threshold ------ 0.016 14.9 0.2 0.04 13.6 0.0 1.8 2 PCMT Protein-L-isoaspartate(D-aspartate) O-methyl 0.055 12.5 0.1 0.092 11.7 0.1 1.3 1 DREV DREV methyltransferase 0.057 14.0 0.0 0.15 12.7 0.0 1.8 1 Methyltransf_12 Methyltransferase domain 0.3 11.2 0.1 0.3 11.2 0.1 2.9 3 UBA_2 Ubiquitin associated domain (UBA) 0.52 10.1 4.4 0.17 11.6 0.8 2.1 2 UBA_3 Fungal ubiquitin-associated domain Domain annotation for each model (and alignments): >> RrnaAD Ribosomal RNA adenine dimethylase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 190.2 3.2 2.2e-59 3.3e-56 2 263 .. 3 261 .. 2 263 .. 0.93 Alignments for each domain: == domain 1 score: 190.2 bits; conditional E-value: 2.2e-59 GCCXXGCG--EE--H.HHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S- CS RrnaAD 2 nklrksygqnflknk.kvineIldklnleekltvleigpGkGilTeeLakrakqvvvlElderLakllqkklskdeklkvvhqD 84 +k +k+ygqnfl+nk ++n+I++ +n+++++++leigpG+G+lT L++r+k+v+ +E+d++L + l+kk+s++e+ ++v D 23726_FusoPortal_Gene_900 3 FKHKKKYGQNFLNNKdEILNKIIEVSNISDNDEILEIGPGQGALTSLLVERVKKVTCVEIDKDLENTLRKKFSSKENYTLVMGD 86 6899**********84699***************************************************************** PP CCGHHHHHHHCTTTCCCCEEEEEEE--TTTHHHHHHHHHHGGGGCEEEEEEEEEHHHHHHHH--TTSTT-SHHHHHHHHHEEEE CS RrnaAD 85 flkfetpeevtekdsihqeflvvanlPyeistqlvkqlleesrfglvkmllwlqkefarrllarpgskkrsrlsvlreaftdvk 168 +l+ ++++ i+q +vvan+Py+i+++++ +++e + + ++++++qke+ +r++a+ g k+r+ l++++e++ +++ 23726_FusoPortal_Gene_900 87 VLEVDFKK------YINQGTKVVANIPYYITSPIINKIIENKDL-IDEAYIMVQKEVGERICAKSG-KERGILTLAVEYYGESE 162 ********......579**********************85555.9*****************866.***************** PP EEEEC-GGCECS--SSSEEEEEEEE.-SS-S.--S-HHHHHHHHHHHHCTTTSBCCCHCCCTS.TCHHHHHCTT..T--TTTCC CS RrnaAD 169 lvakveksifsPppkvdsalveler.kdedl.pqkdlkkyesvvrkllnrkrktlstslksll.pggelqalsst.gindnalv 248 + +++++ f P p+vdsa + ++ kd + ++ + + + + v++++ krk++ ++l +l+ ++++ + + ++ ++ +n ++ 23726_FusoPortal_Gene_900 163 YLFTIPREFFNPIPNVDSAFISIKFyKDDRYkNKISEDLFFKYVKAAFSNKRKNIVNNLATLGySKDKIKEILNQiEVPENERA 246 ***********************972555557778999**********************99636777788878889******* PP CCSGHHHHHHHHHHH CS RrnaAD 249 kklsaeqtldifkel 263 +++s+++++++++++ 23726_FusoPortal_Gene_900 247 ENISIDKFIELINIF 261 **********99987 PP >> Methyltransf_25 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 24.3 0.0 2.4e-08 3.7e-05 1 62 [. 36 94 .. 36 121 .. 0.88 Alignments for each domain: == domain 1 score: 24.3 bits; conditional E-value: 2.4e-08 EEEET-TTSHHHHHHHCTSTSSEEEEEESSHHHHHHHHHHHHHTTTTEEEEESBCCGBSSST CS Methyltransf_25 1 vLdlGcGtGrltlalakrlgkakvtgvDiseemleaareraeeeglnvefvqgdaedlpfed 62 +L++G G+G lt l +r+ +kvt+v+i++++ ++ r+++++++ n ++v gd+ +++f++ 23726_FusoPortal_Gene_900 36 ILEIGPGQGALTSLLVERV--KKVTCVEIDKDLENTLRKKFSSKE-NYTLVMGDVLEVDFKK 94 8*****************9..9*******************8765.*********9999864 PP >> Methyltransf_11 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.8 0.0 2.8e-07 0.00042 1 57 [. 37 93 .. 37 98 .. 0.95 Alignments for each domain: == domain 1 score: 20.8 bits; conditional E-value: 2.8e-07 EEET-TTSHHHHHHHHTTSEEEEEES-HHHHHHHHHHCCHHT-EEETTEECCCTSS- CS Methyltransf_11 1 LdvGcGtGrlaealakrgarvvgvDlspemlklakekaseeglvefvvadaeklpfp 57 L++G G+G+l+ l +r+++v++v++++ + ++ ++k s +++ ++v +d+ +++f+ 23726_FusoPortal_Gene_900 37 LEIGPGQGALTSLLVERVKKVTCVEIDKDLENTLRKKFSSKENYTLVMGDVLEVDFK 93 89***************************************99889****9999887 PP >> NodS Nodulation protein S (NodS) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.1 0.0 2.4e-07 0.00037 43 121 .. 32 112 .. 12 155 .. 0.80 Alignments for each domain: == domain 1 score: 20.1 bits; conditional E-value: 2.4e-07 NodS 43 siadalevGcaaGalterlapyckrltvvdvveeaiararervkesshiswivsdvqqfstdekfd..livvaevlyyied 121 + ++le+G Galt l k++t v++ ++ + r++ + + + ++ dv +++ ++ + + vva++ yyi+ 23726_FusoPortal_Gene_900 32 DNDEILEIGPGQGALTSLLVERVKKVTCVEIDKDLENTLRKKFSSKENYTLVMGDVLEVDFKKYINqgTKVVANIPYYITS 112 34689*************************99988888899999999**********9998876543389********975 PP >> Methyltransf_23 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.4 0.0 1.8e-06 0.0028 20 62 .. 30 91 .. 8 239 .. 0.72 Alignments for each domain: == domain 1 score: 17.4 bits; conditional E-value: 1.8e-06 Methyltransf_23 20 lrargrvLdiGcGtGillellrengfsvtgvDlspeaveiae...................k 62 + ++L+iG+G+G l++ll e++++vt v+++ + ++ 23726_FusoPortal_Gene_900 30 ISDNDEILEIGPGQGALTSLLVERVKKVTCVEIDKDLENTLRkkfsskenytlvmgdvlevD 91 466779***************************99988887666666666555555555551 PP >> Microcin Colicin E1 (microcin) immunity protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.5 1.0 6.7e-06 0.01 22 48 .. 193 219 .. 192 224 .. 0.88 Alignments for each domain: == domain 1 score: 15.5 bits; conditional E-value: 6.7e-06 Microcin 22 ySkyaiEhialkfikkeFfekrknlns 48 y ++ E++++k++k F +krkn+ + 23726_FusoPortal_Gene_900 193 YKNKISEDLFFKYVKAAFSNKRKNIVN 219 667888******************876 PP >> PCMT Protein-L-isoaspartate(D-aspartate) O-methyltransferase (PCMT) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.6 0.0 2.6e-05 0.04 63 119 .. 22 75 .. 8 102 .. 0.77 2 ? -2.5 0.1 2.2 3.4e+03 2 29 .. 213 239 .. 212 242 .. 0.83 Alignments for each domain: == domain 1 score: 13.6 bits; conditional E-value: 2.6e-05 HHHHHHTT--TT-EEEEES-TTSHHHHHHHHHHCTTCEEEEEESBHHHHHHHHHHHH CS PCMT 63 alalellelkkgarvLevGsGsGyltavlarlvgekgrvvsierikelveeakrnlk 119 +++e ++++ +++Le+G G G lt +l + v++++ +e k+l ++ +++ 23726_FusoPortal_Gene_900 22 NKIIEVSNISDNDEILEIGPGQGALTSLLVERVKKVT---CVEIDKDLENTLRKKFS 75 56778888999************99999888888776...67766777666655554 PP == domain 2 score: -2.5 bits; conditional E-value: 2.2 HHHHHHHHHHHTTS---HHHHHHHHHS- CS PCMT 2 kkealieelkkkgviksdkvaealeevd 29 k+++++++l g +dk+ e+l++++ 23726_FusoPortal_Gene_900 213 KRKNIVNNLATLG-YSKDKIKEILNQIE 239 5688999999999.78899999999886 PP >> DREV DREV methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.7 0.1 6.1e-05 0.092 79 135 .. 22 73 .. 6 97 .. 0.73 Alignments for each domain: == domain 1 score: 11.7 bits; conditional E-value: 6.1e-05 DREV 79 rkllaiekespasdqedkklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkk 135 +k++ + s ++ ++l++G G+G +t + + +++v e+ + + + l+kk 23726_FusoPortal_Gene_900 22 NKIIEVS---NIS--DNDEILEIGPGQGALTSLLVERVKKVTCVEIDKDLENTLRKK 73 3444443...334..5669**************************999888888776 PP >> Methyltransf_12 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 0.0001 0.15 1 47 [. 37 81 .. 37 106 .. 0.87 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 0.0001 Methyltransf_12 1 LeiGcGtGsllralleklpdieytgvDispsaiedakeklakrennn 47 LeiG G G+l++ l e+ ++t+v+i+++ + ++k++++en+ 23726_FusoPortal_Gene_900 37 LEIGPGQGALTSLLVER--VKKVTCVEIDKDLENTLRKKFSSKENYT 81 8****************..**************99999999998755 PP >> UBA_2 Ubiquitin associated domain (UBA) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.0 2.2 1.2 1.8e+03 17 37 .. 17 37 .. 16 38 .. 0.72 2 ? -2.7 0.2 4.2 6.3e+03 35 43 .. 124 132 .. 114 133 .. 0.65 3 ? 11.2 0.1 0.00019 0.3 4 30 .. 216 241 .. 215 249 .. 0.79 Alignments for each domain: == domain 1 score: -1.0 bits; conditional E-value: 1.2 UBA_2 17 kdeIyeaLekaekskssneIk 37 kdeI++ + + + + + eI+ 23726_FusoPortal_Gene_900 17 KDEILNKIIEVSNISDNDEIL 37 99**99877666665777776 PP == domain 2 score: -2.7 bits; conditional E-value: 4.2 UBA_2 35 eIkDAYlLm 43 I +AY ++ 23726_FusoPortal_Gene_900 124 LIDEAYIMV 132 578888766 PP == domain 3 score: 11.2 bits; conditional E-value: 0.00019 UBA_2 4 dvvsaLsktMGYdkdeIyeaLekaeks 30 ++v+ L+ t+GY+kd I e L++ e 23726_FusoPortal_Gene_900 216 NIVNNLA-TLGYSKDKIKEILNQIEVP 241 6778885.9******999999976654 PP >> UBA_3 Fungal ubiquitin-associated domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.0 0.1 4.1 6.2e+03 24 28 .. 22 26 .. 14 27 .. 0.59 2 ? 11.6 0.8 0.00011 0.17 10 51 .. 215 256 .. 214 259 .. 0.92 Alignments for each domain: == domain 1 score: -3.0 bits; conditional E-value: 4.1 HHHHH CS UBA_3 24 dkvie 28 +k+ie 23726_FusoPortal_Gene_900 22 NKIIE 26 45555 PP == domain 2 score: 11.6 bits; conditional E-value: 0.00011 HHHHHHHHHHT--HHHHHHHHHHS--SS--SS--HHHHHHHH CS UBA_3 10 rdlVdefvsqGFekdkvievLkrLgiksldpednetankile 51 +++V++++ G+ kdk+ e+L ++++ + ++ +n+ +k +e 23726_FusoPortal_Gene_900 215 KNIVNNLATLGYSKDKIKEILNQIEVPENERAENISIDKFIE 256 689************************999999999999886 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (265 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1123 (0.0671972); expected 334.2 (0.02) Passed bias filter: 554 (0.0331498); expected 334.2 (0.02) Passed Vit filter: 63 (0.00376975); expected 16.7 (0.001) Passed Fwd filter: 11 (0.00065821); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 11 [number of targets reported over threshold] # CPU time: 0.31u 0.27s 00:00:00.58 Elapsed: 00:00:00.35 # Mc/sec: 2201.04 // [ok]