# 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/Drosophila_melanogaster.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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
5.9e-44 150.4 0.0 8.5e-44 149.9 0.0 1.2 1 tr|Q9VKB3|Q9VKB3_DROME Methyltransferase 2, isoform C OS=Dro
Domain annotation for each sequence (and alignments):
>> tr|Q9VKB3|Q9VKB3_DROME Methyltransferase 2, isoform C OS=Drosophila melanogaster OX=7227 GN=Mt2 PE=3 SV=2
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 149.9 0.0 6.2e-48 8.5e-44 1 331 [. 3 340 .. 3 343 .. 0.80
Alignments for each domain:
== domain 1 score: 149.9 bits; conditional E-value: 6.2e-48
-EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--..-SEEEE---GGGTECTSEE-- CS
DNA_methylase 1 lkfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.diDlLvgGfPCqpfSvaGkqkg 82
+++++Lf+GiGG++++++ a+ +++vaa + + a ++Y+ n++ v++ +i+ ++ k++ + + +L+ ++PCqp +G q++
tr|Q9VKB3|Q9VKB3_DROME 3 FRVLELFSGIGGMHYAFNYAQldGQIVAALDVNTVANAVYAHNYGsnLVKTRNIQSLSVKEVTKlQANMLLMSPPCQPHTRQGLQRD 89
799******************99*********************9888899999999999999989********************* PP
TTSGGGHHHHHHHHHHHHH-.-SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGGG--- CS
DNA_methylase 83 fedtrgtLffeiiriikekk.PkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQnReRvfivgirkdleln 168
+ed+r++ ++ + +i e + +++l+ENVkg ++ ++ ++++e+le G++ +l +++vp R R ++++ +
tr|Q9VKB3|Q9VKB3_DROME 90 TEDKRSDALTHLCGLIPECQeLEYILMENVKGFESS---QARNQFIESLERSGFHWREFILTPTQFNVPNTRYRYYCIARKGAD-FP 172
******99999999998876489**********999...568999**********************************88754.44 PP
....-------SS---GGG.....TSGGGTTS-EEBTTTB--GGGSSSTT-EE--.......S---HHHTS--B--TTS...-B--- CS
DNA_methylase 169 ....vlvkvkefrfpkkkdleasktlrdlleekseldeektalsdklvekeekikkttkaseegypkelsrnridkveelkglsatr 251
++ +++ +++ ++++++ ee+ + + + d+++ k +++ +++ ++ + + + + + e+ g t
tr|Q9VKB3|Q9VKB3_DROME 173 faggKIWEEMPGAIAQNQG---LSQIAEIVEEN---VSPDFLVPDDVLTKRVLVMD-IIHPAQSRS-MCFTKGYTHYTEGTGSAYTP 251
5788888888888888875...77899999996...55556666665555555543.333333333.33333233333333333333 PP
-S-STT--B-TTS---EEEETTEEE--T........TTGGG--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHH CS
DNA_methylase 252 rksgrgeklekelaeagklqknkerfkp.......viekkriRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieai 329
+ ++ +++ + ++e+ +++++++++++ ++++ r R +tprE+arl++fP++fef+++++ ++y+++Gns++v+v+ ++
tr|Q9VKB3|Q9VKB3_DROME 252 LSEDESHRIFELVKEIDTSNQDASKSEKilqqrldLLHQVRLRYFTPREVARLMSFPENFEFPPETTnrQKYRLLGNSINVKVVGEL 338
33333334445555555555555544446667777********************************9****************998 PP
HH CS
DNA_methylase 330 ak 331
+k
tr|Q9VKB3|Q9VKB3_DROME 339 IK 340
87 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (335 nodes)
Target sequences: 13786 (7394964 residues searched)
Passed MSV filter: 1458 (0.105759); expected 275.7 (0.02)
Passed bias filter: 480 (0.0348179); expected 275.7 (0.02)
Passed Vit filter: 51 (0.00369941); expected 13.8 (0.001)
Passed Fwd filter: 3 (0.000217612); expected 0.1 (1e-05)
Initial search space (Z): 13786 [actual number of targets]
Domain search space (domZ): 1 [number of targets reported over threshold]
# CPU time: 0.37u 0.02s 00:00:00.39 Elapsed: 00:00:00.12
# Mc/sec: 20644.27
//
[ok]