# 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]