# 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/Mus_musculus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.9e-71 239.1 1.0 7.9e-71 238.9 0.1 1.5 2 sp|Q8CE96|TRM6_MOUSE tRNA (adenine(58)-N(1))-methyltransfera Domain annotation for each sequence (and alignments): >> sp|Q8CE96|TRM6_MOUSE tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Mus musculus OX=10090 G # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 238.9 0.1 3.6e-75 7.9e-71 2 266 .. 21 255 .. 20 304 .. 0.82 2 ? -1.8 0.3 0.073 1.6e+03 90 104 .. 306 350 .. 272 355 .. 0.53 Alignments for each domain: == domain 1 score: 238.9 bits; conditional E-value: 3.6e-75 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeeeee 90 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+ f++ + + + + ++l+e ++e++e + sp|Q8CE96|TRM6_MOUSE 21 IHDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNAIGHSYGSAFDVSSGGSLQLR-------------KKLEEPASETKE---A 91 8899******.777*******************.*******************9954332.............233333333222...3 PP Gcd10p 91 ardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyllekkda 179 ++dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +t+l+++++ l ++++ +++ sp|Q8CE96|TRM6_MOUSE 92 GTDNRNIVDDGKSQKLTQDDIKALKDKGIK-GEEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIVTILKPSTRILSIMYYA-REP 178 89****************************.******************************************************.9** PP Gcd10p 180 qkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepehplkkh 266 +ki ++r +tl+++l+l+n++ag++++v+++++GL+++a++eRmg+ g+i+++++++ p ++ +f++ p+++l+ sp|Q8CE96|TRM6_MOUSE 179 GKINHMRYDTLAQMLTLGNIRAGNKMIVMETCSGLVLGAMMERMGG------FGSIIQLYPGDGPVRAATACFGF----PKSFLSGL 255 **********************************************......***********************....55555433 PP == domain 2 score: -1.8 bits; conditional E-value: 0.073 Gcd10p 90 e.............................ardnreiiddga.rQ 104 + +r+ r+ i +++ rQ sp|Q8CE96|TRM6_MOUSE 306 NivdaaennsgeqrpmeivpgdpenkepkeKRSKRDYIQEKQrRQ 350 111222223333345556666666666666777777777555355 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 22289 (11757923 residues searched) Passed MSV filter: 1232 (0.0552739); expected 445.8 (0.02) Passed bias filter: 459 (0.0205931); expected 445.8 (0.02) Passed Vit filter: 38 (0.00170488); expected 22.3 (0.001) Passed Fwd filter: 1 (4.48652e-05); expected 0.2 (1e-05) Initial search space (Z): 22289 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.28u 0.01s 00:00:00.29 Elapsed: 00:00:00.10 # Mc/sec: 35273.77 // [ok]