# 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/Gcd10p.hmm.txt # target sequence database: proteomes/Trypanosoma_brucei.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Gcd10p [M=300] Accession: PF04189.12 Description: Gcd10p family 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.7e-49 165.2 0.0 1.1e-48 164.8 0.0 1.1 1 tr|Q38A72|Q38A72_TRYB2 tRNA (adenine(58)-N(1))-methyltransfe Domain annotation for each sequence (and alignments): >> tr|Q38A72|Q38A72_TRYB2 tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Trypanosoma brucei br # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 164.8 0.0 1.3e-52 1.1e-48 2 288 .. 31 299 .. 30 312 .. 0.81 Alignments for each domain: == domain 1 score: 164.8 bits; conditional E-value: 1.3e-52 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrv.vpaaeleaeslaeeeaeeeee 87 I+ +++v+++ + +iv+v+ +++++Lg gs++l+++ g ++g + +++ rv vp+++ ++++ e e+ e tr|Q38A72|Q38A72_TRYB2 31 ITLGDYVMIT-GGGMKRIVHVQAGGKLRLGSSGSVQLDKLAGVRFGEVVYYDPRS--------RVfVPTNDYPDLDITTLE-EHIE- 106 8999******.66699******************************954433333........330555544444444443.3333. PP Gcd10p 88 eeeardnreiiddgarQkLtkeeIeeLKke.gasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyl 173 ++rdnr+++d++++Q L++eeI+e+++e g+ + + ++l+e++++f+ KTa+SqeKYl++Kkk+Y +++++++++ ++e++ tr|Q38A72|Q38A72_TRYB2 107 --DGRDNRHLVDENKSQVLSNEEIAEMRREkGV---DVFLNTLVEKSATFHAKTAYSQEKYLRKKKKRYGVLYKIERVTPDGVAETY 188 ..49**************************333...699************************************************ PP Gcd10p 174 lekk.....daq.kilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneq.pnlsllkyfn 253 l + +++ ++l+lr++tl+l+l++++vh+g+r++v+D+t+G l aal+ R+g +g i +i++ ++ pn+ ++++ tr|Q38A72|Q38A72_TRYB2 189 LPTInptdvEPEsRVLRLRADTLALILHHSDVHSGSRVIVYDKTNGHLEAALLTRLGS------DGIIFQIMDRTAqPNMFPSQTMG 269 9863345446778*********************************************......********97766********** PP Gcd10p 254 ydaaepehplkkhlktlswlqllepeedetyeeep 288 + + + +k + +l+ eed+++ee+ tr|Q38A72|Q38A72_TRYB2 270 IE--NVRELWKAVPRNAAFLR---GEEDTENEEKA 299 95..44566665544444443...44444444333 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 388 (0.0453218); expected 171.2 (0.02) Passed bias filter: 175 (0.0204415); expected 171.2 (0.02) Passed Vit filter: 15 (0.00175213); expected 8.6 (0.001) Passed Fwd filter: 1 (0.000116809); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.10u 0.01s 00:00:00.11 Elapsed: 00:00:00.04 # Mc/sec: 32492.72 // [ok]