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