# 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_1395.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1395 [L=270] Description: # 1502661 # 1503470 # 1 # ID=1_1395;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.254 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.2e-08 35.1 0.0 2.8e-08 33.8 0.0 1.6 1 LCM Leucine carboxyl methyltransferase 0.0012 18.0 0.0 0.0017 17.5 0.0 1.2 1 Methyltransf_2 O-methyltransferase Domain annotation for each model (and alignments): >> LCM Leucine carboxyl methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 33.8 0.0 3.4e-12 2.8e-08 52 185 .. 44 170 .. 15 173 .. 0.80 Alignments for each domain: == domain 1 score: 33.8 bits; conditional E-value: 3.4e-12 CCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTTESEEEEET-TT--HHHHHHHTT...TEEEEEEE-HHHHHHHHHHHHHSHH CS LCM 52 ekkekkapelnrgiavRtrffddllkaalsaagirqvViLGAGlDtrayRlpwla...glkvfEvDqpevlefKretlkelga 131 ++ +++ n +a+R +d+++ +++++++i +V+LG Gl t R+ + +++++ +D pev+e++++ e + 23726_FusoPortal_Gene_1395 44 GISKGSFEYTNMAYAARYYNMDKMIIKFIKEHKISNIVLLGIGLETAYDRITQKCglgEVNYYGIDLPEVIEIRKKYFGERK- 125 3445678889999****************889**************999999988899***************998885444. PP HHHEECSEEE--TTSHHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHH CS LCM 132 epkaelryvpvDlrdddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalle 185 + ++++ D+ + +w+e+ +d++ pt+ ++ Gv Y ++++ + ++ 23726_FusoPortal_Gene_1395 126 ----QETLIAGDMFEMKWKEQ-----IDTSIPTLLIVSGVFQYFFEDKIIEFIK 170 ....45678999888888765.....6778*************98887766665 PP >> Methyltransf_2 O-methyltransferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.5 0.0 2e-07 0.0017 64 168 .. 66 176 .. 42 183 .. 0.73 Alignments for each domain: == domain 1 score: 17.5 bits; conditional E-value: 2e-07 HHHHHSS-GGG-SEEEEET-TTSHHHHHHHHH-TT..-EEEEEE-HHHHHHHHHHS-----TTEEEEES-TTTS---.--S.. CS Methyltransf_2 64 kkvleafdleglkslvDvGGGtGalakaivsayPt..ikvivfDlpeVvekAkehlsageeervefvgGdffkdslP.eaD.. 141 k +++ ++ ++++++v +G G+ + i+++ ++ DlpeV+e k++ +ge+++ + ++Gd+f+ + + D 23726_FusoPortal_Gene_1395 66 KMIIKFIKEHKISNIVLLGIGLETAYDRITQKCGLgeVNYYGIDLPEVIEIRKKY--FGERKQETLIAGDMFEMKWKeQIDts 146 4456666788999999999999999999999976511555667************..89999999*******65433222211 PP ...EEEEESSGGGS-HHHHHHHHHHHHHH- CS Methyltransf_2 142 ...aillkrvlhDwsdekcvklLkrvkkal 168 +++++ v + ++k+ +++k++kk 23726_FusoPortal_Gene_1395 147 iptLLIVSGVFQYFFEDKIIEFIKNLKKEF 176 111455566666666677777777776655 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (270 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 508 (0.0303973); expected 334.2 (0.02) Passed bias filter: 318 (0.0190282); expected 334.2 (0.02) Passed Vit filter: 30 (0.00179512); 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.25u 0.24s 00:00:00.49 Elapsed: 00:00:00.32 # Mc/sec: 2452.81 // [ok]