# hmmsearch :: search profile(s) against a sequence 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 HMM file: pfam_list/Methyltransf_5.hmm.txt
# target sequence database: proteomes/Drosophila_melanogaster.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Methyltransf_5 [M=309]
Accession: PF01795.18
Description: MraW methylase family
Scores for complete sequences (score includes all domains):
--- full sequence --- --- best 1 domain --- -#dom-
E-value score bias E-value score bias exp N Sequence Description
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
1.7e-82 277.0 0.0 2.5e-82 276.5 0.0 1.2 1 sp|Q9VGY5|MET15_DROME Probable methyltransferase-like prote
0.00012 21.4 0.0 0.00031 20.0 0.0 1.6 2 sp|Q9VIK9|MET13_DROME Methyltransferase-like protein 13 OS=
0.0095 15.1 0.0 0.025 13.7 0.0 1.6 2 tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransfe
------ inclusion threshold ------
0.12 11.5 0.0 0.16 11.0 0.0 1.2 1 sp|Q9VZD2|BN3D3_DROME Probable RNA methyltransferase CG1134
Domain annotation for each sequence (and alignments):
>> sp|Q9VGY5|MET15_DROME Probable methyltransferase-like protein 15 homolog OS=Drosophila melanogaster OX=7227 GN=CG146
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 276.5 0.0 7.2e-86 2.5e-82 2 308 .. 25 354 .. 24 355 .. 0.88
Alignments for each domain:
== domain 1 score: 276.5 bits; conditional E-value: 7.2e-86
XXX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.......---EEEEES-GGGHHHHHHH CS
Methyltransf_5 2 leHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalaraker..Lkefekrltlvkgnfkelkealae 87
+ HvpVL + ++e+l+p ++g+y D+T+GaggH+++lLe+ +e +++++drD+ a + a+++ +ef+ rl + g+f++l + +e
sp|Q9VGY5|MET15_DROME 25 APHVPVLCDTAIEYLQPVPGGTYFDMTFGAGGHTRRLLEKCPEAKVYALDRDPLAHQLARDMseSEEFKGRLIPLLGKFSDLPKLFKE 112
57***************************************************999999875115689999***************** PP
TT..TT-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-..TTSS--HHHHHHH--HHHHHHHHHHHH--TTHHHHHHHHHHHHH CS
Methyltransf_5 88 an..vtkvdgillDLGvSslqldkaerGFsfaedapLDmRmd..qsaeltAaevlnrlseeeLakilkeygeekfakriakaivekrk 171
++ ++vdg+l+D G Ss+q+d+a rGFs +d+pLDmRmd +s +tAa+vl +++e +L+kil+ ygeek ak+ia+ +v++r+
sp|Q9VGY5|MET15_DROME 113 HGlaKNSVDGMLFDFGCSSMQFDEAVRGFSISRDGPLDMRMDggHSGGVTAADVLANVEEGDLVKILRMYGEEKAAKKIARGLVDARN 200
**87789***********************************3335567**********************************99665 PP
.HS---BHHHHHHHHHHHS-------......-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHHH CS
Methyltransf_5 172 .kkpiettkelaevikkvvpakekrk....kihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkkvfk 254
+iettk+la++i++++ +++ h at++FQA+RI+vN+El+e++ + a e+l+ +grl+ i+FhSled ivk++++
sp|Q9VGY5|MET15_DROME 201 aLFKIETTKQLADIIENIMDGGTRKDklrrPAHSATKTFQAIRIFVNNELNEINYGMVLANEILRVDGRLVTITFHSLEDTIVKRHIN 288
1579*************998665554233368*******************************************************9 PP
HTT------------------.............-EEEEE.E.E---HHHHHH-GGGTT-EEEEEEE CS
Methyltransf_5 255 easkkeevPaglpvaseelka.............klklitkkpilpseeeikeNprsrsakLrviek 308
+ + lp++ +++ a k + + +i+p ++e+++N rsrsakLr++ k
sp|Q9VGY5|MET15_DROME 289 GNVL-AGAANALPLKYSSHYAidepdilesltkkSWKQLHRHVIVPDADEVARNTRSRSAKLRAAVK 354
9887.46677777777666555556644444444667777779********************9865 PP
>> sp|Q9VIK9|MET13_DROME Methyltransferase-like protein 13 OS=Drosophila melanogaster OX=7227 GN=CG2614 PE=2 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -3.9 0.0 1.7 6e+03 212 244 .. 341 373 .. 338 379 .. 0.79
2 ! 20.0 0.0 9.1e-08 0.00031 27 100 .. 475 548 .. 451 559 .. 0.88
Alignments for each domain:
== domain 1 score: -3.9 bits; conditional E-value: 1.7
HHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHH CS
Methyltransf_5 212 vNdEleeleelleaalellakggrlavisFhSl 244
v + l+e+++ l+ ++++l + g i + Sl
sp|Q9VIK9|MET13_DROME 341 VYSTLDEVKQELADSIKNLSPAGLTDQIPYLSL 373
667789999999999999999988877777776 PP
== domain 2 score: 20.0 bits; conditional E-value: 9.1e-08
S--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---.EEEEES-GGGHHHHHHHTTTT-EEEEEEE- CS
Methyltransf_5 27 aTlGaggHsealLeqlselrligidrDksalarakerLkefekr.ltlvkgnfkelkealaeanvtkvdgillDL 100
+ lG+gg l +l++ r++++++D+ +l++a++ ++ ++++ +++v + ++ e ++++++ ++d++l+D
sp|Q9VIK9|MET13_DROME 475 VGLGGGGLCSFLHAALPQARITAVEIDPIMLEVAEQYFELKQDKrFHVVIDDGLDFVERCRNEDI-HFDAVLFDV 548
5699*****999999**********************99888666****999999******9999.9******95 PP
>> tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransferase catalytic subunit TRMT61A OS=Drosophila melanogaste
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -2.8 0.0 0.77 2.6e+03 89 106 .. 9 26 .. 5 30 .. 0.83
2 ! 13.7 0.0 7.4e-06 0.025 12 100 .. 103 192 .. 99 196 .. 0.84
Alignments for each domain:
== domain 1 score: -2.8 bits; conditional E-value: 0.77
TTT-EEEEEEE-S--HHH CS
Methyltransf_5 89 nvtkvdgillDLGvSslq 106
+++k d ++l L vSs++
tr|Q9VBH7|Q9VBH7_DROME 9 HIEKGDVVILYLSVSSMH 26
4889999*********98 PP
== domain 2 score: 13.7 bits; conditional E-value: 7.4e-06
HHHHT---TT-EEEES--TTSHHHHHHHTT--TT.EEEEEES-HHHHHHHTT.....---..EEEEES-GGGHHHHHHHTTTT-EEE CS
Methyltransf_5 12 vvellepkedgvyiDaTlGaggHsealLeqlsel.rligidrDksalarakerLkefekr..ltlvkgnfkelkealaeanvtkvdg 95
+++ le+++++v + + G+g s+++L++l+ + +l ++d + ++a++ ++++ +t+++ +l + +k+d+
tr|Q9VBH7|Q9VBH7_DROME 103 IIHQLEVRPGAVVVESGTGSGSLSHYFLRALKPTgHLHTFDFHEARADQARDEFRRHGIAdfVTVYHRDVCNLGFGAELD--GKADA 187
7999***************************9999*****************9999874335788887776666554433..57888 PP
EEEE- CS
Methyltransf_5 96 illDL 100
++lDL
tr|Q9VBH7|Q9VBH7_DROME 188 VFLDL 192
88887 PP
>> sp|Q9VZD2|BN3D3_DROME Probable RNA methyltransferase CG11342 OS=Drosophila melanogaster OX=7227 GN=CG11342 PE=3 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 11.0 0.0 4.8e-05 0.16 40 116 .. 75 157 .. 48 166 .. 0.77
Alignments for each domain:
== domain 1 score: 11.0 bits; conditional E-value: 4.8e-05
TT--TT.EEEEEES-HHHHHHHTT.....---....EEEEE.S-GGGHHHHHHHTTTT-EEEEEEE-S--HHHHH-GGG---S CS
Methyltransf_5 40 eqlsel.rligidrDksalarakerLkefekr....ltlvk.gnfkelkealaeanvtkvdgillDLGvSslqldkaerGFsf 116
e+l++ +++g+d+D++ ++ra e ++ ++ + ++ + f+++k++++ +n +k+d+i++ ++l+++++G f
sp|Q9VZD2|BN3D3_DROME 75 ERLHRSvKVLGVDIDPRLIQRASEENESPKDVsyacVDVLDdEAFESVKTYMEVNNLEKFDAICCYSITMWIHLNHHDQGLRF 157
445555699*************99776665545555555551568******************98665566788888888776 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (309 nodes)
Target sequences: 13786 (7394964 residues searched)
Passed MSV filter: 459 (0.0332946); expected 275.7 (0.02)
Passed bias filter: 234 (0.0169737); expected 275.7 (0.02)
Passed Vit filter: 22 (0.00159582); expected 13.8 (0.001)
Passed Fwd filter: 4 (0.000290149); expected 0.1 (1e-05)
Initial search space (Z): 13786 [actual number of targets]
Domain search space (domZ): 4 [number of targets reported over threshold]
# CPU time: 0.17u 0.01s 00:00:00.18 Elapsed: 00:00:00.06
# Mc/sec: 38084.06
//
[ok]