# hmmscan :: search sequence(s) against a profile 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 sequence file: AA/23726_FusoPortal_Gene_205.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_205 [L=315] Description: # 218984 # 219928 # 1 # ID=1_205;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.296 Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.5e-118 393.4 0.0 8.5e-118 393.3 0.0 1.0 1 Methyltransf_5 MraW methylase family 5e-05 23.0 0.2 0.0023 17.6 0.2 2.2 1 Methyltransf_31 Methyltransferase domain 0.0019 18.2 0.0 0.016 15.3 0.0 2.0 1 FtsJ FtsJ-like methyltransferase ------ inclusion threshold ------ 0.014 14.7 0.0 0.082 12.2 0.0 1.9 2 Arylsulfotran_2 Arylsulfotransferase (ASST) 0.019 14.5 0.0 2 7.9 0.0 2.3 2 MTS Methyltransferase small domain 0.058 14.0 0.0 0.34 11.5 0.0 2.1 2 Methyltransf_25 Methyltransferase domain Domain annotation for each model (and alignments): >> Methyltransf_5 MraW methylase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 393.3 0.0 3e-121 8.5e-118 3 308 .. 8 311 .. 6 312 .. 0.99 Alignments for each domain: == domain 1 score: 393.3 bits; conditional E-value: 3e-121 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TT.EEEEEES-HHHHHHHTT.....---EEEEES-GGGHHHHH CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlsel.rligidrDksalarakerLkefekrltlvkgnfkelkeal 85 H+pVL e++++l +++dg+yiD+TlG+g Hse +Le+ls++ li+id+D++a+ +k+rL++f + +++kgnf++++++ 23726_FusoPortal_Gene_205 8 YHIPVLYYETLDNLVINPDGTYIDCTLGGGSHSEGILERLSDKgLLISIDQDTNAIVYSKKRLEKFGSKWKVFKGNFENIDTIA 91 6******************************************99*************************************** PP HHTTTT-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTSS--HHHHHHH--HHHHHHHHHHHH--TTHHHHHHHHHHH CS Methyltransf_5 86 aeanvtkvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsaeltAaevlnrlseeeLakilkeygeekfakriakaivek 169 a+v+kvdgil+D+GvSs qld++erGFs++ d LDmRm+++++ +A++v+n++see+L ki++eygee++a++iak ive+ 23726_FusoPortal_Gene_205 92 YMAGVDKVDGILMDIGVSSKQLDDPERGFSYRYDVKLDMRMNTDQKISAYDVVNTYSEEQLSKIIFEYGEERHARKIAKLIVEE 175 ************************************************************************************ PP HHHS---BHHHHHHHHHHHS-------..-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHH CS Methyltransf_5 170 rkkkpiettkelaevikkvvpakekrkkihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkkvf 253 rk +pie+t++l +ik+++p+ + hpa+++FQA+RI+vN+Ele le+++++a+ell+ ggrla+i+FhSledrivk+ f 23726_FusoPortal_Gene_205 176 RKSSPIEKTSDLIALIKRAYPERAS---KHPAKKTFQAIRIEVNRELEVLENAMSKATELLKVGGRLAIITFHSLEDRIVKNKF 256 *********************8776...59****************************************************** PP HHTT-------------------EEEEE.E.E---HHHHHH-GGGTT-EEEEEEE CS Methyltransf_5 254 keaskkeevPaglpvaseelkaklklitkkpilpseeeikeNprsrsakLrviek 308 k++++ +++P+++p++ ++ ++k+++itkkpi+p ++e+k+N+r++s+kLr++e+ 23726_FusoPortal_Gene_205 257 KDLATACKCPKDIPICVCGGVKKFEVITKKPIIPIDDELKNNNRAHSSKLRILER 311 *****************************************************98 PP >> Methyltransf_31 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.6 0.2 8.4e-07 0.0023 4 113 .. 26 245 .. 23 285 .. 0.66 Alignments for each domain: == domain 1 score: 17.6 bits; conditional E-value: 8.4e-07 Methyltransf_31 4 glrvlDlGCGtGeltfelaeelgpkaevvGiDiseeaielArenaqklgldnveFeqgdieelekl.led..ekfDlilsn... 81 + +++D+ G G + + e+l +k+ +++iD + +ai +++++++k g++ ++ +g++e++++ + +k+D il + 23726_FusoPortal_Gene_205 26 DGTYIDCTLGGGSHSEGILERLSDKGLLISIDQDTNAIVYSKKRLEKFGSK-WKVFKGNFENIDTIaYMAgvDKVDGILMDigv 108 5678999999999***9**9******************************8.9***********98777677888877777666 PP Methyltransf_31 82 .................................................................................... 81 23726_FusoPortal_Gene_205 109 sskqlddpergfsyrydvkldmrmntdqkisaydvvntyseeqlskiifeygeerharkiakliveerksspiektsdlialik 192 666666666666666666666666666666666666666666666666666666666665555555555555555554444444 PP Methyltransf_31 82 .....................evlnaladpdkileeiervlkpggilliedpr 113 ev l +++++ +++++lk gg+l+i +++ 23726_FusoPortal_Gene_205 193 rayperaskhpakktfqairiEVNRELEVLENAMSKATELLKVGGRLAIITFH 245 44444444444443332222222244444444444666666666666666666 PP >> FtsJ FtsJ-like methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.3 0.0 5.6e-06 0.016 19 99 .. 23 110 .. 15 266 .. 0.83 Alignments for each domain: == domain 1 score: 15.3 bits; conditional E-value: 5.6e-06 S-TT-EEEEES-TT-HHHHHHTTT...-SEEEEEESS--.........---.TT-EEEE--TTSSSHHHHHHHHHHHSSEEEEE CS FtsJ 19 fkpgktvvDlgaaPGgwsqvllqrlgakgkvvavDlqpm.........apl.pgvttiqgditdpetiekikealserkvdlvv 92 ++p++t +D+ + G+ s+ +l+rl kg ++++D+ + +++ ++ ++++g++ + +ti + kvd ++ 23726_FusoPortal_Gene_205 23 INPDGTYIDCTLGGGSHSEGILERLSDKGLLISIDQDTNaivyskkrlEKFgSKWKVFKGNFENIDTIAYMAGV---DKVDGIL 103 5788999******************999********8886788764443337788999***********99963...58***** PP E------ CS FtsJ 93 sDgapdv 99 +D++ + 23726_FusoPortal_Gene_205 104 MDIGVSS 110 **99532 PP >> Arylsulfotran_2 Arylsulfotransferase (ASST) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 0.3 0.0 0.12 3.4e+02 94 115 .. 49 70 .. 21 81 .. 0.71 2 ? 12.2 0.0 3e-05 0.082 189 252 .. 111 186 .. 83 194 .. 0.80 Alignments for each domain: == domain 1 score: 0.3 bits; conditional E-value: 0.12 Arylsulfotran_2 94 ldsliyeidietnevlfrwSal 115 +l++ id +tn++++++ l 23726_FusoPortal_Gene_205 49 DKGLLISIDQDTNAIVYSKKRL 70 3456777888888888776555 PP == domain 2 score: 12.2 bits; conditional E-value: 3e-05 Arylsulfotran_2 189 fel..dnstfsrQHdariveetedkivislfnNentdt..........nerepttglildvdlqnkevtllrklkd 252 +l ++ fs+ +d+++ + +++++is+++ nt + ++e++ +i ++ ++++++++++k++d 23726_FusoPortal_Gene_205 111 KQLddPERGFSYRYDVKLDMRMNTDQKISAYDVVNTYSeeqlskiifeYGEERHARKIAKLIVEERKSSPIEKTSD 186 455778889***********************97777666777788887888888899999999999988888875 PP >> MTS Methyltransferase small domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 4.2 0.0 0.0092 26 26 89 .. 20 84 .. 13 102 .. 0.83 2 ? 7.9 0.0 0.0007 2 112 160 .. 216 263 .. 210 267 .. 0.82 Alignments for each domain: == domain 1 score: 4.2 bits; conditional E-value: 0.0092 HCHCCTCSEEEEET-TTSHHHHHHHHHTTC.E.EEEEESBHHHHHHHHHHHHHTTGCEEEEE-ST CS MTS 26 tlpkelsgkvlDlGcGaGvlgaalakespeae.vtlvdvdarAlesakenleanklevevvasdv 89 +l + +g+ +D G G + ++++ ++ + +d+d++A+ +k+ le+ + + +v++ + 23726_FusoPortal_Gene_205 20 NLVINPDGTYIDCTLGGGSHSEGILERLSDKGlLISIDQDTNAIVYSKKRLEKFGSKWKVFKGNF 84 4455578888999999999998888877776658899***************9999988888765 PP == domain 2 score: 7.9 bits; conditional E-value: 0.0007 HCHHHHHHHHHHHHHHHEEEEEEEEEEEETTSCHHHHHHHHHS--EEEE CS MTS 112 atelavaeeliaeAkeaLkegGeLliVanrklgyeslleelfgnvevva 160 + el+v e+ + +A e Lk gG+L i++ + l+ +++++++f+++++++ 23726_FusoPortal_Gene_205 216 NRELEVLENAMSKATELLKVGGRLAIITFHSLE-DRIVKNKFKDLATAC 263 5678899999*****************999887.567788888887766 PP >> Methyltransf_25 Methyltransferase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 0.00012 0.34 3 60 .. 31 89 .. 30 121 .. 0.83 2 ? -0.4 0.0 0.63 1.8e+03 77 97 .] 217 237 .. 192 237 .. 0.71 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 0.00012 EET-TTSHHHHHHHCTS.TSSEEEEEESSHHHHHHHHHHHHHTTTTEEEEESBCCGBSS CS Methyltransf_25 3 dlGcGtGrltlalakrl.gkakvtgvDiseemleaareraeeeglnvefvqgdaedlpf 60 d+ G G +++ +++rl k+ ++ +D + +++ + ++r+e+ g ++++++g++e+++ 23726_FusoPortal_Gene_205 31 DCTLGGGSHSEGILERLsDKGLLISIDQDTNAIVYSKKRLEKFGSKWKVFKGNFENIDT 89 77889***********9549********************99999*********99873 PP == domain 2 score: -0.4 bits; conditional E-value: 0.63 SSHHHHHHHHHHHHHHEEEEE CS Methyltransf_25 77 lseedleaalreiarvLkpgG 97 + e le+a++++ ++Lk gG 23726_FusoPortal_Gene_205 217 RELEVLENAMSKATELLKVGG 237 344578899999999999998 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (315 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 656 (0.0392532); expected 334.2 (0.02) Passed bias filter: 467 (0.027944); expected 334.2 (0.02) Passed Vit filter: 38 (0.00227382); expected 16.7 (0.001) Passed Fwd filter: 6 (0.000359023); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.27u 0.23s 00:00:00.50 Elapsed: 00:00:00.31 # Mc/sec: 2953.92 // [ok]