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