# 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/Ciona_intestinalis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.9e-47 160.3 0.0 1.4e-46 159.3 0.0 1.5 1 tr|F6XUV4|F6XUV4_CIOIN Uncharacterized protein OS=Ciona inte 5e-46 157.5 0.0 6e-46 157.2 0.0 1.1 1 tr|F6YI25|F6YI25_CIOIN tRNA (Cytosine(38)-C(5))-methyltransf 0.00043 19.6 0.1 0.00063 19.1 0.1 1.4 1 tr|F6WB92|F6WB92_CIOIN methyltransferase-like protein 5 isof ------ inclusion threshold ------ 0.055 12.7 1.3 0.074 12.3 1.3 1.1 1 tr|H2XRT1|H2XRT1_CIOIN Uncharacterized protein OS=Ciona inte 6.5 5.9 18.3 7.4 5.7 18.3 1.1 1 tr|F6YRL7|F6YRL7_CIOIN Uncharacterized protein OS=Ciona inte Domain annotation for each sequence (and alignments): >> tr|F6XUV4|F6XUV4_CIOIN Uncharacterized protein OS=Ciona intestinalis OX=7719 GN=zf(cxxc)-5 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 159.3 0.0 4.1e-50 1.4e-46 2 335 .] 796 1249 .. 795 1249 .. 0.75 Alignments for each domain: == domain 1 score: 159.3 bits; conditional E-value: 4.1e-50 EEEEES-TT-HHHHHHHHTT.-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-............GS--...-SEEEE-- CS DNA_methylase 2 kfidLfaGiGGfrlgleaag.fecvaaseidkkaaktYeanfk..kvsigDitkidkk............klpe..diDlLvgGf 69 + +d f+G+GG++ g+++ag e +a+e ++ aa+ Y+ n++ v ++D + + + lp+ d++lL+gG+ tr|F6XUV4|F6XUV4_CIOIN 796 RSLDVFSGCGGLSEGFHQAGiAEPSYAIELWEPAAQAYRLNNPgaTVFTEDCNVLLEMvmngeerskcgqRLPQkgDVELLCGGP 880 679*****************99**********************855566666555445566677889999998667******** PP -GGGTECTSE..E--TTSGGGHHHHHHHHHHHHH--SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB- CS DNA_methylase 70 PCqpfSvaGk..qkgfedtrgtLffeiiriikekkPkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQn 152 PCq+fS ++ + +++L + + + + ++P +fllENV++ ++ ++ +lk l +l +GY++++ +l+a +yGv+Q tr|F6XUV4|F6XUV4_CIOIN 881 PCQGFSGMNRfnSREYSRFKNSLVVSYLSYCDYYRPRFFLLENVRNFVSFKNCMVLKLSLAALLRMGYQCTFGVLQAGHYGVSQT 965 ******998833445555556**************************************************************** PP -EEEEEEEEEGGG---...............-------SS---GGG..............TSGGGTTS-EEBTTTB--GGGSSST CS DNA_methylase 153 ReRvfivgirkdleln...............vlvkvkefrfpkkkd.........leasktlrdlleekseldeektalsdklve 213 R+R +i++ +++l ++ +v ++r++++ + ++ +++ dl++ + + + + ++ + + tr|F6XUV4|F6XUV4_CIOIN 966 RRRAIILAAAPGEQLPlypeplhtfstrggsLSAQVGDTRYTNNIKwstsapyrtITVHDAMSDLPKIPNGHAKLEMPYTGEPLT 1050 ************9987767777777666666555555555555555666665555556666666666554466666666666666 PP T-EE--.......S---HHHTS--.................................B--TTS...-B----S-STT-....... CS DNA_methylase 214 keekikkttkaseegypkelsrnr.................................idkveelkglsatrrksgrge....... 258 + +++++ k ++ +++ + ++ +k l+ ++++++g+ tr|F6XUV4|F6XUV4_CIOIN 1051 HFQRMMR-GKDYQGVVRD-----HickdmsplvaarmaliplipgsdwrdlpnkphkLSDGTTAKVLRYEYHDKKQGRsstggmr 1129 6666663.2223333333.....24444488889999999999999999999996663333333333333333333336677777 PP ...............................................-B-TTS---EEEETTEEE--T.TTGGG--B-BHHHHHH CS DNA_methylase 259 ...............................................klekelaeagklqknkerfkpviekkriRrLtprEaar 296 ++ +++ ++++k+++ v++++++R ++rE+ar tr|F6XUV4|F6XUV4_CIOIN 1130 gvcscaegracdamdrqfntlipwclphtgnrhnnwaglygrlcwdgFFSTTVTNPEPMGKQGR----VLHPQQHRVVSVRECAR 1210 7777778888887777777666666666666666666666665555545555555555555555....9**************** PP HHT--TT-----SBHHHHHHHHTS--HHHHHHHHHHHHH CS DNA_methylase 297 lqGfPddfefgskksaaykqiGnsVpvpvieaiakeikk 335 qGfPd++ f ++ ++ + +Gn+Vp+p+ +ai+ +ikk tr|F6XUV4|F6XUV4_CIOIN 1211 SQGFPDSYRFFGSILDKHREVGNAVPPPMSKAIGLQIKK 1249 **********************************99986 PP >> tr|F6YI25|F6YI25_CIOIN tRNA (Cytosine(38)-C(5))-methyltransferase-like OS=Ciona intestinalis OX=7719 GN=LOC100183289 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 157.2 0.0 1.8e-49 6e-46 1 331 [. 5 335 .. 5 339 .. 0.79 Alignments for each domain: == domain 1 score: 157.2 bits; conditional E-value: 1.8e-49 -EEEEES-TT-HHHHHHHHTT...-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--..-SEEEE---GGGTECTSEE- CS DNA_methylase 1 lkfidLfaGiGGfrlgleaag...fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.diDlLvgGfPCqpfSvaGkqk 81 lk+++L +GiGG++++l a+ e+v++++i + a +Y+ nf+ k+ + i+ ++ k++ + + L+ ++PCqpf G qk tr|F6YI25|F6YI25_CIOIN 5 LKVLELYSGIGGMHYALLGANlknCEVVCSVDISPAASLVYKHNFPgtKHWERSIEGFSAKDFDNmGFNTLMMSPPCQPFTRVGLQK 91 589*************99999989****************************************977******************** PP -TTSGGGHHHHHHHHHHHHH-..-SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGGG- CS DNA_methylase 82 gfedtrgtLffeiiriikekk..PkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQnReRvfivgirkdle 166 +++d+r f++ ++++ + + P+++l+ENVkg ++ k+ + +l+ le+l Y+ ll +++Gvp +R R ++++ rk+l+ tr|F6YI25|F6YI25_CIOIN 92 DINDPRTRSFIYLMKVLPQLSnpPTYILMENVKGFENS---KAHDMFLQVLEQLEYSTAQFLLTPKQFGVPNSRLRYYLLAKRKPLQ 175 *****************987777***********9888...67999999999*********************************99 PP ---------SS---GGG.....TSGGGTTS-EEBTTTB--GGGSSSTT-EE--.......S---H.HHTS--B--TTS...-B---- CS DNA_methylase 167 lnvlvkvkefrfpkkkdleasktlrdlleekseldeektalsdklvekeekikkttkaseegypk.elsrnridkveelkglsatrr 252 +v+++++ ++ + + +k l++ ++ +e++++++ k++e++ + ++ t + +k l + ++ ++ + +s+t tr|F6YI25|F6YI25_CIOIN 176 FPDEVNTPNEVIETMPQ-ALQKCLNSNAG-----KEQNEKHTIKTLENYLEEAESTLKEHALTEKtLLRYLQVMDIVTPRHQSSTCF 256 87777777777777764.33444444444.....45566666666666666665555555444443444444333333344444443 PP S-STT-..-B-TTS---EEEETTEEE--T.TTGGG--B-BHHHHHHHHT--TT.-----SBH..HHHHHHHTS--HHHHHHHHH CS DNA_methylase 253 ksgrge..klekelaeagklqknkerfkpviekkriRrLtprEaarlqGfPdd.fefgskks..aaykqiGnsVpvpvieaiak 331 +++ g + + ++ ++g +k ++ +++ R +t rE+a +++fP+d f f+++ + ++yk++Gns++v v+ ++k tr|F6YI25|F6YI25_CIOIN 257 TKSYGYyaEGTGSVLNMGGEIKMSN-----PLHNNLRYFTAREVANIMCFPKDkFHFPENFTrkQKYKLLGNSLNVYVVSCLLK 335 3333221122244555555555555.....456899***************9999999887777**************988876 PP >> tr|F6WB92|F6WB92_CIOIN methyltransferase-like protein 5 isoform X4 OS=Ciona intestinalis OX=7719 GN=LOC100184601 PE= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.1 0.1 1.9e-07 0.00063 3 70 .. 52 130 .. 50 160 .. 0.76 Alignments for each domain: == domain 1 score: 19.1 bits; conditional E-value: 1.9e-07 EEEES-TT-HHHHHHHHTT-EEEEEEE--TTTHHHHHHH-.........-SEEEES-TTTS-...GS--.-SEEEE--- CS DNA_methylase 3 fidLfaGiGGfrlgleaagfecvaaseidkkaaktYeanfk........kvsigDitkidkk...klpediDlLvgGfP 70 + dL +G+G +++g + g ++++ ++id++a ++++ n + + ++ Di ++ + l++ +D ++ +P tr|F6WB92|F6WB92_CIOIN 52 VCDLGCGCGILSIGSSLLGANHCLGIDIDEDALEIFQSNCEayelnnvvECIQADIARFSPSrnmILAKRFDTVLMNPP 130 78******************************99999999866665556666666666666664323324666666555 PP >> tr|H2XRT1|H2XRT1_CIOIN Uncharacterized protein OS=Ciona intestinalis OX=7719 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.3 1.3 2.2e-05 0.074 123 255 .. 65 217 .. 61 238 .. 0.61 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 2.2e-05 HHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGGG---..........-------SS---GGG.....TSGGGTTS-EE CS DNA_methylase 123 lkvilekleelGYevsykllnasdyGvpQnReRvfivgirkdleln..........vlvkvkefrfpkkkdleasktlrdlleekse 199 lk i e+l ++ Y+ +y+ +sd + + eRv ++++++ l+ + + +k k+ e+ k+++++ ++ + tr|H2XRT1|H2XRT1_CIOIN 65 LKLIKEELIDAVYSDTYTTSSSSDGRNSTQLERVPYFALLQRLHEQrsestelmegK--------INKVKQ-ENEKIMKEIRDKDVS 142 78899*************************************999899555444442........444443.444444444444444 PP ..........BTTTB--GGGSSSTT-EE--.......S---H.........HHTS--B--TTS...-B----S-S CS DNA_methylase 200 ..........ldeektalsdklvekeekikkttkaseegypk.........elsrnridkveelkglsatrrksg 255 l++++ +l++ +++ke++ + +++ +++ + e+ ++ ++ + +l+g+s + rk++ tr|H2XRT1|H2XRT1_CIOIN 143 iaelqalvesLKQNEVELNEVILSKEAELDDINSKVKDIKVHaddktqayeEVLKELKQIIGNLRGESGNLRKYK 217 444222222244444444444444444444433333333332333333333555555666777777777777766 PP >> tr|F6YRL7|F6YRL7_CIOIN Uncharacterized protein OS=Ciona intestinalis OX=7719 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 5.7 18.3 0.0022 7.4 165 299 .. 8 149 .. 1 169 [. 0.59 Alignments for each domain: == domain 1 score: 5.7 bits; conditional E-value: 0.0022 G---.-------SS---GGG.....TSGGGTTS-EEBTTTB--GGGSSSTT-EE--.......S---HHHTS--B--TTS...-B-- CS DNA_methylase 165 leln.vlvkvkefrfpkkkdleasktlrdlleekseldeektalsdklvekeekikkttkaseegypkelsrnridkveelkglsat 250 ++++ v kvk+ +++kk+ + ++ r+++++ s++ +++++ s + +++e+k+kk++ + ++ +++ s++r++k+ + k+ +++ tr|F6YRL7|F6YRL7_CIOIN 8 KDVReVVHKVKKAKHEKKS---EDSNKRSYFSSTSKKRKKRYSSSSCSTDSERKVKKKKRKRTTSEDTSSSEDRTQKTLKRKKAKHK 91 3333455555555555554...67888999999888999999999999999999987666655555555555554444444444444 PP --S-STT-........-B-TTS---EEEETTEEE--T..TTGGG--B-BHHHHHHHHT CS DNA_methylase 251 rrksgrge........klekelaeagklqknkerfkp.viekkriRrLtprEaarlqG 299 +++++++ + e++ +++ +l knke+ + +R++t++E+ + q+ tr|F6YRL7|F6YRL7_CIOIN 92 KKTKKKKMkkektiktTDENSPKNIKSLKKNKEETADpDPRALAMRPMTQEEWEKRQS 149 3333332234444433333333333444444443332566667788888888877665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (335 nodes) Target sequences: 16678 (5370518 residues searched) Passed MSV filter: 1418 (0.0850222); expected 333.6 (0.02) Passed bias filter: 485 (0.0290802); expected 333.6 (0.02) Passed Vit filter: 44 (0.00263821); expected 16.7 (0.001) Passed Fwd filter: 5 (0.000299796); expected 0.2 (1e-05) Initial search space (Z): 16678 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.20u 0.01s 00:00:00.21 Elapsed: 00:00:00.06 # Mc/sec: 29985.39 // [ok]