# 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/DNA_methylase.hmm.txt
# target sequence database: proteomes/Dictyostelium_discoideum.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: DNA_methylase [M=335]
Accession: PF00145.17
Description: C-5 cytosine-specific DNA methylase
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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
3.6e-42 144.4 3.6 4.3e-42 144.2 3.6 1.1 1 sp|Q54JH6|CMT1_DICDI DNA (cytosine-5)-methyltransferase OS=D
------ inclusion threshold ------
0.081 11.8 5.9 0.22 10.4 0.0 2.7 3 sp|Q54WD6|TRM5_DICDI tRNA (guanine(37)-N1)-methyltransferase
Domain annotation for each sequence (and alignments):
>> sp|Q54JH6|CMT1_DICDI DNA (cytosine-5)-methyltransferase OS=Dictyostelium discoideum OX=44689 GN=dnmA PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 144.2 3.6 6.7e-46 4.3e-42 2 333 .. 5 359 .. 4 361 .. 0.73
Alignments for each domain:
== domain 1 score: 144.2 bits; conditional E-value: 6.7e-46
EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-..-SEEEES-TTTS-GS--..-SEEEE---GGGTECTSEE--TTSG CS
DNA_methylase 2 kfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk.kvsigDitkidkkklpe.diDlLvgGfPCqpfSvaGkqkgfedt 86
+++++ +GiGG+++gl++ g f+++ + +i+ +a Y+ f+ + s++ i+ + ++l++ + + ++PCqpf G qk+ +d+
sp|Q54JH6|CMT1_DICDI 5 RVLEFYSGIGGMHYGLQESGvdFQVIQSFDINTNANLNYKYTFNeDSSQKSIESYSVEELEGfKANAWLMSPPCQPFTRLGLQKDDQDN 93
899*****************************************999**************999************************* PP
GGHHHHHHHHHHHHH-..-SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGGG---..... CS
DNA_methylase 87 rgtLffeiiriikekk..PkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQnReRvfivgirkdleln..... 168
r + ff ++ ++k P+++l+ENV g k+ +++t +++l++l +++Y+++ l +++G + +R R f+++ r+++ ln
sp|Q54JH6|CMT1_DICDI 94 RTNSFFHLLDVLTKIKdpPTYILIENVFGFAKKGSSNTRDHLLDTLIKMNYSFQEFHLSPQQFGLANQRLRYFCIAKRNGK-LNfkkeq 181
**************9999************************************999999*******************97.3334333 PP
..........-------SS---GGG.....TSGGGTTS-EE.....BTTTB--GGGSSSTT-EE--.......S---HHHTS--..... CS
DNA_methylase 169 ..........vlvkvkefrfpkkkdleasktlrdlleekse.....ldeektalsdklvekeekikkttkaseegypkelsrnr..... 237
+l++ ++++ +++k+ ++ k+l +++ + de ++++++l+++e + +
sp|Q54JH6|CMT1_DICDI 182 dkhnekvdenKLNNNSNNNNEQNKY-DNLKILDHIPGYDFHttleeCDEISNYFDKDLTDDELYE------------------Kykvph 251
3333333333222222222222222.333333333333222333333333333333333333333..................033344 PP
.............B--TTS...-B----S-STT--B-TTS---EEEETTEEE--T.......TTGGG--B-BHHHHHHHHT--TT---- CS
DNA_methylase 238 .............idkveelkglsatrrksgrgeklekelaeagklqknkerfkp......viekkriRrLtprEaarlqGfPddfefg 307
+++++ ++ +++ k + +g+ ++ +++fk+ ++ + + R ++p+E rl GfP++f+f+
sp|Q54JH6|CMT1_DICDI 252 nlllskgmlfdikQKDSKTSNCVTKSYGKFIE---------GTGSIIQMDNNFKAdindnkSLIPLKLRYFSPKEITRLHGFPEEFKFS 331
44444444443333334444444444444443.........3333333333333333333389************************** PP
-SBH..HHHHHHHTS--HHHHHHHHHHH CS
DNA_methylase 308 skks..aaykqiGnsVpvpvieaiakei 333
+k + + y++iGns++v+++ +++k +
sp|Q54JH6|CMT1_DICDI 332 PKLTtiQCYRLIGNSLNVKIVSELLKVL 359
**999******************99976 PP
>> sp|Q54WD6|TRM5_DICDI tRNA (guanine(37)-N1)-methyltransferase OS=Dictyostelium discoideum OX=44689 GN=trmt5 PE=3 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -1.6 0.0 0.14 9.2e+02 128 156 .. 90 118 .. 34 172 .. 0.62
2 ? 10.4 0.0 3.4e-05 0.22 4 55 .. 218 272 .. 216 281 .. 0.81
3 ? 1.5 2.7 0.017 1.1e+02 208 264 .. 394 446 .. 346 460 .] 0.47
Alignments for each domain:
== domain 1 score: -1.6 bits; conditional E-value: 0.14
HHHHHHTEEEEEEEEEGGGTT-SB--EEE CS
DNA_methylase 128 ekleelGYevsykllnasdyGvpQnReRv 156
+ +++ Ye +k l +d G+p + eR+
sp|Q54WD6|TRM5_DICDI 90 LNYNNFSYEQVLKTLLPKDVGIPFSFERI 118
23334444444666666666666666665 PP
== domain 2 score: 10.4 bits; conditional E-value: 3.4e-05
EEES-TT-HHHHHHHHTT-EEEEEEE--TTTHHHHHHH-....-SEEEES-TTTS CS
DNA_methylase 4 idLfaGiGGfrlgleaagfecvaaseidkkaaktYeanfk...kvsigDitkidk 55
d+faG+G f+l ++++ v+a++ +++ +k + n k +s +i+++d
sp|Q54WD6|TRM5_DICDI 218 CDMFAGVGPFALPAAKIKKCKVYANDLNPSSVKYMKENAKtnrLESKVEISNLDA 272
7*****************999****************998654444445555554 PP
== domain 3 score: 1.5 bits; conditional E-value: 0.017
GGSSSTT-EE--.......S---H..HHTS--B--TTS...-B----S-STT--B-TTS CS
DNA_methylase 208 sdklvekeekikkttkaseegypk..elsrnridkveelkglsatrrksgrgeklekel 264
++ ++++ +++++t+++s ++ + + + n+ ++ve+ + +++++++ ++ e+
sp|Q54WD6|TRM5_DICDI 394 TTPTTSNTNTSTTTSTTSTSTTTTesTNTNNSANNVEDKNNKKRNSTEDSN------ET 446
333333333333333333322222212222223333333333333333332......22 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (335 nodes)
Target sequences: 12739 (6856946 residues searched)
Passed MSV filter: 4988 (0.391553); expected 254.8 (0.02)
Passed bias filter: 793 (0.0622498); expected 254.8 (0.02)
Passed Vit filter: 527 (0.041369); expected 12.7 (0.001)
Passed Fwd filter: 5 (0.000392495); expected 0.1 (1e-05)
Initial search space (Z): 12739 [actual number of targets]
Domain search space (domZ): 2 [number of targets reported over threshold]
# CPU time: 0.87u 0.02s 00:00:00.89 Elapsed: 00:00:00.25
# Mc/sec: 9188.31
//
[ok]