# 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/Homo_sapiens.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 9.3e-72 243.9 0.5 1.7e-71 243.1 0.4 1.5 1 sp|Q9UJA5-3|TRM6_HUMAN Isoform 3 of tRNA (adenine(58)-N(1)) 2.5e-71 242.5 0.8 4.3e-71 241.8 0.5 1.6 1 sp|Q9UJA5|TRM6_HUMAN tRNA (adenine(58)-N(1))-methyltransf 3.4e-71 242.1 1.0 3.4e-71 242.1 1.0 1.6 2 sp|Q9UJA5-2|TRM6_HUMAN Isoform 2 of tRNA (adenine(58)-N(1)) 4.7e-21 77.5 0.8 8e-21 76.7 0.2 1.7 2 tr|B4DUV6|B4DUV6_HUMAN cDNA FLJ51319, highly similar to tRN ------ inclusion threshold ------ 0.018 16.5 0.3 0.026 15.9 0.3 1.2 1 tr|H0YJ17|H0YJ17_HUMAN Neuroguidin (Fragment) OS=Homo sapie 0.2 13.0 0.1 0.39 12.1 0.1 1.4 1 tr|Q5T0V7|Q5T0V7_HUMAN Dystonin (Fragment) OS=Homo sapiens 0.41 12.0 0.3 0.51 11.7 0.3 1.2 1 tr|B7ZKW0|B7ZKW0_HUMAN ATP1B4 protein OS=Homo sapiens GN=AT 0.94 10.8 1.0 1.3 10.3 1.0 1.2 1 sp|Q9UN42|AT1B4_HUMAN Protein ATP1B4 OS=Homo sapiens GN=AT 1 10.7 1.1 1.4 10.2 1.1 1.2 1 sp|Q9UN42-2|AT1B4_HUMAN Isoform B of Protein ATP1B4 OS=Homo Domain annotation for each sequence (and alignments): >> sp|Q9UJA5-3|TRM6_HUMAN Isoform 3 of tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Homo sap # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 243.1 0.4 1.7e-75 1.7e-71 2 266 .. 20 254 .. 19 322 .. 0.83 Alignments for each domain: == domain 1 score: 243.1 bits; conditional E-value: 1.7e-75 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+ fe+ + + + k+ +++e +ae++ sp|Q9UJA5-3|TRM6_HUMAN 20 IRDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNVIGHSYGTAFEVTSGGSLQPKK---------KREEPTAETK------- 88 999*******.777*******************.********************9654433.........3333333332....... PP Gcd10p 89 eeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleylle 175 e+++dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +tv +++++ l ++++ sp|Q9UJA5-3|TRM6_HUMAN 89 EAGTDNRNIVDDGKSQKLTQDDIKALKDKGIK-GEEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIITVVKPSTRILSIMYYA 174 2389****************************.****************************************************** PP Gcd10p 176 kkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepehp 262 ++++ki ++r +tl+++l+l+n++ag++++v+++++GL+++a++eRmg+ g+i++++++ p ++ +f++ p+++ sp|Q9UJA5-3|TRM6_HUMAN 175 -REPGKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAMMERMGG------FGSIIQLYPGGGPVRAATACFGF----PKSF 250 .9************************************************......***********************....6665 PP Gcd10p 263 lkkh 266 l+ sp|Q9UJA5-3|TRM6_HUMAN 251 LSGL 254 5533 PP >> sp|Q9UJA5|TRM6_HUMAN tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Homo sapiens GN=TRMT6 P # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 241.8 0.5 4.4e-75 4.3e-71 2 266 .. 20 254 .. 19 325 .. 0.83 Alignments for each domain: == domain 1 score: 241.8 bits; conditional E-value: 4.4e-75 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeeeee 90 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+ fe+ + + + k+ +++e +ae++ e+ sp|Q9UJA5|TRM6_HUMAN 20 IRDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNVIGHSYGTAFEVTSGGSLQPKK---------KREEPTAETK-------EA 90 999*******.777*******************.********************9654433.........3333333332.......23 PP Gcd10p 91 ardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyllekkda 179 ++dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +tv +++++ l ++++ +++ sp|Q9UJA5|TRM6_HUMAN 91 GTDNRNIVDDGKSQKLTQDDIKALKDKGIK-GEEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIITVVKPSTRILSIMYYA-REP 177 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|Q9UJA5|TRM6_HUMAN 178 GKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAMMERMGG------FGSIIQLYPGGGPVRAATACFGF----PKSFLSGL 254 **********************************************......***********************....66655533 PP >> sp|Q9UJA5-2|TRM6_HUMAN Isoform 2 of tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Homo sap # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 242.1 1.0 3.4e-75 3.4e-71 2 266 .. 20 254 .. 19 322 .. 0.83 2 ? -0.9 0.7 0.35 3.4e+03 75 101 .. 293 319 .. 253 368 .. 0.53 Alignments for each domain: == domain 1 score: 242.1 bits; conditional E-value: 3.4e-75 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+ fe+ + + + k+ +++e +ae++ sp|Q9UJA5-2|TRM6_HUMAN 20 IRDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNVIGHSYGTAFEVTSGGSLQPKK---------KREEPTAETK------- 88 999*******.777*******************.********************9654433.........3333333332....... PP Gcd10p 89 eeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleylle 175 e+++dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +tv +++++ l ++++ sp|Q9UJA5-2|TRM6_HUMAN 89 EAGTDNRNIVDDGKSQKLTQDDIKALKDKGIK-GEEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIITVVKPSTRILSIMYYA 174 2389****************************.****************************************************** PP Gcd10p 176 kkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepehp 262 ++++ki ++r +tl+++l+l+n++ag++++v+++++GL+++a++eRmg+ g+i++++++ p ++ +f++ p+++ sp|Q9UJA5-2|TRM6_HUMAN 175 -REPGKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAMMERMGG------FGSIIQLYPGGGPVRAATACFGF----PKSF 250 .9************************************************......***********************....6665 PP Gcd10p 263 lkkh 266 l+ sp|Q9UJA5-2|TRM6_HUMAN 251 LSGL 254 5533 PP == domain 2 score: -0.9 bits; conditional E-value: 0.35 Gcd10p 75 eslaeeeaeeeeeeeeardnreiiddg 101 e+ +++e+e+e++++ea + + d++ sp|Q9UJA5-2|TRM6_HUMAN 293 EEKQASEQENEDSMAEAPESNHPEDQE 319 222222222222222222222222222 PP >> tr|B4DUV6|B4DUV6_HUMAN cDNA FLJ51319, highly similar to tRNA (adenine-N(1)-)-methyltransferase non-catalytic subunit # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 76.7 0.2 8.1e-25 8e-21 173 266 .. 2 84 .. 1 152 [. 0.74 2 ? -1.5 1.5 0.54 5.3e+03 75 101 .. 123 149 .. 84 187 .. 0.50 Alignments for each domain: == domain 1 score: 76.7 bits; conditional E-value: 8.1e-25 Gcd10p 173 llekkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaep 259 +++ ++++ki ++r +tl+++l+l+n++ag++++v+++++GL+++a++eRmg+ g+i++++++ p ++ +f++ p tr|B4DUV6|B4DUV6_HUMAN 2 YYA-REPGKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAMMERMGG------FGSIIQLYPGGGPVRAATACFGF----P 77 566.8************************************************......***********************....6 PP Gcd10p 260 ehplkkh 266 +++l+ tr|B4DUV6|B4DUV6_HUMAN 78 KSFLSGL 84 6655533 PP == domain 2 score: -1.5 bits; conditional E-value: 0.54 Gcd10p 75 eslaeeeaeeeeeeeeardnreiiddg 101 e+ +++e+e+e++++ea + + d++ tr|B4DUV6|B4DUV6_HUMAN 123 EEKQASEQENEDSMAEAPESNHPEDQE 149 222222222222222222222222222 PP >> tr|H0YJ17|H0YJ17_HUMAN Neuroguidin (Fragment) OS=Homo sapiens GN=NGDN PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 15.9 0.3 2.6e-06 0.026 48 125 .. 61 168 .. 46 171 .. 0.77 Alignments for each domain: == domain 1 score: 15.9 bits; conditional E-value: 2.6e-06 Gcd10p 48 ltfeildkreeeeks....rlrvvpaaeleaeslaeeeaeeeeeeee..........................ardnreiiddgarQ 104 lt+ ildk + + + lr+v+ ++l +e+ +e+eae++++e++ a+ ++ ++ ++r tr|H0YJ17|H0YJ17_HUMAN 61 LTHLILDKASGGSLQghdaVLRLVEIRTLSSEDEEEDEAEDDQSEASgkksvkgvskkyvpprlvpvhydeteAEREKKRLERAKRR 147 5555666655544445555677777766666666666655555544467789*******************99888999999***** PP Gcd10p 105 kLtkeeIeeLKkegasagkei 125 +L+++ I eLK++ +a +ei tr|H0YJ17|H0YJ17_HUMAN 148 ALSSSVIRELKEQYSDAPEEI 168 *************99999987 PP >> tr|Q5T0V7|Q5T0V7_HUMAN Dystonin (Fragment) OS=Homo sapiens GN=DST PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.1 0.1 3.9e-05 0.39 170 212 .. 632 674 .. 618 682 .. 0.84 Alignments for each domain: == domain 1 score: 12.1 bits; conditional E-value: 3.9e-05 Gcd10p 170 leyllekkdaqkilelreetlglllslanvhaggryLvvDdtg 212 ++y+ +++ +q++ +lr+e ++l + v ++gr+L ++t+ tr|Q5T0V7|Q5T0V7_HUMAN 632 GKYYQADQLVQRVAKLRDEIMALRNECSSVYSKGRILTTEQTK 674 5788888999******************************997 PP >> tr|B7ZKW0|B7ZKW0_HUMAN ATP1B4 protein OS=Homo sapiens GN=ATP1B4 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.7 0.3 5.2e-05 0.51 46 90 .. 15 71 .. 4 140 .. 0.74 Alignments for each domain: == domain 1 score: 11.7 bits; conditional E-value: 5.2e-05 Gcd10p 46 ygltfeildkreeeeks............rlrvvpaaeleaeslaeeeaeeeeeeee 90 y +++ d +e ++++ r+ vvp++e+e+e+ ++ee+eeee+eee tr|B7ZKW0|B7ZKW0_HUMAN 15 YSYRYRLDDPDEANQNYladeeeeaeeeaRVTVVPKSEEEEEEEEKEEEEEEEKEEE 71 55555555555444444555555555556*******999999999888887777665 PP >> sp|Q9UN42|AT1B4_HUMAN Protein ATP1B4 OS=Homo sapiens GN=ATP1B4 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.3 1.0 0.00013 1.3 46 90 .. 15 71 .. 4 108 .. 0.70 Alignments for each domain: == domain 1 score: 10.3 bits; conditional E-value: 0.00013 Gcd10p 46 ygltfeildkreeeeks............rlrvvpaaeleaeslaeeeaeeeeeeee 90 y +++ d +e ++++ r+ vvp++e+e+e+ ++ee+eeee+eee sp|Q9UN42|AT1B4_HUMAN 15 YSYRYRLDDPDEANQNYladeeeeaeeeaRVTVVPKSEEEEEEEEKEEEEEEEKEEE 71 55555555555444444555555555556*******999999999888887777655 PP >> sp|Q9UN42-2|AT1B4_HUMAN Isoform B of Protein ATP1B4 OS=Homo sapiens GN=ATP1B4 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.2 1.1 0.00015 1.4 46 90 .. 15 71 .. 4 105 .. 0.70 Alignments for each domain: == domain 1 score: 10.2 bits; conditional E-value: 0.00015 Gcd10p 46 ygltfeildkreeeeks............rlrvvpaaeleaeslaeeeaeeeeeeee 90 y +++ d +e ++++ r+ vvp++e+e+e+ ++ee+eeee+eee sp|Q9UN42-2|AT1B4_HUMAN 15 YSYRYRLDDPDEANQNYladeeeeaeeeaRVTVVPKSEEEEEEEEKEEEEEEEKEEE 71 55555555555444444555555555556*******999999999888887776655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 88473 (35070517 residues searched) Passed MSV filter: 4768 (0.0538921); expected 1769.5 (0.02) Passed bias filter: 1798 (0.0203226); expected 1769.5 (0.02) Passed Vit filter: 150 (0.00169543); expected 88.5 (0.001) Passed Fwd filter: 9 (0.000101726); expected 0.9 (1e-05) Initial search space (Z): 88473 [actual number of targets] Domain search space (domZ): 9 [number of targets reported over threshold] # CPU time: 0.87u 0.03s 00:00:00.90 Elapsed: 00:00:00.28 # Mc/sec: 36279.85 // [ok]