# 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/Methyltr_RsmB-F.hmm.txt
# target sequence database: proteomes/Drosophila_melanogaster.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Methyltr_RsmB-F [M=200]
Accession: PF01189.16
Description: 16S rRNA methyltransferase RsmB/F
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.4e-70 236.7 0.3 4.1e-70 235.1 0.0 2.0 2 tr|A1Z909|A1Z909_DROME LD11307p OS=Drosophila melanogaster O
6.9e-39 133.2 0.0 1.2e-38 132.4 0.0 1.4 1 tr|Q8IPS4|Q8IPS4_DROME FI01826p OS=Drosophila melanogaster O
7.2e-36 123.3 0.0 1.1e-34 119.5 0.0 2.2 1 sp|Q9W4M9|NSUN2_DROME tRNA (cytosine(34)-C(5))-methyltransf
1.2e-30 106.3 0.1 2.4e-30 105.3 0.1 1.5 1 tr|Q9VDZ4|Q9VDZ4_DROME LD16340p OS=Drosophila melanogaster O
3.4e-27 95.0 0.0 6.1e-27 94.2 0.0 1.5 1 tr|Q9VPX3|Q9VPX3_DROME LD40271p OS=Drosophila melanogaster O
------ inclusion threshold ------
0.02 14.2 0.0 0.049 12.9 0.0 1.6 1 sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillar
Domain annotation for each sequence (and alignments):
>> tr|A1Z909|A1Z909_DROME LD11307p OS=Drosophila melanogaster OX=7227 GN=CG8545-RA PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 235.1 0.0 1.8e-73 4.1e-70 1 200 [] 428 636 .. 428 636 .. 0.94
2 ? -2.6 0.3 1.2 2.8e+03 143 185 .. 649 690 .. 647 699 .. 0.66
Alignments for each domain:
== domain 1 score: 235.1 bits; conditional E-value: 1.8e-73
HHHH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC--EEEEEEEE CS
Methyltr_RsmB-F 1 allLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlggelfDriLlDa 87
+++Lapqe+e+ilD+c+aPGgK +hia +mkn+g+++A+d n+ R+k++++n +rLg+ n+ v++ D++k+ ++++g fDr+LlDa
tr|A1Z909|A1Z909_DROME 428 VMALAPQENERILDMCSAPGGKGSHIASIMKNTGVLFANDSNRDRIKAIVANFHRLGIVNAVVSCEDGTKFRNIMTG--FDRVLLDA 512
589*************************************************************************9..******** PP
--CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHE....EEEEEEEEEES---CTCTHHHHHHHHHH-TTEEEESSBSTTS CS
Methyltr_RsmB-F 88 PCsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidll....kkgGvlvYsTCsvlveEneevieeaLkkredvelvetelsep 170
PC+gtgv+++dp+vk++++e d++++ +lq++ll +aid + +gG++vYsTCsvl+eEne+vi++aLkkr +v+lv+t l +
tr|A1Z909|A1Z909_DROME 513 PCTGTGVVSKDPSVKTTKSEVDVQRCYNLQRKLLLTAIDCTdaksSTGGYIVYSTCSVLPEENEWVIDYALKKR-NVKLVPTGLDFG 598
***************************************99998789**************************9.9******66544 PP
TTSSECEEC........T-EEE-TTTTS.-SEEEEEEE CS
Methyltr_RsmB-F 171 kglklgeea........gtlrllPhtenkdgffiaklk 200
++ ++ + t+r++Pht+n+dgf++aklk
tr|A1Z909|A1Z909_DROME 599 VEGFTKFRQhrfhpslnLTKRYYPHTHNMDGFYVAKLK 636
333333333689999999******************96 PP
== domain 2 score: -2.6 bits; conditional E-value: 1.2
CTCTHHHHHHHHHH-TTEEEESSBSTTSTTSSECEECT-EEE- CS
Methyltr_RsmB-F 143 veEneevieeaLkkredvelvetelsepkglklgeeagtlrll 185
eE+e++++ea+++ + +e +e ++ ++e+g ++ll
tr|A1Z909|A1Z909_DROME 649 QEEDEKQLDEAIEADATTEAAAE-SEEVEEEVASKEKGPRKLL 690
68999****99999877776554.4454555555555556655 PP
>> tr|Q8IPS4|Q8IPS4_DROME FI01826p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG11109 PE=2 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 132.4 0.0 5.3e-42 1.2e-38 2 200 .] 221 437 .. 220 437 .. 0.77
Alignments for each domain:
== domain 1 score: 132.4 bits; conditional E-value: 5.3e-42
HHH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC........... CS
Methyltr_RsmB-F 2 llLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg........... 77
+L+pq++e+ilD+caaPG KT+hiaelm ++g vvA+d++++R++++ l + v++ + +
tr|Q8IPS4|Q8IPS4_DROME 221 RVLDPQPGERILDMCAAPGNKTSHIAELMGDRGSVVALDNSASRVRSMLPKLGHYKSITAHVFNSTEAVAADAPSApvgeftgppfp 307
58*******************************************999776655544444444433222222222223347779*** PP
.--EEEEEEEE--CCGGGTTT-TTHH.HH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-TTEEE CS
Methyltr_RsmB-F 78 .elfDriLlDaPCsgtgvirrdpdvk.lrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkredvel 162
e fDriLlDaPCsg+g +p++ ++ k ++++ + q++l++ a++l+++gG+lvYsTC+v+ +E e ++ aL+k + +l
tr|Q8IPS4|Q8IPS4_DROME 308 cESFDRILLDAPCSGLGN---RPQLScSIKQVKVLSSYPHIQRRLFEQAVQLVRPGGILVYSTCTVTEDECECLVAWALRKFVELRL 391
***************995...344430457888999**********************************************99999 PP
ESSB..STTSTTSSECEEC....T-EEE-TTTTS.-..SEEEEEEE CS
Methyltr_RsmB-F 163 vete..lsepkglklgeea....gtlrllPhtenkd..gffiaklk 200
++ + +p g ea +r+ P n d gffiak++
tr|Q8IPS4|Q8IPS4_DROME 392 TDATprWGGPGLPLPGFEAcksrLLQRFGPSGANADtvGFFIAKFQ 437
8753113333333344444555477788899888876689999985 PP
>> sp|Q9W4M9|NSUN2_DROME tRNA (cytosine(34)-C(5))-methyltransferase OS=Drosophila melanogaster OX=7227 GN=Nsun2 PE=2 SV
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 119.5 0.0 4.6e-38 1.1e-34 2 199 .. 170 429 .. 169 430 .. 0.86
Alignments for each domain:
== domain 1 score: 119.5 bits; conditional E-value: 4.6e-38
HHH---TT-EEEECS-TT-HHHHHHHHHTTTT........SEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC.... CS
Methyltr_RsmB-F 2 llLapqeeetilDlcaaPGgKTthiaelmknk........gtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg.... 77
++L+ ++ +++lD+caaPG+KT+ + e ++ + g v+A+d++++R +++ krL+ + v++ D + +++ ++
sp|Q9W4M9|NSUN2_DROME 170 IVLDVRPTDKVLDMCAAPGSKTAQLIEALHAApeehkippGFVLANDVDNNRCYMLVHQAKRLNSPCLLVTNHDSSVFPNLVTTkpdg 257
7899*******************99887543212223333789***********************************9999888999 PP
....--EEEEEEEE--CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-T.TE CS
Methyltr_RsmB-F 78 ....elfDriLlDaPCsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkre.dv 160
+fD+iL D+PCsg g++r++pd+ l+ + ++ +l+ q +++ +++l gG lvYsTCs+++ Ene+v+++++k+ +
sp|Q9W4M9|NSUN2_DROME 258 skaiLKFDKILCDVPCSGDGTLRKNPDIWLKWNLAQAYNLHGIQYRIVRRGAEMLEVGGRLVYSTCSLNPIENEAVLQRIIKDADgAL 345
999889**********************7788889999*******************************************9986234 PP
EEESSBSTTSTTSSECEEC..............................................T-EEE-TTTTS.-SEEEEEE CS
Methyltr_RsmB-F 161 elvetelsepkglklgeea..............................................gtlrllPhtenkdgffiakl 199
elv+ glk++ + + lr+lPh ++ gff+a l
sp|Q9W4M9|NSUN2_DROME 346 ELVDAGHL-VPGLKYKPGMtdwklatkevdqiftrfeevpeslhtiirpgmfplpademakigleKCLRVLPHLQDSGGFFVAVL 429
66666332.2344444444559999999999999999999999999999999999999999999999999999999999999987 PP
>> tr|Q9VDZ4|Q9VDZ4_DROME LD16340p OS=Drosophila melanogaster OX=7227 GN=Nsun5 PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 105.3 0.1 1e-33 2.4e-30 3 199 .. 222 424 .. 220 425 .. 0.90
Alignments for each domain:
== domain 1 score: 105.3 bits; conditional E-value: 1e-33
HH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCC.HHHHCC--EEEEEEEE- CS
Methyltr_RsmB-F 3 lLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarki.dkvlggelfDriLlDaP 88
lLap + t+lD+caaPG KT hi ++m+nkg+++ v++++ R + + e +k g+ ++ + +Da + ++ ++ iL+D
tr|Q9VDZ4|Q9VDZ4_DROME 222 LLAPPSGATVLDMCAAPGMKTVHICNVMQNKGCIYSVEQDHVRYNTLCEITKDAGCDIVKPILGDALNLtPERFP--DVEYILVDPS 306
799******************************************************999999998776144444..599******* PP
-CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-TTEEEESSBST.TSTTSS CS
Methyltr_RsmB-F 89 CsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkredvelvetels.epkglk 174
Csg g+ +r +++k +++l+ lq ++ls a+ + + Y+TCs+ eEne+v+++ L+ ++ +l + + k +
tr|Q9VDZ4|Q9VDZ4_DROME 307 CSGSGMQNRMTVCDEPKEDKRLQKLQGLQIKILSHAMGAFPNVKRIAYCTCSLWKEENEQVVQRCLQLNPSFKLLSCKKAlRNKWHN 393
****************************************9999*****************************99886441445566 PP
ECEEC......T-EEE-TTTTS.-SEEEEEE CS
Methyltr_RsmB-F 175 lgeea......gtlrllPhtenkdgffiakl 199
+g + + l P+++ +dg+f+a +
tr|Q9VDZ4|Q9VDZ4_DROME 394 VGDKDypnigkNVLYCQPDSDLTDGIFLALF 424
6666689987777888899*********987 PP
>> tr|Q9VPX3|Q9VPX3_DROME LD40271p OS=Drosophila melanogaster OX=7227 GN=l(2)10685 PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 94.2 0.0 2.6e-30 6.1e-27 3 158 .. 285 446 .. 283 498 .. 0.83
Alignments for each domain:
== domain 1 score: 94.2 bits; conditional E-value: 2.6e-30
HH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHH........HTT--SEEEEES-CCCCHHHHCC--EE CS
Methyltr_RsmB-F 3 lLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlk........rLgvtntivleaDarkidkvlggelfD 81
+L+ +++e++lD+ca+PGgK ++++ + ++v +d++++Rl+++ + ++ r + +i+++ Dar+ d++ fD
tr|Q9VPX3|Q9VPX3_DROME 285 FLQVKPGERVLDACASPGGKSLLMLQT-LHLDHLVCNDIQESRLNKLRKVMQeylfdykeRWAGKRLILSQSDARNLDQYEQ---FD 367
79999********************99.6889999**********9987765222222226677799**********99866...** PP
EEEEEE--CCGGGTTT-TT..HHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-T CS
Methyltr_RsmB-F 82 riLlDaPCsgtgvirrdpd..vklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkre 158
+iL+D+PC++ + ++d + + k+ ++ +lq +l + ll++gG lvYsTCs+ + n v++ aL+k
tr|Q9VPX3|Q9VPX3_DROME 368 KILVDVPCTTDRHVLNEQDnnIFKPTRIKERLRIPELQAGILANCLRLLRPGGSLVYSTCSLSPIQNDGVVHMALQKVF 446
*********9765554444015444455677889******************************************965 PP
>> sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillarin OS=Drosophila melanogaster OX=7227 GN=Fib PE=2 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 12.9 0.0 2.1e-05 0.049 10 87 .. 186 262 .. 179 278 .. 0.77
Alignments for each domain:
== domain 1 score: 12.9 bits; conditional E-value: 2.1e-05
-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC.--EEEEEEEE CS
Methyltr_RsmB-F 10 etilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg.elfDriLlDa 87
++l l aa G +h+++++ ++g v+Av+ +++ + +++ k+ tn+i + Dar+++k+ ++D i D+
sp|Q9W1V3|FBRL_DROME 186 SKVLYLGAASGTTVSHVSDVVGPEGLVYAVEFSHRSGRDLINVAKK--RTNIIPIIEDARHPHKYRMLvGMVDTIFADV 262
6789999************************987665555544433..3456778889999998754435888888876 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (200 nodes)
Target sequences: 13786 (7394964 residues searched)
Passed MSV filter: 548 (0.0397505); expected 275.7 (0.02)
Passed bias filter: 344 (0.0249529); expected 275.7 (0.02)
Passed Vit filter: 43 (0.00311911); expected 13.8 (0.001)
Passed Fwd filter: 6 (0.000435224); expected 0.1 (1e-05)
Initial search space (Z): 13786 [actual number of targets]
Domain search space (domZ): 6 [number of targets reported over threshold]
# CPU time: 0.13u 0.01s 00:00:00.14 Elapsed: 00:00:00.06
# Mc/sec: 24649.88
//
[ok]