# 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_1930.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1930 [L=323] Description: # 2072239 # 2073207 # -1 # ID=1_1930;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.292 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 9.3e-149 494.7 0.8 1e-148 494.6 0.8 1.0 1 ALAD Delta-aminolevulinic acid dehydratase ------ inclusion threshold ------ 0.022 14.0 0.2 0.055 12.7 0.0 1.7 2 MTHFR Methylenetetrahydrofolate reductase Domain annotation for each model (and alignments): >> ALAD Delta-aminolevulinic acid dehydratase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 494.6 0.8 1.2e-152 1e-148 2 317 .] 3 319 .. 2 319 .. 0.99 Alignments for each domain: == domain 1 score: 494.6 bits; conditional E-value: 1.2e-152 --GGGGGSSHHHHHHC-S----GGGEEEEEEEETT-SS-EE-TTSTT-EEEEGGGHHHHHHHHHHTT--EEEEEEE-SS-BSS CS ALAD 2 krlrrlrksealrelvretrlsasdlilplfvkegeeekeeiksmpgvyrlsvdelvkeveelvklgikavllFgvpeekdee 84 +r+rrlr++ +relv++ ++++s+li+plfv++ge+ k+ei+smp+++r+s+d+l +e+++l+klgi+++llFg+p++kde 23726_FusoPortal_Gene_1930 3 TRTRRLRRNFLTRELVKNISIEKSSLIYPLFVCDGENIKSEIESMPQQFRYSLDRLNEELDDLLKLGINNILLFGIPNHKDEL 85 79********************************************************************************* PP -GGGGSTTSHHHHHHHHHHHH-.TTSEEEEEE-STTTBSS-SSSEEESS.EEHHHHHHHHHHHHHHHHHCT-SEEEE-S--TT CS ALAD 85 gsealnedglvqrairaikeev.pelvvitDvclceytshGhcgileekevdndetlellakiavsqaeaGadivapSdmmDg 166 gs+a++e+g+vqr++r+i++++ ++++++tDvc+ceytshGhcgil++++vdndetle++akia+s+a+aGadi+apSdmmDg 23726_FusoPortal_Gene_1930 86 GSQAYDENGIVQRTVRQIRRDYqDKFLIVTDVCMCEYTSHGHCGILHNHDVDNDETLEYIAKIALSHAKAGADIIAPSDMMDG 168 ********************9945599******************************************************** PP HHHHHHHHHHHTT-TT-EEB---EEB--CGGHHHHHHHT---S.S--TTTSB-TT-HHHHHHHHHHHHHTT-SBEEEES-TTG CS ALAD 167 rvaaireaLdeagfeevailsYsaKyaSafYgPFrdaaesapsfgdrktYqmdpansrealreveldleeGADivmvKPalay 249 r+++ire Ld+++f++++i+sYs+Ky+Sa+YgPFrdaa+sapsfgdrktYqmd+++ ++ +reve d++eGAD++mvKPa+ay 23726_FusoPortal_Gene_1930 169 RIGKIREVLDKNNFKDIPIMSYSVKYSSAYYGPFRDAADSAPSFGDRKTYQMDFRSYNNFYREVEADTQEGADFIMVKPAMAY 251 *********************************************************************************** PP HHHHHHHHTT-SS-EEEEE-HHHHHHHHHHHHTTSS-HHHHHHHHHHHHHHTT-SEEEETTHHHHHH. CS ALAD 250 lDiirevkeefelpvaayqvSGEYamikaaaengwldekkvvlEsllsikRAGadliltYfakevakl 317 lD+i++v++ ++lp++ay+vSGEY+m+kaaa+n+w+de+k+v+E++ +ikRAGad+i+tY+ak++ak+ 23726_FusoPortal_Gene_1930 252 LDVIKSVSKISNLPIVAYNVSGEYSMVKAAAKNNWIDEQKIVMENMYAIKRAGADIIITYHAKDIAKW 319 ******************************************************************98 PP >> MTHFR Methylenetetrahydrofolate reductase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 6.6e-06 0.055 82 166 .. 53 146 .. 41 163 .. 0.68 2 ? -2.9 0.0 0.36 3e+03 170 178 .. 302 310 .. 277 318 .. 0.77 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 6.6e-06 MTHFR 82 askeeldkalekakelGirnilALrGDppkdedkweekeee.fkyaldLvrlireeygdlf....digvaaYP.evhp..... 153 +s+ l+++l++ +lGi nil L G p ++ d++ ++ ++ + vr ir++y d+f d+++ Y ++h+ 23726_FusoPortal_Gene_1930 53 YSLDRLNEELDDLLKLGINNIL-LFGIPNHK-DELGSQAYDeNGIVQRTVRQIRRDYQDKFlivtDVCMCEYTsHGHCgilhn 133 6778999**************9.77876665.56777776636777889************3333345556766777788833 PP MTHFR 154 eaksleadlkylk 166 ++ + ++ l+y+ 23726_FusoPortal_Gene_1930 134 HDVDNDETLEYIA 146 3333334455543 PP == domain 2 score: -2.9 bits; conditional E-value: 0.36 MTHFR 170 daGAdfiit 178 +aGAd iit 23726_FusoPortal_Gene_1930 302 RAGADIIIT 310 688888888 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (323 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 652 (0.0390139); expected 334.2 (0.02) Passed bias filter: 485 (0.0290211); expected 334.2 (0.02) Passed Vit filter: 37 (0.00221398); 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.23s 00:00:00.49 Elapsed: 00:00:00.30 # Mc/sec: 3129.90 // [ok]