# 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/MafB19-deam.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MafB19-deam [M=138] Accession: PF14437.5 Description: MafB19-like deaminase 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.9e-23 83.6 0.2 3.2e-23 83.4 0.2 1.1 1 tr|A2F5I8|A2F5I8_TRIVA Cytidine and deoxycytidylate deaminas 0.0025 19.0 0.6 0.0059 17.8 0.6 1.6 1 tr|A2DS60|A2DS60_TRIVA Cytidine and deoxycytidylate deaminas Domain annotation for each sequence (and alignments): >> tr|A2F5I8|A2F5I8_TRIVA Cytidine and deoxycytidylate deaminase zinc-binding region family protein OS=Trichomonas vagi # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 83.4 0.2 1.3e-27 3.2e-23 2 106 .. 17 123 .. 16 158 .. 0.87 Alignments for each domain: == domain 1 score: 83.4 bits; conditional E-value: 1.3e-27 MafB19-deam 2 dkekafrkalaevekatkikedkigavivkdkkviargsnekelnakataHaEiaalqqAakklk..skrlddatlyvtlepCdmCa 86 +ek+++ al+ + +a++++e+ +g+vi+++++via+g n+++ +++at HaE ++ + + ++ l++ tlyvt epC mCa tr|A2F5I8|A2F5I8_TRIVA 17 MHEKYMQLALKAADDAFDAGEVAVGCVIIHNGQVIASGGNQTNAKNDATRHAELVTFKHLKESRSdyKQILKESTLYVTCEPCIMCA 103 599***************************************************88877666555114559**************** PP MafB19-deam 87 gaivqaglkklvygagnpks 106 ai g++++vyg+ n + tr|A2F5I8|A2F5I8_TRIVA 104 SAIKMMGIPRVVYGCLNDRF 123 **************997765 PP >> tr|A2DS60|A2DS60_TRIVA Cytidine and deoxycytidylate deaminase zinc-binding region family protein OS=Trichomonas vagi # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.8 0.6 2.3e-07 0.0059 46 114 .. 171 245 .. 139 261 .. 0.73 Alignments for each domain: == domain 1 score: 17.8 bits; conditional E-value: 2.3e-07 MafB19-deam 46 nakataHaEiaalqqAakklkskrlddat...lyvtlepCdmCagaivqaglkklvygagnpks...lnkLkvke 114 + +H+ i+a+ A+ ++ +t +y epC mC a+v +++ +l + npk +n+ ++++ tr|A2DS60|A2DS60_TRIVA 171 CECCFEHGVIDAVGEASRWAIQNDSYLCTnldVYCYYEPCCMCTMAMVHSRVGRLFFIEPNPKYggvMNQAHINQ 245 445678999999999988777766655554448*************************99999977777666654 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (138 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1599 (0.0318589); expected 1003.8 (0.02) Passed bias filter: 921 (0.0183503); expected 1003.8 (0.02) Passed Vit filter: 50 (0.000996214); expected 50.2 (0.001) Passed Fwd filter: 2 (3.98486e-05); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.22u 0.01s 00:00:00.23 Elapsed: 00:00:00.11 # Mc/sec: 21396.52 // [ok]