# 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/Ciona_intestinalis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.5e-25 88.1 0.0 9.8e-25 87.5 0.0 1.2 1 tr|F6X2Y1|F6X2Y1_CIOIN tRNA wybutosine-synthesizing protein Domain annotation for each sequence (and alignments): >> tr|F6X2Y1|F6X2Y1_CIOIN tRNA wybutosine-synthesizing protein 4 isoform X2 OS=Ciona intestinalis OX=7719 GN=LOC1001819 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 87.5 0.0 5.9e-29 9.8e-25 24 187 .. 15 200 .. 5 201 .. 0.86 Alignments for each domain: == domain 1 score: 87.5 bits; conditional E-value: 5.9e-29 CCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT...ESEEEEET-TT--HHHHHHHTT CS LCM 24 laeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saag...irqvViLGAGlDtrayRlpwla 106 + +s+vs+ ++++++ + d+ l+++ ++k ++++ +nrg+++R ++++++l+ +l s+ ++qvV LGAG+Dtr + + a tr|F6X2Y1|F6X2Y1_CIOIN 15 NDSSIVSKLSMATCNYIDDQYLKYF--VKKPQRRSALINRGYYLRCKAIEHVLELFLsSN-LecnVKQVVSLGAGFDTRYFWARKYA 98 5566777777776666666666665..6788899***********************543.2479***************9888888 PP .....TEEEEEEE-HHHHHHHHHHHHHSHHHHH..............EECSEEE--TTS..HHHHHHHHHHTT-TTSEEEEEEESSG CS LCM 107 .....glkvfEvDqpevlefKretlkelgaepk..............aelryvpvDlrd..ddwpeaLaaagfdpekptawvaeGvl 172 +l+++E+D+pev+++K+++ e g ++ ++ ++ Dlr+ ++ +aL + g + ++pt++++e+vl tr|F6X2Y1|F6X2Y1_CIOIN 99 qklkkELVYVEIDFPEVMKRKISLGTEFGLYQSdfsneqsqniilrtENYAAIGHDLRNmgPQLHSALSEVGVNFSQPTLFISEVVL 185 88999******************99999998889*****9999997555666******988899999999***************** PP GGS-HHHHHHHHHHH CS LCM 173 vYLtpeasdalleai 187 +YLt +s++++++i tr|F6X2Y1|F6X2Y1_CIOIN 186 AYLTGSQSSKVISWI 200 *************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 16678 (5370518 residues searched) Passed MSV filter: 427 (0.0256026); expected 333.6 (0.02) Passed bias filter: 307 (0.0184075); expected 333.6 (0.02) Passed Vit filter: 25 (0.00149898); expected 16.7 (0.001) Passed Fwd filter: 1 (5.99592e-05); expected 0.2 (1e-05) Initial search space (Z): 16678 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.09u 0.01s 00:00:00.10 Elapsed: 00:00:00.09 # Mc/sec: 11218.42 // [ok]