# 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/LCM.hmm.txt # target sequence database: proteomes/Trypanosoma_brucei.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: LCM [M=188] Accession: PF04072.14 Description: Leucine carboxyl methyltransferase 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-26 92.0 0.2 6.7e-26 90.3 0.0 1.9 2 tr|Q38BF5|Q38BF5_TRYB2 Leucine carboxyl methyltransferase 1 3.1e-22 78.4 0.0 6.8e-19 67.5 0.0 2.5 2 tr|Q38DD6|Q38DD6_TRYB2 Uncharacterized protein OS=Trypanosom Domain annotation for each sequence (and alignments): >> tr|Q38BF5|Q38BF5_TRYB2 Leucine carboxyl methyltransferase 1 OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.6 0.1 1 4.3e+03 88 94 .. 5 11 .. 3 12 .. 0.87 2 ! 90.3 0.0 1.6e-29 6.7e-26 2 185 .. 30 204 .. 29 206 .. 0.84 Alignments for each domain: == domain 1 score: -3.6 bits; conditional E-value: 1 EEEET-T CS LCM 88 vViLGAG 94 vV LG+G tr|Q38BF5|Q38BF5_TRYB2 5 VVTLGCG 11 9*****9 PP == domain 2 score: 90.3 bits; conditional E-value: 1.6e-29 HHHHHHHHHHHHHH..HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT.. CS LCM 2 lgvaaaralesrra..adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag.. 84 l+++a a++ r ++ +++Dp+a+ +v ++++ ++p +nrg+++Rt +++ l ++++ ++ tr|Q38BF5|Q38BF5_TRYB2 30 LQQTAYSACSRRVHcvSKTYLNDPFASCFVMDSTVM---------------------NSPLMNRGTWLRTEAIERSLLHFARGQCge 95 566666666666655678**************9664.....................489*********************544468 PP .ESEEEEET-TT--HHHHHHHTT....TEEEEEEE-HHHHHHHHHHHHHSHHHHH...EECSEEE--TTS.HHHHHHHHHHTT-TTS CS LCM 85 .irqvViLGAGlDtrayRlpwla....glkvfEvDqpevlefKretlkelgaepk...aelryvpvDlrd.ddwpeaLaaagfdpek 162 qv+ G+G+Dt+ +Rl++ + k++E+D+p+++++Kr+++++ +a+++ e ++v++Dlr+ d+++e L+++ ++ ++ tr|Q38BF5|Q38BF5_TRYB2 96 vGLQVISFGSGMDTLYFRLKKDQpdihIAKYIELDFPDLVARKRRVIDDTEALTRytgPEYELVACDLRKtDEMRELLKKSALN-NA 181 745****************99998898899**********************99953335556*******67777777777665.** PP EEEEEEESSGGGS-HHHHHHHHH CS LCM 163 ptawvaeGvlvYLtpeasdalle 185 pt+++ae + vYL++ s+ ll tr|Q38BF5|Q38BF5_TRYB2 182 PTVILAEMIFVYLEERVSSTLLT 204 **************999988875 PP >> tr|Q38DD6|Q38DD6_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 8.1 0.0 0.00026 1.1 1 80 [. 61 120 .. 61 133 .. 0.77 2 ! 67.5 0.0 1.6e-22 6.8e-19 87 188 .] 149 274 .. 141 274 .. 0.92 Alignments for each domain: == domain 1 score: 8.1 bits; conditional E-value: 0.00026 HHHHHHHHHHHHHHH..HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH CS LCM 1 algvaaaralesrra..adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal 80 a++ ++ + + s+r+ ++g++dDp+ + +v++ ++ ++p +nrg+++R ++ d++++++ tr|Q38DD6|Q38DD6_TRYB2 61 AVQRTNDDSVVSKRSavHKGYFDDPYLRFFVKKLSR----------------------RSPLINRGYYLRMLVITDVIERCI 120 577777887777774489999999999999999855......................478888888888888888888877 PP == domain 2 score: 67.5 bits; conditional E-value: 1.6e-22 EEEEET-TT--HHHHHHHTT...TEEEEEEE-HHHHHHHHHHHHHSHHHHH.................EECSEEE--TTS..HHHHH CS LCM 87 qvViLGAGlDtrayRlpwla...glkvfEvDqpevlefKretlkelgaepk.................aelryvpvDlrd..ddwpe 151 qv+ LGAG+Dt+a+Rl++ +++++EvD+p+v++ K ++k ++ + ++ v+ Dlr+ d+ + tr|Q38DD6|Q38DD6_TRYB2 149 QVISLGAGYDTLAMRLKQRPdygNVHFYEVDFPAVMQSKSMLVKMAPFGSFpedivadpggelvklygNNYWAVGTDLRStnRDLVT 235 9****************998899*********************999887789********999887666667******88889999 PP HHHHH..TT-TTSEEEEEEESSGGGS-HHHHHHHHHHHH CS LCM 152 aLaaa..gfdpekptawvaeGvlvYLtpeasdalleaia 188 L++ +f+++ pt++ ae+v+ Y++p a+++l+++ia tr|Q38DD6|Q38DD6_TRYB2 236 CLREVspQFSTDNPTVLYAECVMQYMPPVAASHLIKQIA 274 9999999*****************************996 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 289 (0.0337577); expected 171.2 (0.02) Passed bias filter: 187 (0.0218432); expected 171.2 (0.02) Passed Vit filter: 13 (0.00151851); expected 8.6 (0.001) Passed Fwd filter: 2 (0.000233618); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.07u 0.00s 00:00:00.07 Elapsed: 00:00:00.04 # Mc/sec: 20362.11 // [ok]