# 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_1151.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1151 [L=236] Description: # 1253621 # 1254328 # 1 # ID=1_1151;partial=00;start_type=TTG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.251 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.8e-57 193.6 1.5 3.1e-57 193.4 1.5 1.0 1 Methyltrans_RNA RNA methyltransferase 0.0029 17.4 4.7 0.35 10.7 0.3 2.7 3 CryBP1 CryBP1 protein Domain annotation for each model (and alignments): >> Methyltrans_RNA RNA methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 193.4 1.5 3.7e-61 3.1e-57 3 223 .. 17 229 .. 15 232 .. 0.96 Alignments for each domain: == domain 1 score: 193.4 bits; conditional E-value: 3.7e-61 EEHHHHHHHHTTTT--TT-EEEEEESSSEEEEEEEEEEESSEEEEEEEEEEE---SSSSEEEEEEE--SSHHHHHHHHHHHHT CS Methyltrans_RNA 3 ltdeeahHlvkVlRlkegdeiklfdgdgkeyeaeieeiskkkvtvkilekeeenkespekitlaqalpkgdklelilqkatEl 85 ++ +++H+++ R +gd+i+++d +ey++eie+i++k++++kilek e + + +++ ++++kgdk++l++qk+tEl 23726_FusoPortal_Gene_1151 17 IDAGDINHIKNAFRKVKGDKIRAVD-GANEYLCEIEKIEDKEIKLKILEKIEDKFSLDIEVDAGISILKGDKMDLTIQKLTEL 98 677899*******************.7789***************************************************** PP T-SEEEEEE-TTSSS...HHHHHHHHHHHHHHHHHHHHHHT-SS--EEEEEEEHHHHHHHCSCCEEEEEEETTSSEBGGGSSH CS Methyltrans_RNA 86 GvseivplisersvvkldekkkdkklerwqkilieAaeQsgrtrlpevkepislkelleelkeadlklildeeekksleelse 168 G+++i+p+ +r+vvkl dkk++rw++i +eA++Q++ + ++e++++++l +k++dl+l+ +e+ek+ ++ 23726_FusoPortal_Gene_1151 99 GINKIIPISVKRCVVKL-----DKKKDRWDTISKEALKQCQGVVPTVIDEIKKINKLN--FKDYDLILVPYENEKEIFIKEIL 174 **************999.....9*********************************96..7899999**********999999 PP CCHTHHCSEEEEEE--TT---HHHHHHHHHTTCEEEB-SSS---HHHHHHHHHHH CS Methyltrans_RNA 169 kkkelkevkvlliiGPEGGfsekEiellkeagfkpvsLgpriLRteTAalaalsa 223 + ++k k+l+iiG EGGf+++Ei+ lke+g+k++sLg+riLR+eTAa+++ ++ 23726_FusoPortal_Gene_1151 175 RDLKIKPSKILYIIGAEGGFEKEEIDYLKENGAKIISLGKRILRAETAAIVTGGV 229 99999999******************************************98665 PP >> CryBP1 CryBP1 protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 8.2 0.3 0.00024 2 58 95 .. 11 48 .. 4 59 .. 0.75 2 ! 10.7 0.3 4.2e-05 0.35 3 65 .. 55 115 .. 53 120 .. 0.90 3 ? -1.2 0.0 0.19 1.6e+03 98 145 .. 128 179 .. 116 184 .. 0.72 Alignments for each domain: == domain 1 score: 8.2 bits; conditional E-value: 0.00024 CryBP1 58 vdvcvvkivgaieyilnaypavrgdkygrgegdievip 95 d ++v g i++i na+ v+gdk +g e + 23726_FusoPortal_Gene_1151 11 DDFILVIDAGDINHIKNAFRKVKGDKIRAVDGANEYLC 48 5677778899****************988887655544 PP == domain 2 score: 10.7 bits; conditional E-value: 4.2e-05 CryBP1 3 dteveftcvtcvpesfavdiesleeirfifntdklyccldktteecnipgigdievdvcvvki 65 d e+++ ++ + + f++die i i++ dk+ ++k te i i i v cvvk+ 23726_FusoPortal_Gene_1151 55 DKEIKLKILEKIEDKFSLDIEVDAGIS-ILKGDKMDLTIQKLTE-LGINKIIPISVKRCVVKL 115 789999999999*********988886.889**********996.6899************97 PP == domain 3 score: -1.2 bits; conditional E-value: 0.19 CryBP1 98 eanvccegt..tfvdnv..vgafavddpeliedpceeldfcdiegvltdisi 145 ea c+g t++d++ + + d +li p e+ + i+ +l d+ i 23726_FusoPortal_Gene_1151 128 EALKQCQGVvpTVIDEIkkINKLNFKDYDLILVPYENEKEIFIKEILRDLKI 179 5556666643367776633556777888899999888888888888888776 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (236 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 765 (0.0457755); expected 334.2 (0.02) Passed bias filter: 419 (0.0250718); expected 334.2 (0.02) Passed Vit filter: 27 (0.00161561); expected 16.7 (0.001) Passed Fwd filter: 2 (0.000119674); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.26u 0.24s 00:00:00.50 Elapsed: 00:00:00.32 # Mc/sec: 2143.94 // [ok]