# 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_1029.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1029 [L=183] Description: # 1138551 # 1139099 # 1 # ID=1_1029;partial=00;start_type=ATG;rbs_motif=AGxAGG/AGGxGG;rbs_spacer=5-10bp;gc_cont=0.271 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-58 197.0 0.0 2.2e-58 196.8 0.0 1.0 1 Cons_hypoth95 Conserved hypothetical protein 95 4.4e-08 33.0 0.2 5.6e-08 32.6 0.2 1.1 1 Methyltransf_31 Methyltransferase domain 5.7e-07 29.1 0.0 6.5e-07 28.9 0.0 1.1 1 MTS Methyltransferase small domain 2.6e-06 27.2 0.5 0.028 14.0 0.1 2.3 2 MethyltransfD12 D12 class N6 adenine-specific DNA methyltran 3.2e-05 23.6 0.0 4.4e-05 23.2 0.0 1.2 1 Methyltransf_15 RNA cap guanine-N2 methyltransferase 5.4e-05 22.9 0.0 0.00084 19.0 0.0 2.1 1 UPF0020 Putative RNA methylase family UPF0020 0.00022 21.0 0.0 0.00029 20.6 0.0 1.3 1 Met_10 Met-10+ like-protein 0.00045 21.1 0.0 0.00068 20.5 0.0 1.2 1 Methyltransf_24 Methyltransferase domain 0.00045 20.0 0.0 0.00055 19.7 0.0 1.2 1 N6-adenineMlase Probable N6-adenine methyltransferase 0.00053 20.6 0.0 0.00073 20.1 0.0 1.3 1 Methyltransf_25 Methyltransferase domain 0.0015 18.2 0.0 0.016 14.8 0.0 2.1 2 N6_N4_Mtase DNA methylase 0.0034 16.6 0.0 0.0037 16.4 0.0 1.1 1 RsmJ Ribosomal RNA large subunit methyltransferas ------ inclusion threshold ------ 0.02 14.2 0.0 0.024 14.0 0.0 1.2 1 PrmA Ribosomal protein L11 methyltransferase (Prm 0.082 12.5 0.2 5.8 6.5 0.1 2.2 2 DUF4869 Domain of unknown function (DUF4869) Domain annotation for each model (and alignments): >> Cons_hypoth95 Conserved hypothetical protein 95 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 196.8 0.0 1.8e-61 2.2e-58 1 181 [. 1 180 [. 1 181 [. 0.99 Alignments for each domain: == domain 1 score: 196.8 bits; conditional E-value: 1.8e-61 EEXXSCCCTTXEEEXTTXXXXXSSSHHHHHHHHHHHCCCHCC-EEEETT-TTSHHHHHHHHTT-SEEEEEES-HHHHHHHHHH CS Cons_hypoth95 1 lriigGkakgrkLkvlegkstRPttdrvrealFnilaaeiegarvLDlfaGsGalglEAlSRgakkavlvekdkeavkilken 83 +rii+G+ak+r +k+ +g +tRPt ++v+e+lF+i+a++ieg+ +LDlf+GsG+++lEA+SRgak+av++ekd ea+k + en 23726_FusoPortal_Gene_1029 1 MRIIAGEAKNRIIKTRKGFDTRPTLESVKESLFSIIAPYIEGSIFLDLFSGSGSISLEAISRGAKRAVMIEKDGEALKYIIEN 83 9********************************************************************************** PP HCCCT--TTTEEEEES-HHHHCTS--SS-EEEEEEE--SSSTTCH.HHHHHHHHHTT-EEEEEEEEEEEECCSSS--.-TTEE CS Cons_hypoth95 84 leelkleeetaelkedaleallklakkkeefdivflDPPyakellieevlelleekkllkenalvvvetekekelpespenle 166 +++l+++++++++k+d +a++ l++k+e+fdi+f+DPPy+++ + ++vl+++ + ++l+e++l+++e++ + l ++++ ++ 23726_FusoPortal_Gene_1029 84 IDNLGFSDRCRAYKNDVIRAIEILGRKNEKFDIIFMDPPYQDNVT-KKVLKAIDKANILAEDGLIICEHHLLEDLEDNIASFR 165 *******************************************97.************************************* PP EEEEEEETTEEEEEE CS Cons_hypoth95 167 lvkekkygettlaly 181 +++e+ky++++l++y 23726_FusoPortal_Gene_1029 166 KTDERKYNKKILTFY 180 **************9 PP >> Methyltransf_31 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 32.6 0.2 4.7e-11 5.6e-08 7 124 .. 45 163 .. 41 180 .. 0.80 Alignments for each domain: == domain 1 score: 32.6 bits; conditional E-value: 4.7e-11 Methyltransf_31 7 vlDlGCGtGeltfelaeelgpkaevvGiDiseeaielArenaqklgld.nveFeqgdiee.lekl.ledekfDlilsnevl.. 84 +lDl +G+G +++e+ ++ ++ v i+ + ea+++ en ++lg++ +++ + d+ +e l + ekfD+i+ + + 23726_FusoPortal_Gene_1029 45 FLDLFSGSGSISLEAISR--GAKRAVMIEKDGEALKYIIENIDNLGFSdRCRAYKNDVIRaIEILgRKNEKFDIIFMDPPYqd 125 79**************65..346699*********************99***********999988888********866655 PP Methyltransf_31 85 naladpdkileeiervlkpggilliedprllaelpakvee 124 n ++ k++ + + l ++g +++e++ l++l +++ + 23726_FusoPortal_Gene_1029 126 NVTKKVLKAIDK-ANILAEDGLIICEHHL-LEDLEDNIAS 163 555555555556.6789999999988887.7777654444 PP >> MTS Methyltransferase small domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 28.9 0.0 5.5e-10 6.5e-07 35 133 .. 45 145 .. 14 174 .. 0.83 Alignments for each domain: == domain 1 score: 28.9 bits; conditional E-value: 5.5e-10 EEEET-TTSHHHHHHHHHTTC.EEEEEESBHHHHHHHHHHHHHTTGC..EEEEE-STTTTST.....TS-EEEEEE---STTT CS MTS 35 vlDlGcGaGvlgaalakespeaevtlvdvdarAlesakenleankle..vevvasdvleeve.....ekkfdlivsNPPfhag 110 lDl +G G +++++ ++ + + +++++d +Al+ en+ + +++ ++++++dv+ ++e ++kfd+i + PP++++ 23726_FusoPortal_Gene_1029 45 FLDLFSGSGSISLEAISRGAK-RAVMIEKDGEALKYIIENIDNLGFSdrCRAYKNDVIRAIEilgrkNEKFDIIFMDPPYQDN 126 69***************9988.789*********************999*********98865555599***********988 PP CHCHHHHHHHHHHHHHHHEEEEE CS MTS 111 katelavaeeliaeAkeaLkegG 133 ++v ++ i +A ++L+e+G 23726_FusoPortal_Gene_1029 127 VT--KKVLKA-IDKA-NILAEDG 145 76..444444.3333.5666666 PP >> MethyltransfD12 D12 class N6 adenine-specific DNA methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.3 0.0 0.00016 0.19 5 41 .. 29 62 .. 28 71 .. 0.83 2 ! 14.0 0.1 2.4e-05 0.028 149 189 .. 85 125 .. 71 139 .. 0.79 Alignments for each domain: == domain 1 score: 11.3 bits; conditional E-value: 0.00016 SCTTHHHHHCCS-..--SEEEETT-TTSCCHHCS... CS MethyltransfD12 5 KrkllpeIkelipkdkgktfvepFaGggaVffeakkq 41 K++l + I +i+ g+ f ++F+G+g++++ea ++ 23726_FusoPortal_Gene_1029 29 KESLFSIIAPYIE---GSIFLDLFSGSGSISLEAISR 62 7888888888888...689*************99764 PP == domain 2 score: 14.0 bits; conditional E-value: 2.4e-05 HHHHHHHHHHHCCEEEC-CHHHHTT--TTE.EEEE--S-T. CS MethyltransfD12 149 elklrseklknakiecgdfdavllkakieddivYlDPPYnp 189 + +s + + +k+++ ++ ++l +++++ di+++DPPY+ 23726_FusoPortal_Gene_1029 85 DNLGFSDRCRAYKNDVIRAIEILGRKNEKFDIIFMDPPYQD 125 445566777779999999999999***************97 PP >> Methyltransf_15 RNA cap guanine-N2 methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 23.2 0.0 3.7e-08 4.4e-05 3 81 .. 44 123 .. 42 139 .. 0.91 Alignments for each domain: == domain 1 score: 23.2 bits; conditional E-value: 3.7e-08 Methyltransf_15 3 vvlDlfgGaGGnviq.fanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPw 81 ++lDlf+G+G + + + + ++ + i+++ + lk + +N g sd+ + dv +++ l ++ k+d++f+ pP 23726_FusoPortal_Gene_1029 44 IFLDLFSGSGSISLEaISRGAKRAVMIEKDGEALKYIIENIDNLGFSDRCRAYKNDVIRAIEILGRKNEKFDIIFMDPPY 123 78********86654145558999****************************************999999*********5 PP >> UPF0020 Putative RNA methylase family UPF0020 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.0 0.0 7.1e-07 0.00084 27 138 .. 40 130 .. 31 137 .. 0.88 Alignments for each domain: == domain 1 score: 19.0 bits; conditional E-value: 7.1e-07 -TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SSHHHHHHHHHHHH--EEEEES-HHHHHHHHHHHHHTT-GGGEEEEE CS UPF0020 27 kpgepllDPmCGsGtilIEaallganiapgllrefvlelkaeaeeearaelklygsDldrrvvqgareNaekagvgdliefsq 109 +g+ +lD ++GsG+i +Ea ga++a + +++d ++++ eN + g +d+ +++ 23726_FusoPortal_Gene_1029 40 IEGSIFLDLFSGSGSISLEAISRGAKRA-------------------------VMIEKDGEALKYIIENIDNLGFSDRCRAYK 97 57899***********************.........................679*************************99 PP --G.GG--...STTS---EEEE----------H CS UPF0020 110 ada.akLr...lkegevdvivtnpPYGerlgsk 138 d+ ++++ +k++++d+i +pPY + +++k 23726_FusoPortal_Gene_1029 98 NDViRAIEilgRKNEKFDIIFMDPPYQDNVTKK 130 997345667789999***********9887765 PP >> Met_10 Met-10+ like-protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.6 0.0 2.4e-07 0.00029 96 177 .. 37 122 .. 1 149 [. 0.77 Alignments for each domain: == domain 1 score: 20.6 bits; conditional E-value: 2.4e-07 TT--TT-EEEETT-TTTCCHHHCCHHTTSSEEEEEES-HHHHHHHHHHHHHTT-TTCEEEEES-GGG---.....TT-EEEEE CS Met_10 96 klveegevVvdvfAGiGpfaipiakyskakrvlaidlnPesvkylkeniklNkvekvvkpileDgrevir.....eqladrvv 173 + eg + +d+f G G ++ + + akr + i+++ e+ ky+ eni +++ +++++D+ i+ ++++d + 23726_FusoPortal_Gene_1029 37 APYIEGSIFLDLFSGSGSISLEAISRG-AKRAVMIEKDGEALKYIIENIDNLGFSDRCRAYKNDVIRAIEilgrkNEKFDIIF 118 44568999*************887755.555677******************************9776655666667777777 PP E--T CS Met_10 174 mnlp 177 m+ p 23726_FusoPortal_Gene_1029 119 MDPP 122 7655 PP >> Methyltransf_24 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.5 0.0 5.7e-07 0.00068 15 103 .. 57 149 .. 48 151 .. 0.80 Alignments for each domain: == domain 1 score: 20.5 bits; conditional E-value: 5.7e-07 Methyltransf_24 15 lralrrngagrlvavdhdpgaeeaaalare...aglddrvrlikgdseealpsla..agsiDllfiDgghaaeavledleal. 91 l+a+ r ++ r v++++d e+ + ++++ +g++dr r k+d a++ l ++++D++f+D + + + +++ l+a+ 23726_FusoPortal_Gene_1029 57 LEAISRGAK-RAVMIEKDG--EALKYIIENidnLGFSDRCRAYKNDVIRAIEILGrkNEKFDIIFMDPPYQDNVTKKVLKAId 136 566777775.899999999..344444444777**************888887777679********9999999999999995 PP Methyltransf_24 92 .lprlapggvvlf 103 ++la++g++++ 23726_FusoPortal_Gene_1029 137 kANILAEDGLIIC 149 5579*******99 PP >> N6-adenineMlase Probable N6-adenine methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.7 0.0 4.6e-07 0.00055 84 150 .. 110 175 .. 78 182 .. 0.69 Alignments for each domain: == domain 1 score: 19.7 bits; conditional E-value: 4.6e-07 N6-adenineMlase 84 lkgkfdviviDPPF...lseeclkklaetaklllk.klilctgevmenlekllgvkktkfepehenklsne 150 ++kfd+i +DPP+ +++++lk ++++ l + li+c ++ e+le+++ ++f+ + e+k +++ 23726_FusoPortal_Gene_1029 110 KNEKFDIIFMDPPYqdnVTKKVLKAIDKANILA-EdGLIICEHHLLEDLEDNI----ASFRKTDERKYNKK 175 489***********8877777777777744444.4479999988888766544....45555555544443 PP >> Methyltransf_25 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.1 0.0 6.1e-07 0.00073 2 73 .. 46 121 .. 45 145 .. 0.80 Alignments for each domain: == domain 1 score: 20.1 bits; conditional E-value: 6.1e-07 EEET-TTSHHHHHHHCTSTSSEEEEEESSHHHHHHHHHHHHHTTT..TEEEEESBC.CG...BSSSTTTEEEEEEEST CS Methyltransf_25 2 LdlGcGtGrltlalakrlgkakvtgvDiseemleaareraeeegl..nvefvqgda.ed...lpfedgsfDlvvssga 73 Ldl G+G ++l++ r g+++ ++++ + e+l+++ e++++ g+ +++ ++ d+ ++ l ++++fD++++ 23726_FusoPortal_Gene_1029 46 LDLFSGSGSISLEAISR-GAKRAVMIEKDGEALKYIIENIDNLGFsdRCRAYKNDViRAieiLGRKNEKFDIIFMD-P 121 9****************.69********************99999******999997552225689********93.3 PP >> N6_N4_Mtase DNA methylase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.8 0.0 1.3e-05 0.016 190 229 .. 40 80 .. 4 83 .. 0.84 2 ? 1.5 0.0 0.16 1.9e+02 2 15 .. 115 128 .. 114 177 .. 0.74 Alignments for each domain: == domain 1 score: 14.8 bits; conditional E-value: 1.3e-05 -TT-EEEETT-TTTHHHHHHHHTT-E.EEEEESSHHHHHHH CS N6_N4_Mtase 190 kkgdivlDpFaGSGtTavaakelgrk.fiGieieeeyvela 229 +g+i lD F GSG++ +a++ g k ++ ie++ e+++ + 23726_FusoPortal_Gene_1029 40 IEGSIFLDLFSGSGSISLEAISRGAKrAVMIEKDGEALKYI 80 469***************999998766999***99998854 PP == domain 2 score: 1.5 bits; conditional E-value: 0.16 EEEEE--SSSSSCS CS N6_N4_Mtase 2 dlivtdpPYntgkd 15 d+i+ dpPY + 23726_FusoPortal_Gene_1029 115 DIIFMDPPYQDNVT 128 789999**986654 PP >> RsmJ Ribosomal RNA large subunit methyltransferase D, RlmJ # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.4 0.0 3.1e-06 0.0037 62 161 .. 46 148 .. 5 179 .. 0.72 Alignments for each domain: == domain 1 score: 16.4 bits; conditional E-value: 3.1e-06 --EEE-HHHHHHH.HS-TT-EEEEE---TTHHHHHHHHTTT...-TTEEEESS-HHHHHHHS-S-TTS-EEEEE----SSCCH CS RsmJ 62 lryYpGspllarr.llReqDrlvlteLhpkdaalLkaelak...dkrvrvkkadGfaalkaklppkekRglvliDPPyeakde 140 l++++Gs + + + R + r+v++e + + + +++ + r r++k+d a++ ++ek ++++DPPy + 23726_FusoPortal_Gene_1029 46 LDLFSGSGSISLEaISRGAKRAVMIEKDGEALKYIIENIDNlgfSDRCRAYKNDVIRAIEILGRKNEKFDIIFMDPPYQDN-V 127 78999999997651679********9777777766655443222789*********99999999999***********755.5 PP HHHHHHHHHHHHHH-TTSEEE CS RsmJ 141 ydrvvaalrealkRlaqavla 161 ++v++a+++a + ++ + 23726_FusoPortal_Gene_1029 128 TKKVLKAIDKANILAEDGLII 148 68999**99998777777655 PP >> PrmA Ribosomal protein L11 methyltransferase (PrmA) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.0 0.0 2e-05 0.024 161 231 .. 41 118 .. 24 167 .. 0.82 Alignments for each domain: == domain 1 score: 14.0 bits; conditional E-value: 2e-05 PrmA 161 egeevlDvGtGSGiLaiaaakLgakkvvavDiDevAvkvakeNaelNkveaklevalseell........keekadviv 231 eg+ +lD+ GSG +++ a gak++v + D A+k eN++ + +++++ + +++++ k+ek+d+i+ 23726_FusoPortal_Gene_1029 41 EGSIFLDLFSGSGSISLEAISRGAKRAVMIEKDGEALKYIIENIDNLGFSDRCRAY-KNDVIraieilgrKNEKFDIIF 118 567889***************************************99999999774.4555533434444489999996 PP >> DUF4869 Domain of unknown function (DUF4869) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 3.8 0.0 0.033 39 46 73 .. 53 80 .. 46 97 .. 0.80 2 ? 6.5 0.1 0.0048 5.8 8 38 .. 114 144 .. 107 157 .. 0.83 Alignments for each domain: == domain 1 score: 3.8 bits; conditional E-value: 0.033 DUF4869 46 GkispkklsgGvktLlllanedekifnv 73 G+is + +s G+k ++++++ e + + 23726_FusoPortal_Gene_1029 53 GSISLEAISRGAKRAVMIEKDGEALKYI 80 89****************9998865444 PP == domain 2 score: 6.5 bits; conditional E-value: 0.0048 DUF4869 8 fkntYedeWltdelakeiikdiDksevlsee 38 f+ + d +d+++k+++k iDk+++l+e+ 23726_FusoPortal_Gene_1029 114 FDIIFMDPPYQDNVTKKVLKAIDKANILAED 144 666777777899*************998764 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (183 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 437 (0.0261489); expected 334.2 (0.02) Passed bias filter: 325 (0.0194471); expected 334.2 (0.02) Passed Vit filter: 36 (0.00215414); expected 16.7 (0.001) Passed Fwd filter: 14 (0.000837721); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 14 [number of targets reported over threshold] # CPU time: 0.26u 0.25s 00:00:00.51 Elapsed: 00:00:00.33 # Mc/sec: 1612.08 // [ok]