# 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/Nematostella_vectensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.5e-52 179.5 0.0 2.9e-52 178.5 0.0 1.5 1 tr|A7SNZ9|A7SNZ9_NEMVE Cytosine-specific methyltransferase ( 1.6e-25 90.6 0.0 1.3e-19 71.2 0.0 2.1 2 tr|A7SQ72|A7SQ72_NEMVE Predicted protein (Fragment) OS=Nemat 1.1e-20 74.7 0.0 2.9e-20 73.3 0.0 1.6 1 tr|A7SUR9|A7SUR9_NEMVE Predicted protein (Fragment) OS=Nemat ------ inclusion threshold ------ 0.016 15.0 0.0 0.023 14.5 0.0 1.1 1 tr|A7RKB2|A7RKB2_NEMVE Predicted protein OS=Nematostella vec 0.44 10.3 6.2 1.3 8.7 6.2 1.6 1 tr|A7SQT6|A7SQT6_NEMVE Predicted protein OS=Nematostella vec Domain annotation for each sequence (and alignments): >> tr|A7SNZ9|A7SNZ9_NEMVE Cytosine-specific methyltransferase (Fragment) OS=Nematostella vectensis OX=45351 GN=v1g12549 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 178.5 0.0 5.8e-56 2.9e-52 2 335 .] 790 1243 .. 789 1243 .. 0.73 Alignments for each domain: == domain 1 score: 178.5 bits; conditional E-value: 5.8e-56 EEEEES-TT-HHHHHHHHTT.-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-............GS--...-SEEEE-- CS DNA_methylase 2 kfidLfaGiGGfrlgleaag.fecvaaseidkkaaktYeanfk..kvsigDitkidkk............klpe..diDlLvgGf 69 + +d faG+GG++ gl++ag e +a+e+++ aa+ Y+ n++ v + D +++ k k+p+ +++lL+gG+ tr|A7SNZ9|A7SNZ9_NEMVE 790 RSLDVFAGCGGLSEGLHQAGvAESLWAIEKEEPAAQAYRLNNPgcTVFTDDCNTLLKLamegeatnstgqKIPQrgEVELLCGGP 874 679*****************99**********************85555555554443555677777777***9766******** PP -GGGTECTSE..E--TTSGGGHHHHHHHHHHHHH--SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB- CS DNA_methylase 70 PCqpfSvaGk..qkgfedtrgtLffeiiriikekkPkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQn 152 PCq+fS + +++L + + + + ++P +f+lENV++ ++ +k+ +lk l l +GY++++ +l+a yGvpQ tr|A7SNZ9|A7SNZ9_NEMVE 875 PCQGFSGMNRfnTREYSLFKNSLVVSYLSYCDFYRPRFFILENVRNFVSFKKSMVLKLTLRCLLRMGYQCTFGVLQAGCYGVPQT 959 ******8775334445555679*************************************************************** PP -EEEEEEEEEGGG---.................-------SS---GGG..............TSGGGTTS-EEBTTTB--GGGSS CS DNA_methylase 153 ReRvfivgirkdleln.................vlvkvkefrfpkkkd.........leasktlrdlleekseldeektalsdkl 211 R+R +i++ ++el v v e++f+++ ++++++ dl+e + + ++ tr|A7SNZ9|A7SNZ9_NEMVE 960 RRRAIIMAAAPGEELPlypepthcfsprtcqltV--MVDEKKFTSNITrtasapyrtITVRDAMSDLPEI-----RNGASAAESS 1037 ************9998999999998888888553..3333333333334467744443333344444443.....2333333333 PP STT-EE--.......S---HHHTS--.................................B--TTS...-B----S-STT-..... CS DNA_methylase 212 vekeekikkttkaseegypkelsrnr.................................idkveelkglsatrrksgrge..... 258 +e e+ ++ +++ +++y++ l r+ + ++ ++k+l+ t++++++ + tr|A7SNZ9|A7SNZ9_NEMVE 1038 YEGEAISHFQRQIRGNQYQP-LLRDHickemsalvearmrhiplapgsdwrdlpnkeirLSDGTYSKKLQYTHHDKKNEKgssqr 1121 44444444444444455553.2222233444555577777778888777777777777734444444444444444443355667 PP .................................................-B-TTS---EEEETTEEE--T.TTGGG--B-BHHHH CS DNA_methylase 259 .................................................klekelaeagklqknkerfkpviekkriRrLtprEa 294 ++ +++ ++++k+++ v++++++R ++rE+ tr|A7SNZ9|A7SNZ9_NEMVE 1122 lrgvcscaegrpcepadrqfntlipwclphtgnrhnhwaglygrlewdgYFSTTITNPEPMGKQGR----VLHPEQHRVVSVREC 1202 777777777777777777777777777777777777777777766666555556666666666666....9************** PP HHHHT--TT-----SBHHHHHHHHTS--HHHHHHHHHHHHH CS DNA_methylase 295 arlqGfPddfefgskksaaykqiGnsVpvpvieaiakeikk 335 ar qGfPd++ f ++ ++ +q+Gn+Vp+p++ ai++eikk tr|A7SNZ9|A7SNZ9_NEMVE 1203 ARSQGFPDTYRFYGSILDKHRQVGNAVPPPLAAAIGREIKK 1243 **************************************986 PP >> tr|A7SQ72|A7SQ72_NEMVE Predicted protein (Fragment) OS=Nematostella vectensis OX=45351 GN=v1g127267 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 71.2 0.0 2.6e-23 1.3e-19 2 124 .. 338 474 .. 337 480 .. 0.91 2 ! 17.4 0.0 6.1e-07 0.003 254 332 .. 530 615 .. 514 618 .. 0.78 Alignments for each domain: == domain 1 score: 71.2 bits; conditional E-value: 2.6e-23 EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--...-SEEEE---GGGTECTSE.E- CS DNA_methylase 2 kfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.d.iDlLvgGfPCqpfSvaGk.qk 81 +++ Lf+Gi +l++ g + ++++eid++a+++ + n++ +++gDi+++++++++e +Dl++gG PCq+ S+a ++ tr|A7SQ72|A7SQ72_NEMVE 338 RVLALFDGIATGLQALNELGivSDKYYSAEIDEQAIQVTKVNHGdrITHLGDIKDLTESQIRElGpFDLVIGGSPCQDLSIANPaRR 424 7889****************99999*******************9999***********9987646****************8549* PP -TTSGGGHHHHHHHHHHHHH-........-SEEEEEEEGGGGGGCCSHHHH CS DNA_methylase 82 gfedtrgtLffeiiriikekk........PkvfllENVkgLlkhdkgktlk 124 g+ + +g Lffe+ r++ + k P + l+ENV g+ +dk kt+ tr|A7SQ72|A7SQ72_NEMVE 425 GIFEGSGRLFFEFFRLLMHAKpsrtcpsrPFFWLFENVVGMRAEDK-KTIS 474 ********************9888998889999*********9988.4444 PP == domain 2 score: 17.4 bits; conditional E-value: 6.1e-07 -STT--B-TTS---EEEETTEEE--T......TTGGG--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHHHHH CS DNA_methylase 254 sgrgeklekelaeagklqknkerfkp.....viekkriRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieaiake 332 + r++k+++ +++a++l+++k++ p + ++r+ L E+ rl GfP+ ++ + ++ +++Gn+ +vpv+++++ tr|A7SQ72|A7SQ72_NEMVE 530 KARFTKVQTITTNANSLTQTKKNILPvavndDGGQEREDILWCTEMERLFGFPSHYTDVNNMGrtQRQRLLGNAWSVPVVRHLLSP 615 4445566777777777777777666677766777888888999*************977665556*****************9865 PP >> tr|A7SUR9|A7SUR9_NEMVE Predicted protein (Fragment) OS=Nematostella vectensis OX=45351 GN=v1g132811 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 73.3 0.0 6e-24 2.9e-20 1 175 [. 7 184 .. 7 219 .. 0.82 Alignments for each domain: == domain 1 score: 73.3 bits; conditional E-value: 6e-24 -EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-.-SEEEES-TTTS-.....GS--.-SEEEE---GGGTECTSEE CS DNA_methylase 1 lkfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfkkvsigDitkidkk.....klpediDlLvgGfPCqpfSvaGkq 80 ++++++ +GiGG++++l+ ++ e+vaa ei +a ++Y nf+ ++i + + + ++p+ +v ++PCqp+ G q tr|A7SUR9|A7SUR9_NEMVE 7 FRVVEFYSGIGGMHYALKGCKknAEVVAALEISTTANTVYGHNFPTTKIWNCNIEVCElcnvtTMPA--IYMVMSPPCQPYTWVGLQ 91 799***************99999*********************95555554444444454458885..67899************* PP --TTSGGGHHHHHHHHHHHHH-..-SEEEEEEEGGGGGGCCS.HHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGG CS DNA_methylase 81 kgfedtrgtLffeiiriikekk..PkvfllENVkgLlkhdkg.ktlkvilekleelGYevsykllnasdyGvpQnReRvfivgirkd 164 + +d+r f i ++k + Pk+ l+ENVkg ++d++ l +++++ ++ + ++G+p +R R ++++ r++ tr|A7SUR9|A7SUR9_NEMVE 92 GASKDPRALSFLHILSLLKRLQhpPKYWLIENVKGFETSDTRfYILLAFCNSF----------IVSSPQFGIPNSRLRYYLLAKRHP 168 ********99**********9999***************99766777777766..........89*********************9 PP G---.....------- CS DNA_methylase 165 leln.....vlvkvke 175 l+ + + + +++ tr|A7SUR9|A7SUR9_NEMVE 169 LTFStamgnKFSYYEQ 184 9866444443333333 PP >> tr|A7RKB2|A7RKB2_NEMVE Predicted protein OS=Nematostella vectensis OX=45351 GN=v1g84090 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.5 0.0 4.7e-06 0.023 3 44 .. 151 201 .. 150 221 .. 0.74 Alignments for each domain: == domain 1 score: 14.5 bits; conditional E-value: 4.7e-06 EEEES-TT-HHHHHHH.HTT-EEEEEEE--TTTHHHHHHH-.........- CS DNA_methylase 3 fidLfaGiGGfrlgle.aagfecvaaseidkkaaktYeanfk........k 44 ++dLfaGiG f l + a+ v a e++++a++ +++n + tr|A7RKB2|A7RKB2_NEMVE 151 IVDLFAGIGYFVLPFLvHANAAFVHACEWNPHAVEALRRNLVqnkvnekcE 201 69***********8751666********************97543222222 PP >> tr|A7SQT6|A7SQT6_NEMVE Predicted protein OS=Nematostella vectensis OX=45351 GN=v1g215965 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.7 6.2 0.00028 1.3 185 278 .. 14 122 .. 2 144 .. 0.56 Alignments for each domain: == domain 1 score: 8.7 bits; conditional E-value: 0.00028 ....TSGGGTTS-EE....BTTTB--GGGSSSTT-EE--.......S---HHHTS--................B--TTS...-B--- CS DNA_methylase 185 easktlrdlleekse....ldeektalsdklvekeekikkttkaseegypkelsrnr................idkveelkglsatr 251 ++++ ++++e+s+ +++kt + +k++ k +++ ++ k++++gy ++++ ++ +++++++++++ tr|A7SQT6|A7SQT6_NEMVE 14 NKTGYHPKIMTEASDsvrqPKHNKTGYHTKIMTKASTSVRQPKHNKTGYHPKIMTGAirqpkhnktgyhpkimTEASDSVRQPKHNK 100 455555566666666666789999999999999999999999999999999666666444444444444444433333333333333 PP -S-STT--B-TTS---EEEETTEEE-- CS DNA_methylase 252 rksgrgeklekelaeagklqknkerfk 278 + ++ k +++a++++++ +++k tr|A7SQT6|A7SQT6_NEMVE 101 TGYHT-----KIMTKASTSVRQPKHNK 122 33332.....33333333333333333 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (335 nodes) Target sequences: 24428 (8257945 residues searched) Passed MSV filter: 1625 (0.066522); expected 488.6 (0.02) Passed bias filter: 723 (0.0295972); expected 488.6 (0.02) Passed Vit filter: 65 (0.00266088); expected 24.4 (0.001) Passed Fwd filter: 6 (0.00024562); expected 0.2 (1e-05) Initial search space (Z): 24428 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.27u 0.01s 00:00:00.28 Elapsed: 00:00:00.09 # Mc/sec: 30737.91 // [ok]