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