# 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/MafB19-deam.hmm.txt # target sequence database: proteomes/Anolis_carolinensis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MafB19-deam [M=138] Accession: PF14437.5 Description: MafB19-like deaminase 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.3e-20 73.6 0.5 1.8e-20 73.1 0.5 1.2 1 tr|R4GBK4|R4GBK4_ANOCA Uncharacterized protein OS=Anolis car 1.1e-14 54.4 0.0 1.3e-14 54.2 0.0 1.1 1 tr|H9GK09|H9GK09_ANOCA dCMP deaminase OS=Anolis carolinensis 5.8e-10 39.1 0.0 1.1e-06 28.5 0.0 2.2 2 tr|H9GW10|H9GW10_ANOCA Adenosine deaminase, tRNA specific 3 ------ inclusion threshold ------ 0.067 13.0 0.0 0.11 12.3 0.0 1.4 1 tr|G1KZ85|G1KZ85_ANOCA Apolipoprotein B mRNA editing enzyme Domain annotation for each sequence (and alignments): >> tr|R4GBK4|R4GBK4_ANOCA Uncharacterized protein OS=Anolis carolinensis OX=28377 GN=adat2 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 73.1 0.5 4e-24 1.8e-20 9 107 .. 125 235 .. 119 270 .. 0.85 Alignments for each domain: == domain 1 score: 73.1 bits; conditional E-value: 4e-24 MafB19-deam 9 kalaevekatkikedkigavivkdkkviargsnekelnakataHaEiaalqqAakklkskrldd......atlyvtlepCdmCagai 89 + +ae+++a++ +e+++g++ v +++vi++g+ne + ++at HaE+ a+ qA + ++++ + ++lyvt+epC mCa a+ tr|R4GBK4|R4GBK4_ANOCA 125 ETAAEAKEALENGEVPVGCLMVYNNEVIGKGKNEVNETKNATRHAEMVAIDQAIDWCHKQKKRMeevflhTVLYVTVEPCIMCAAAL 211 56899************************************************9988888776434444489*************** PP MafB19-deam 90 vqaglkklvygagnpks......l 107 +++ +vyg++n+ + l tr|R4GBK4|R4GBK4_ANOCA 212 RMMKIPLVVYGCQNERFggcgsvL 235 *************99887774432 PP >> tr|H9GK09|H9GK09_ANOCA dCMP deaminase OS=Anolis carolinensis OX=28377 GN=DCTD PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 54.2 0.0 2.8e-18 1.3e-14 3 103 .. 31 147 .. 29 182 .. 0.84 Alignments for each domain: == domain 1 score: 54.2 bits; conditional E-value: 2.8e-18 MafB19-deam 3 kekafrkalaevekatkikedkigaviv.kdkkviargsn...................ekelnakataHaEiaalqqAakklkskr 69 + ++f+ +++ +++k + ++ga iv ++k++++g+n + +++ HaE +a++ +++ tr|H9GK09|H9GK09_ANOCA 31 WPEYFMAVAFLSAQRSKDPHSQVGACIVnSENKIVGIGYNgmpngcsddslpwtrtaesKLDTKYPYVCHAELNAIMNK---NSA-D 113 678999******************9666589************************99776666777889*******974...444.4 PP MafB19-deam 70 lddatlyvtlepCdmCagaivqaglkklvygagn 103 +++ +yv l pC+ Ca++i+qag+k++++ +++ tr|H9GK09|H9GK09_ANOCA 114 VKGCSMYVALFPCNECAKLIIQAGIKEVIFMSDK 147 79***************************99984 PP >> tr|H9GW10|H9GW10_ANOCA Adenosine deaminase, tRNA specific 3 OS=Anolis carolinensis OX=28377 GN=ADAT3 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 8.4 0.0 0.00037 1.7 6 41 .. 178 215 .. 173 252 .. 0.81 2 ! 28.5 0.0 2.3e-10 1.1e-06 75 124 .. 269 327 .. 263 335 .. 0.81 Alignments for each domain: == domain 1 score: 8.4 bits; conditional E-value: 0.00037 MafB19-deam 6 afrkalaevekatkikedkigaviv..kdkkviargsn 41 ++ a++ +++++k++ +++gav+v +kv+a+g+ tr|H9GW10|H9GW10_ANOCA 178 HMELAIQAAQQGAKQGMRPVGAVVVdpSSGKVLAVGHD 215 56677888999***********999556789*****97 PP == domain 2 score: 28.5 bits; conditional E-value: 2.3e-10 MafB19-deam 75 lyvtlepCdmCagaivqaglkklvygagnpks.....lnkLkvkekngrvelv....ed 124 +y t+epC mCa +v +++ ++ yg+ p + l + n+r+e++ e tr|H9GW10|H9GW10_ANOCA 269 MYLTREPCAMCAMSLVHSRIERVFYGVPSPHGalgtaLHIHSRRDLNHRYEVFrgvlEG 327 9***************************9999986666666666789999998555444 PP >> tr|G1KZ85|G1KZ85_ANOCA Apolipoprotein B mRNA editing enzyme catalytic subunit 2 OS=Anolis carolinensis OX=28377 GN=A # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.3 0.0 2.4e-05 0.11 47 92 .. 92 137 .. 64 157 .. 0.81 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 2.4e-05 MafB19-deam 47 akataHaEiaalqqAakklkskrlddatlyvtlepCdmCagaivqa 92 ++a aHaE a ++ + ++t yv pC Ca iv+a tr|G1KZ85|G1KZ85_ANOCA 92 EHAAAHAEEAFFNTILPTCEPGFRYEVTWYVSSSPCVSCAERIVKA 137 467788887666655555555556699****************998 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (138 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 600 (0.0323834); expected 370.6 (0.02) Passed bias filter: 356 (0.0192142); expected 370.6 (0.02) Passed Vit filter: 31 (0.00167314); expected 18.5 (0.001) Passed Fwd filter: 4 (0.000215889); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.06 # Mc/sec: 22532.07 // [ok]