# 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/Trichoplax_adhaerens.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.1e-34 117.8 0.0 4.2e-23 81.6 0.0 2.1 2 tr|B3RL43|B3RL43_TRIAD Uncharacterized protein OS=Trichoplax 9.1e-09 34.5 0.2 0.00057 18.7 0.0 2.2 2 tr|B3SFP3|B3SFP3_TRIAD Uncharacterized protein (Fragment) OS Domain annotation for each sequence (and alignments): >> tr|B3RL43|B3RL43_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_49657 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 81.6 0.0 7.3e-27 4.2e-23 1 123 [. 9 138 .. 9 146 .. 0.92 2 ! 34.5 0.0 1.5e-12 8.7e-09 285 332 .. 202 251 .. 156 254 .] 0.87 Alignments for each domain: == domain 1 score: 81.6 bits; conditional E-value: 7.3e-27 -EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--..-SEEEE---GGGTECTSEE-- CS DNA_methylase 1 lkfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.diDlLvgGfPCqpfSvaGkqkg 82 ++++++ +GiGG++++l++ + ++ aa++i+ a ++Y+ nf+ v + +i ki+ k+l+e + Dl + ++PCqpf Gk+ + tr|B3RL43|B3RL43_TRIAD 9 MQVVEFYSGIGGMHYALQESNinAKILAAIDINTVANNVYRHNFGntPVWQREIGKISLKELQElNGDLYTMSPPCQPFTRLGKKAD 95 789***************9998899********************8877889999999999998889******************** PP TTSGGGHHHHHHHHHHHHH-..-SEEEEEEEGGGGGGCCSHHH CS DNA_methylase 83 fedtrgtLffeiiriikekk..PkvfllENVkgLlkhdkgktl 123 +d r + f +i ++ + + Pk++llENVkg ++ ++++ tr|B3RL43|B3RL43_TRIAD 96 VNDARTSSFLHVIDLLIKMEnpPKYILLENVKGFETSAARNEF 138 **************99766666************999999887 PP == domain 2 score: 34.5 bits; conditional E-value: 1.5e-12 G--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHHHHH CS DNA_methylase 285 riRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieaiake 332 + R +tprE+a + +fP f+f+++ + + y+++Gns++v v+ +++k tr|B3RL43|B3RL43_TRIAD 202 KLRYFTPREVANIHCFPVHFNFPENATkkQCYRLLGNSLNVHVASELLKL 251 5799*********************9999****************99885 PP >> tr|B3SFP3|B3SFP3_TRIAD Uncharacterized protein (Fragment) OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_63030 PE=4 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 18.7 0.0 9.9e-08 0.00057 2 39 .. 5 44 .. 4 54 .. 0.93 2 ! 13.6 0.1 3.4e-06 0.02 64 80 .. 121 137 .. 111 140 .] 0.88 Alignments for each domain: == domain 1 score: 18.7 bits; conditional E-value: 9.9e-08 EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHH CS DNA_methylase 2 kfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYe 39 k++dLfaG GG+ g+ + + fe+++++e+d +a t++ tr|B3SFP3|B3SFP3_TRIAD 5 KIVDLFAGPGGLGEGFLSLKdaFEICVSAEMDTHARSTLR 44 79*****************99**************99876 PP == domain 2 score: 13.6 bits; conditional E-value: 3.4e-06 EEEE---GGGTECTSEE CS DNA_methylase 64 lLvgGfPCqpfSvaGkq 80 +L+gG+PCq++S G tr|B3SFP3|B3SFP3_TRIAD 121 VLIGGPPCQAYSLVGRA 137 79************965 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (335 nodes) Target sequences: 11518 (5222075 residues searched) Passed MSV filter: 1198 (0.104011); expected 230.4 (0.02) Passed bias filter: 328 (0.0284772); expected 230.4 (0.02) Passed Vit filter: 32 (0.00277826); expected 11.5 (0.001) Passed Fwd filter: 3 (0.000260462); expected 0.1 (1e-05) Initial search space (Z): 11518 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.18u 0.01s 00:00:00.19 Elapsed: 00:00:00.06 # Mc/sec: 29156.59 // [ok]