# 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/GCD14.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: GCD14 [M=247] Accession: PF08704.9 Description: tRNA methyltransferase complex GCD14 subunit 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4e-107 357.1 0.0 5.1e-107 356.8 0.0 1.1 1 tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransfe 9.4e-05 21.8 0.0 0.00014 21.3 0.0 1.2 1 sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillar Domain annotation for each sequence (and alignments): >> tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransferase catalytic subunit TRMT61A OS=Drosophila melanogaste # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 356.8 0.0 7.4e-111 5.1e-107 1 246 [. 72 309 .. 72 310 .. 0.99 Alignments for each domain: == domain 1 score: 356.8 bits; conditional E-value: 7.4e-111 GCD14 1 gfvyvlaptPeLWtlaLphrtqilytaDisliimkLelkpGsvvvesGtGsgslshaiartvaPtGhlftfefheqraekareefre 87 g+++vl+ptPeLWt++Lphrtqi+yt+Dis+ii++Le++pG+vvvesGtGsgslsh+++r+++PtGhl+tf+fhe+ra++ar+efr+ tr|Q9VBH7|Q9VBH7_DROME 72 GWAHVLQPTPELWTQTLPHRTQIIYTPDISMIIHQLEVRPGAVVVESGTGSGSLSHYFLRALKPTGHLHTFDFHEARADQARDEFRR 158 8************************************************************************************** PP GCD14 88 helsqlvtvthrDvckeGfdeevsgkvDavfLDlPaPweaiphaakalkkeggrlCsfsPCieqvqrtceaLrelgfteietlevll 174 h+++++vtv+hrDvc+ Gf +e++gk+DavfLDlPaP+ a+pha+kalk +ggr+CsfsPCieq qr++++L++lgf+ei +levl+ tr|Q9VBH7|Q9VBH7_DROME 159 HGIADFVTVYHRDVCNLGFGAELDGKADAVFLDLPAPDLAVPHAFKALKLSGGRFCSFSPCIEQSQRCIQELTKLGFNEIVSLEVLQ 245 *************************************************************************************** PP GCD14 175 rtyevrkskleeidl.aekrleenerkriekakrrkaskdakvedkssksstevllkPvselkghtgyLtfat 246 +++ v++++l++idl + k +++e++++++a ++ka k+ k++ +ss+ P ++l ghtg+Ltfat tr|Q9VBH7|Q9VBH7_DROME 246 QENVVKTRTLPVIDLeFLKLPKTDETTAATNALESKAPKEVKKYLTSSN--------P-QTLPGHTGFLTFAT 309 ***************999*******************************........7.99***********8 PP >> sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillarin OS=Drosophila melanogaster OX=7227 GN=Fib PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 21.3 0.0 2e-08 0.00014 35 150 .. 179 295 .. 173 313 .. 0.68 Alignments for each domain: == domain 1 score: 21.3 bits; conditional E-value: 2e-08 GCD14 35 kLelkpGsvvvesGtGsgslshaiartvaPtGhlftfefheqraekareefrehelsqlvtvthrDvc.keGfdeevsgkvDavfLDlP 122 ++++ pGs+v+ G+ sg+ + ++ v+P G ++++ef ++ + + + +++ +++ + D+ + + + + g vD++f D+ sp|Q9W1V3|FBRL_DROME 179 QIHMPPGSKVLYLGAASGTTVSHVSDVVGPEGLVYAVEFSHRSGRD-LINVAKKRT-NIIPI-IEDARhPHKY-RMLVGMVDTIFADVA 263 67889**********************************7544433.333333332.33333.4455524555.5677*********** PP GCD14 123 aPwea..ip.haakalkkeggrlCsf.sPCie 150 +P++ + +a + lk+ g + s + Ci+ sp|Q9W1V3|FBRL_DROME 264 QPDQGriVAlNAQHFLKNGGHFVISIkASCID 295 **985223313333344444444443245665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (247 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 1089 (0.0789932); expected 275.7 (0.02) Passed bias filter: 610 (0.0442478); expected 275.7 (0.02) Passed Vit filter: 50 (0.00362687); expected 13.8 (0.001) Passed Fwd filter: 8 (0.000580299); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.25u 0.01s 00:00:00.26 Elapsed: 00:00:00.08 # Mc/sec: 22831.95 // [ok]