# 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/MT-A70.hmm.txt # target sequence database: proteomes/Caenorhabditis_elegans.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MT-A70 [M=172] Accession: PF05063.13 Description: MT-A70 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.1e-49 167.9 0.1 4.4e-49 166.9 0.1 1.5 1 sp|Q09956|DAMT1_CAEEL DNA N6-methyl methyltransferase OS=Cae Domain annotation for each sequence (and alignments): >> sp|Q09956|DAMT1_CAEEL DNA N6-methyl methyltransferase OS=Caenorhabditis elegans OX=6239 GN=damt-1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 166.9 0.1 2.2e-53 4.4e-49 1 172 [] 185 349 .. 185 349 .. 0.96 Alignments for each domain: == domain 1 score: 166.9 bits; conditional E-value: 2.2e-53 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHH.HHHHHHHCT-EEEEEEEEEEE-TTSCBS- CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellel.krellkkWglelveeltwlkinklgepia 87 fd+ii+dPpw++k+v+rk ++++++ee l++L+ip++l++++l+++W+tnr ++e+ + e+++kWg+e+v+++++lki+++g+p++ sp|Q09956|DAMT1_CAEEL 185 FDLIIADPPWFSKSVKRK---RTYQMDEEVLDCLDIPVILTHDALIAFWITNRIGIEEeMIERFDKWGMEVVATWKLLKITTQGDPVY 269 89***************9...699********************************999***************************** PP BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTSSBSEEEE..E-CTTTSS-THHHHHHHHHC.TCS-EEEET--G--TTEEEEETT CS MT-A70 88 tlesgavlkkskEhllvgvkgnvkrstdkiinknvdtdliis..perehsrKPpelfellehll.esarkLelFaRslrpgwltvGne 172 + + ++++k ++E+l+++ k++++r+++ +++++++ p+++hs+KPp l++ll+h+ e+++ LelFaRsl p++++vG+e sp|Q09956|DAMT1_CAEEL 270 DFD-NQKHKVPFESLMLAKKKDSMRKFE------LPENFVFAsvPMSVHSHKPP-LLDLLRHFGiEFTEPLELFARSLLPSTHSVGYE 349 855.679*****************9993......5999999999**********.99*******9*********************86 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 20008 (8229220 residues searched) Passed MSV filter: 594 (0.0296881); expected 400.2 (0.02) Passed bias filter: 397 (0.0198421); expected 400.2 (0.02) Passed Vit filter: 20 (0.0009996); expected 20.0 (0.001) Passed Fwd filter: 1 (4.998e-05); expected 0.2 (1e-05) Initial search space (Z): 20008 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.07 # Mc/sec: 20220.37 // [ok]