Converged !!!

E-value threshold for inclusion in PSI-Blast iteration 2: 0.001 
E-value threshold for inclusion in PSI-Blast iteration 3: 




 Results of PSI-Blast iteration 1

Legend:
         - means that  the alignment score was below the threshold on the previous iteration 			
         - means that  the alignment was checked on the previous iteration 			
			
Sequences with E-value BETTER than threshold 

                                                                     Score    E
Sequences producing significant alignments:                          (bits) Value


 sp|P41819|DIM1_YEAST  DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLM...   201  2e-51

 gnl|PID|e214702  (Z68294) dimethylase [Kluyveromyces lactis]          190  3e-48

 sp|P78697|DIM1_KLULA  DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLM...   190  3e-48

 gi|4050050  (AF102147) putative dimethyladenosine transferase [H...   186  5e-47

 gi|3859994  (AF091078) unknown [Homo sapiens]                         186  5e-47

 gnl|PID|e1321529  (AJ009761) putative dimethyladenosine transfer...   186  5e-47

 gi|2529685  (AC002535) putative dimethyladenosine transferase [A...   170  3e-42

 sp|Q09522|YQN1_CAEEL  HYPOTHETICAL 34.1 KD PROTEIN E02H1.1 IN CH...   168  1e-41


Sequences with E-value WORSE than threshold 

 gi|2648762  (AE000979) dimethyladenosine transferase (ksgA) [Arc...    33  0.72

 gnl|PID|d1031885  (AP000007) 290aa long hypothetical dimethylade...    33  0.72

 sp|P44749|KSGA_HAEIN  DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLM...    33  0.94

 sp|P37468|KSGA_BACSU  DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLM...    33  0.94

 sp|Q03986|ERMD_BACLI  RRNA ADENINE N-6-METHYLTRANSFERASE (MACROL...    32  1.6

 sp|P51659|DHB4_HUMAN  ESTRADIOL 17 BETA-DEHYDROGENASE 4 (17-BETA...    32  1.6

 sp|P06992|KSGA_ECOLI  DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLM...    31  2.1

 gnl|PID|d1017401  (D90899) dimethyladenosine transferase [Synech...    31  2.1

 sp|Q04720|ERMJ_BACAN  RRNA ADENINE N-6-METHYLTRANSFERASE (MACROL...    31  3.6

 pir||A47697  rRNA methyltransferase - Bacillus anthracis               31  3.6

 gi|987265  (U31089) Abl binding protein 3 [Homo sapiens]               31  3.6

 gi|2984099  (AE000757) dimethyladenosine transferase [Aquifex ae...    31  3.6

 gi|915311  (U23435) Abl interactor 2 [Homo sapiens]                    30  4.8

 gnl|PID|e224171  (X95632) Arg protein tyrosine kinase-binding pr...    30  4.8

 gnl|PID|e224172  (X95677) argBPIB gene product [Homo sapiens]          30  4.8

 gnl|PID|e307817  (Y11786) dimethyladenosine transferase [Rickett...    30  4.8

 gi|4155935  (AE001555) DIMETHYLADENOSINE TRANSFERASE [Helicobact...    30  4.8

 sp|P45438|ERMK_BACLI  RRNA ADENINE N-6-METHYLTRANSFERASE (MACROL...    30  6.2

 pir||B48474  glycoprotein B - feline herpesvirus 1 (fragment) >g...    30  6.2

 pir||A56602  glycoprotein B homolog precursor - feline herpesvir...    30  6.2

 gi|471723  (U01887) gB homolog [Meleagrid herpesvirus type 1]          30  6.2

 prf||1911192A  glycoprotein B [Feline herpesvirus 1]                   30  6.2

 gi|2688507  (AE001160) dimethyladenosine transferase (ksgA) [Bor...    29  8.2

 gnl|PID|e293782  (Y10669) lef-8 [Spodoptera littoralis nuclear p...    29  8.2

 gi|471720  (U01886) gB homolog [Gallid herpesvirus 2]                  29  8.2

 gi|1369852  (L46702) kinesin heavy chain-like protein [Solanum t...    29  8.2

 gi|2314603  (AE000643) 16S rRNA (adenosine-N6,N6-)-dimethyltrans...    29  8.2


Alignments
 sp|P41819|DIM1_YEAST DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLMETHIONINE-6-N',
           N'-ADENOSYL(RRNA) DIMETHYLTRANSFERASE) (18S RRNA
           DIMETHYLASE) >gi|1078347|pir||S47985 rRNA
           (adenine-N6,N6-)-dimethyltransferase (EC 2.1.1.-) DIM1 -
           yeast (Saccharomyces cerevisiae) >gi|601876 (L26480)
           dimethyladenosine transferase [Saccharomyces cerevisiae]
           >gi|1370549|gnl|PID|e246972 (Z73622) ORF YPL266w
           [Saccharomyces cerevisiae]
           Length = 318
           
 Score =  201 bits (505), Expect = 2e-51
 Identities = 120/120 (100%), Positives = 120/120 (100%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT
Sbjct: 198 PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 257

Query: 61  FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
           FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF
Sbjct: 258 FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 317


 gnl|PID|e214702 (Z68294) dimethylase [Kluyveromyces lactis]
           Length = 320
           
 Score =  190 bits (477), Expect = 3e-48
 Identities = 89/121 (73%), Positives = 103/121 (84%), Gaps = 1/121 (0%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+VESSVVR+EIKNPRPQVD+NEWDGLLRIVFVRKNRTI+AGFKSTT ++ILEKNYK 
Sbjct: 199 PPPKVESSVVRIEIKNPRPQVDFNEWDGLLRIVFVRKNRTIAAGFKSTTDLEILEKNYKA 258

Query: 61  FLAMNNEMVDDTKG-SMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIH 119
           FLA  + +   + G S+ + VKEKI+ VL ET L +KRAGKCDQ DFL+LLY FHQVGIH
Sbjct: 259 FLATQSAVPTTSSGDSLINEVKEKIEQVLSETGLAEKRAGKCDQTDFLKLLYGFHQVGIH 318

Query: 120 F 120
           F
Sbjct: 319 F 319


 sp|P78697|DIM1_KLULA DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLMETHIONINE-6-N',
           N'-ADENOSYL(RRNA) DIMETHYLTRANSFERASE) (18S RRNA
           DIMETHYLASE)
           Length = 320
           
 Score =  190 bits (477), Expect = 3e-48
 Identities = 89/121 (73%), Positives = 103/121 (84%), Gaps = 1/121 (0%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+VESSVVR+EIKNPRPQVD+NEWDGLLRIVFVRKNRTI+AGFKSTT ++ILEKNYK 
Sbjct: 199 PPPKVESSVVRIEIKNPRPQVDFNEWDGLLRIVFVRKNRTIAAGFKSTTDLEILEKNYKA 258

Query: 61  FLAMNNEMVDDTKG-SMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIH 119
           FLA  + +   + G S+ + VKEKI+ VL ET L +KRAGKCDQ DFL+LLY FHQVGIH
Sbjct: 259 FLATQSAVPTTSSGDSLINEVKEKIEQVLSETGLAEKRAGKCDQTDFLKLLYGFHQVGIH 318

Query: 120 F 120
           F
Sbjct: 319 F 319


 gi|4050050 (AF102147) putative dimethyladenosine transferase [Homo sapiens]
           Length = 313
           
 Score =  186 bits (467), Expect = 5e-47
 Identities = 59/120 (49%), Positives = 83/120 (69%), Gaps = 5/120 (4%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+VESSVVR+E KNP P +++ EWDGL+RI FVRKN+T+SA FKS+ V  +LEKNY+ 
Sbjct: 198 PPPKVESSVVRIEPKNPPPPINFQEWDGLVRITFVRKNKTLSAAFKSSAVQQLLEKNYRI 257

Query: 61  FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
             +++N ++ +        + +KI  +L  T   DKRA   D +DF+RLL+ F+  GIHF
Sbjct: 258 HCSVHNIIIPEDFS-----IADKIQQILTSTGFSDKRARSMDIDDFIRLLHGFNAEGIHF 312


 gi|3859994 (AF091078) unknown [Homo sapiens]
           Length = 125
           
 Score =  186 bits (467), Expect = 5e-47
 Identities = 59/120 (49%), Positives = 83/120 (69%), Gaps = 5/120 (4%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+VESSVVR+E KNP P +++ EWDGL+RI FVRKN+T+SA FKS+ V  +LEKNY+ 
Sbjct: 10  PPPKVESSVVRIEPKNPPPPINFQEWDGLVRITFVRKNKTLSAAFKSSAVQQLLEKNYRI 69

Query: 61  FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
             +++N ++ +        + +KI  +L  T   DKRA   D +DF+RLL+ F+  GIHF
Sbjct: 70  HCSVHNIIIPEDFS-----IADKIQQILTSTGFSDKRARSMDIDDFIRLLHGFNAEGIHF 124


 gnl|PID|e1321529 (AJ009761) putative dimethyladenosine transferase [Homo sapiens]
           Length = 157
           
 Score =  186 bits (467), Expect = 5e-47
 Identities = 59/120 (49%), Positives = 83/120 (69%), Gaps = 5/120 (4%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+VESSVVR+E KNP P +++ EWDGL+RI FVRKN+T+SA FKS+ V  +LEKNY+ 
Sbjct: 42  PPPKVESSVVRIEPKNPPPPINFQEWDGLVRITFVRKNKTLSAAFKSSAVQQLLEKNYRI 101

Query: 61  FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
             +++N ++ +        + +KI  +L  T   DKRA   D +DF+RLL+ F+  GIHF
Sbjct: 102 HCSVHNIIIPEDFS-----IADKIQQILTSTGFSDKRARSMDIDDFIRLLHGFNAEGIHF 156


 gi|2529685 (AC002535) putative dimethyladenosine transferase [Arabidopsis
           thaliana]
           Length = 353
           
 Score =  170 bits (426), Expect = 3e-42
 Identities = 61/157 (38%), Positives = 83/157 (52%), Gaps = 37/157 (23%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+V+SSVVR+E + P PQV+  EWDG LR+ F+RKN+T+ + FK  +V+ +LEKN+KT
Sbjct: 196 PPPKVDSSVVRIEPRRPGPQVNKKEWDGFLRVCFIRKNKTLGSIFKQKSVLSMLEKNFKT 255

Query: 61  FLAM-----------------------------NNEMVDDT--------KGSMHDVVKEK 83
             A+                             +NEM DD          G      KEK
Sbjct: 256 LQAVLASLQNNGEPALNTTSMDLGDQSMGMEDDDNEMDDDDMEMDEGEGDGGETSEFKEK 315

Query: 84  IDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
           +  VLKE    +KR+ K  Q +FL LL  F++ GIHF
Sbjct: 316 VMNVLKEGGFEEKRSSKLSQQEFLYLLSLFNKSGIHF 352


 sp|Q09522|YQN1_CAEEL HYPOTHETICAL 34.1 KD PROTEIN E02H1.1 IN CHROMOSOME II
           >gi|3875432|gnl|PID|e1344981 (Z47075) similar to rRNA
           methyltransferase [Caenorhabditis elegans]
           Length = 304
           
 Score =  168 bits (420), Expect = 1e-41
 Identities = 46/120 (38%), Positives = 73/120 (60%), Gaps = 9/120 (7%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           PPP+V+S+VVR+  KNP P     EW+GLLR+ F+RKN+T+ A F+ + V++++E N++ 
Sbjct: 192 PPPKVDSAVVRIAPKNPPPP----EWEGLLRLCFMRKNKTLMAIFRLSNVIEVIEDNFRK 247

Query: 61  FLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
             +  N+ +          +K+ I+  L  +  G+ RA K    DFL LL AF++  IHF
Sbjct: 248 VCSFKNKPIPKDLN-----MKKVIEETLTASGYGESRARKMRVEDFLALLLAFNKADIHF 302


 gi|2648762 (AE000979) dimethyladenosine transferase (ksgA) [Archaeoglobus
           fulgidus]
           Length = 244
           
 Score = 33.0 bits (73), Expect = 0.72
 Identities = 9/44 (20%), Positives = 21/44 (47%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAG 44
           PPP+VES++VR+  +      +   ++  +   F  + + +   
Sbjct: 163 PPPKVESAIVRIVPEPEVFVENRELFEKFVTFAFSMRRKRMGKI 206


 gnl|PID|d1031885 (AP000007) 290aa long hypothetical dimethyladenosine transferase
           [Pyrococcus horikoshii]
           Length = 290
           
 Score = 33.0 bits (73), Expect = 0.72
 Identities = 16/94 (17%), Positives = 43/94 (45%), Gaps = 8/94 (8%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKT 60
           P P+V+S++V +E +  +  ++    + L++ +F  + +T+    K +  + +L  +   
Sbjct: 200 PRPKVDSALVLIEPRKDKIVLN----ENLVKALFQHRRKTVPRALKDS--IHMLGVSKDE 253

Query: 61  FLAMNNEMVDDTKG--SMHDVVKEKIDTVLKETD 92
              + N +    K    ++    + I+  LK+  
Sbjct: 254 IRGIINNVPHSNKRVFQLYPEEVKDIEEYLKKHG 287


 sp|P44749|KSGA_HAEIN DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLMETHIONINE-6-N',
           N'-ADENOSYL(RRNA) DIMETHYLTRANSFERASE) (16S RRNA
           DIMETHYLASE) (HIGH LEVEL KASUGAMYCIN RESISTANCE PROTEIN
           KSGA) (KASUGAMYCIN DIMETHYLTRANSFERASE)
           >gi|1073774|pir||I64076 16S rRNA
           (adenosine-N6,N6-)-dimethyltransferase (ksgA) homolog -
           Haemophilus influenzae (strain Rd KW20) >gi|1573534
           (U32736) dimethyladenosine transferase (ksgA)
           [Haemophilus influenzae Rd]
           Length = 287
           
 Score = 32.6 bits (72), Expect = 0.94
 Identities = 15/45 (33%), Positives = 23/45 (50%), Gaps = 6/45 (13%)

Query: 1   PPPQVESSVVRLEI-KNPRPQVDYNEWDGLLRIV---FVRKNRTI 41
           P P+V+S+VVRL   K     V    W  L R+    F ++ +T+
Sbjct: 188 PAPKVDSAVVRLIPHKELPHPVKDLYW--LNRVCSQAFNQRRKTL 230


 sp|P37468|KSGA_BACSU DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLMETHIONINE-6-N',
           N'-ADENOSYL(RRNA) DIMETHYLTRANSFERASE) (16S RRNA
           DIMETHYLASE) (HIGH LEVEL KASUGAMYCIN RESISTANCE PROTEIN
           KSGA) (KASUGAMYCIN DIMETHYLTRANSFERASE)
           >gi|2126951|pir||S66071 high level kasgamycin resistance
           protein - Bacillus subtilis >gi|467431|gnl|PID|d1005819
           (D26185) high level kasgamycin resistance [Bacillus
           subtilis] >gi|2632309|gnl|PID|e1181975 (Z99104)
           dimethyladenosine transferase [Bacillus subtilis]
           Length = 292
           
 Score = 32.6 bits (72), Expect = 0.94
 Identities = 27/120 (22%), Positives = 45/120 (37%), Gaps = 30/120 (25%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNE---WDGLLRIVFVRKNRTISAGFKSTTVMDILEKN 57
           P P V+S+V+RL I    P VD      +  L++  F ++ +T+                
Sbjct: 197 PQPNVDSAVIRL-ILRDGPAVDVENESFFFQLIKASFAQRRKTLLNNL------------ 243

Query: 58  YKTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAG-KCDQNDFLRLLYAFHQV 116
                 +NN              K  I+ VL+ET++  KR G      +F  L    ++ 
Sbjct: 244 ------VNNLP-------EGKAQKSTIEQVLEETNIDGKRRGESLSIEEFAALSNGLYKA 290


 sp|Q03986|ERMD_BACLI RRNA ADENINE N-6-METHYLTRANSFERASE
           (MACROLIDE-LINCOSAMIDE-STREPTOGRAMIN B RESISTANCE
           PROTEIN) (ERYTHROMYCIN RESISTANCE PROTEIN)
           >gi|2126802|pir||I39885
           macrolide-lincosamide-streptogramin B resistance
           methylase - Bacillus licheniformis >gi|143201 (M29832)
           macrolide-lincosamide-streptogramin B resistance
           methylase [Bacillus licheniformis]
           Length = 287
           
 Score = 31.8 bits (70), Expect = 1.6
 Identities = 25/104 (24%), Positives = 43/104 (41%), Gaps = 8/104 (7%)

Query: 1   PPPQVESSVVRLEIKNPRP--QVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNY 58
           PPP+V+S++VR+  K   P     Y  + GL        N  +    +       + K+ 
Sbjct: 182 PPPKVDSAMVRITRKKDAPLSHKHYIAFRGLAEYALKEPNIPLCVALRGIFTPRQM-KHL 240

Query: 59  KTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDK--RAGK 100
           +  L +NNE    T G++ +     I   + +  +  K  RA K
Sbjct: 241 RKSLKINNE---KTVGTLTENQWAVIFNTMTQYVMHHKWPRANK 281


 sp|P51659|DHB4_HUMAN ESTRADIOL 17 BETA-DEHYDROGENASE 4 (17-BETA-HSD 4)
           (17-BETA-HYDROXYSTEROID DEHYDROGENASE 4)
           >gi|2134658|pir||S59136 17beta-hydroxysteroid
           dehydrogenase IV - human >gi|1050517 (X87176)
           17beta-hydroxysteroid dehydrogenase [Homo sapiens]
           >gi|4165049 (AF057740) 17-beta-hydroxysteroid
           dehydrogenase IV [Homo sapiens]
           Length = 736
           
 Score = 31.8 bits (70), Expect = 1.6
 Identities = 9/69 (13%), Positives = 22/69 (31%), Gaps = 3/69 (4%)

Query: 35  VRKNRTISAGFKSTT---VMDILEKNYKTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKET 91
           +R  RT+ A  +        + ++ N+K      N     +       + E +  +  E 
Sbjct: 247 LRWERTLGAIVRQKNHPMTPEAVKANWKKICDFENASKPQSIQESTGSIIEVLSKIDSEG 306

Query: 92  DLGDKRAGK 100
            +      +
Sbjct: 307 GVSANHTSR 315


 sp|P06992|KSGA_ECOLI DIMETHYLADENOSINE TRANSFERASE (S-ADENOSYLMETHIONINE-6-N',
           N'-ADENOSYL(RRNA) DIMETHYLTRANSFERASE) (16S RRNA
           DIMETHYLASE) (HIGH LEVEL KASUGAMYCIN RESISTANCE PROTEIN
           KSGA) (KASUGAMYCIN DIMETHYLTRANSFERASE)
           >gi|66453|pir||XYECRO rRNA
           (adenine-N6,N6-)-dimethyltransferase (EC 2.1.1.-) -
           Escherichia coli >gi|40917 (X04711) put. 16S rRNA
           methyltransferase (ksgA) (AA 1-273) [Escherichia coli]
           >gi|146571 (M11054) methyltransferase (ksgA)
           [Escherichia coli] >gi|216476|gnl|PID|d1001799 (D10483)
           rRNA (adenosine-N6,N6-)-dimethyltransferase [Escherichia
           coli] >gi|1786236 (AE000115)
           S-adenosylmethionine-6-N',N'-adenosyl (rRNA)
           dimethyltransferase [Escherichia coli]
           >gi|224747|prf||1112172A methyltransferase [Escherichia
           coli]
           Length = 273
           
 Score = 31.4 bits (69), Expect = 2.1
 Identities = 15/43 (34%), Positives = 22/43 (50%), Gaps = 2/43 (4%)

Query: 1   PPPQVESSVVRLEI--KNPRPQVDYNEWDGLLRIVFVRKNRTI 41
           PPP+V+S+VVRL      P P  D      +    F ++ +TI
Sbjct: 183 PPPKVDSAVVRLVPHATMPHPVKDVRVLSRITTEAFNQRRKTI 225


 gnl|PID|d1017401 (D90899) dimethyladenosine transferase [Synechocystis sp.]
           Length = 284
           
 Score = 31.4 bits (69), Expect = 2.1
 Identities = 14/47 (29%), Positives = 23/47 (48%), Gaps = 1/47 (2%)

Query: 1   PPPQVESSVVRLEI-KNPRPQVDYNEWDGLLRIVFVRKNRTISAGFK 46
           PPP+V+S+V+RL      +   D    D LL + F  + + +    K
Sbjct: 193 PPPKVDSAVIRLTPYPVEQLPGDRRLLDQLLCLGFANRRKMLRNNLK 239


 sp|Q04720|ERMJ_BACAN RRNA ADENINE N-6-METHYLTRANSFERASE
           (MACROLIDE-LINCOSAMIDE-STREPTOGRAMIN B RESISTANCE
           PROTEIN) (ERYTHROMYCIN RESISTANCE PROTEIN) >gi|143198
           (L08389) rRNA methyltransferase [Bacillus anthracis]
           Length = 287
           
 Score = 30.7 bits (67), Expect = 3.6
 Identities = 25/104 (24%), Positives = 43/104 (41%), Gaps = 8/104 (7%)

Query: 1   PPPQVESSVVRLEIKNPRP--QVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNY 58
           PPP+V+S++VR+  K   P     Y  + GL        N  +    +       + K+ 
Sbjct: 182 PPPKVDSAMVRITRKKDAPLSHKHYIAFRGLAEYALKEPNIPLCVRLRGIFTPRQM-KHL 240

Query: 59  KTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDK--RAGK 100
           +  L +NNE    T G++ +     I   + +  +  K  RA K
Sbjct: 241 RKSLKINNE---KTVGTLTENQWAVIFNTMTQYVMHHKWPRANK 281


 pir||A47697 rRNA methyltransferase - Bacillus anthracis
           Length = 287
           
 Score = 30.7 bits (67), Expect = 3.6
 Identities = 25/104 (24%), Positives = 43/104 (41%), Gaps = 8/104 (7%)

Query: 1   PPPQVESSVVRLEIKNPRP--QVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNY 58
           PPP+V+S++VR+  K   P     Y  + GL        N  +    +       + K+ 
Sbjct: 182 PPPKVDSAMVRITRKKDAPLSHKHYIAFRGLAEYALKEPNIPLCVRLRGIFTPRQM-KHL 240

Query: 59  KTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDK--RAGK 100
           +  L +NNE    T G++ +     I   + +  +  K  RA K
Sbjct: 241 RKSLKINNE---KTVGTLTENQWAVIFNTMTQYVMHHKWPRANK 281


 gi|987265 (U31089) Abl binding protein 3 [Homo sapiens]
           Length = 390
           
 Score = 30.7 bits (67), Expect = 3.6
 Identities = 11/30 (36%), Positives = 14/30 (46%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLL 30
           PPP VE  V       P P  DY E + ++
Sbjct: 281 PPPPVEEPVFDESPPPPPPPEDYEEEEAVV 310


 gi|2984099 (AE000757) dimethyladenosine transferase [Aquifex aeolicus]
           Length = 248
           
 Score = 30.7 bits (67), Expect = 3.6
 Identities = 11/41 (26%), Positives = 21/41 (50%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTI 41
           PPP+V+S+V++L      P  D   +   L  +F  + + +
Sbjct: 168 PPPKVQSAVIKLVKNEKFPVKDLKNYKKFLTKIFQNRRKVL 208


 gi|915311 (U23435) Abl interactor 2 [Homo sapiens]
           Length = 401
           
 Score = 30.3 bits (66), Expect = 4.8
 Identities = 11/30 (36%), Positives = 13/30 (42%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLL 30
           PPP VE  V       P P  DY E +  +
Sbjct: 292 PPPPVEEPVFDESPPPPPPPEDYEEEEAAV 321


 gnl|PID|e224171 (X95632) Arg protein tyrosine kinase-binding protein [Homo sapiens]
           Length = 475
           
 Score = 30.3 bits (66), Expect = 4.8
 Identities = 11/30 (36%), Positives = 13/30 (42%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLL 30
           PPP VE  V       P P  DY E +  +
Sbjct: 366 PPPPVEEPVFDESPPPPPPPEDYEEEEAAV 395


 gnl|PID|e224172 (X95677) argBPIB gene product [Homo sapiens]
           Length = 299
           
 Score = 30.3 bits (66), Expect = 4.8
 Identities = 11/30 (36%), Positives = 13/30 (42%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLL 30
           PPP VE  V       P P  DY E +  +
Sbjct: 190 PPPPVEEPVFDESPPPPPPPEDYEEEEAAV 219


 gnl|PID|e307817 (Y11786) dimethyladenosine transferase [Rickettsia prowazekii]
           >gi|3861209|gnl|PID|e1342953 (AJ235272)
           DIMETHYLADENOSINE TRANSFERASE (ksgA) [Rickettsia
           prowazekii]
           Length = 268
           
 Score = 30.3 bits (66), Expect = 4.8
 Identities = 11/47 (23%), Positives = 26/47 (54%), Gaps = 1/47 (2%)

Query: 1   PPPQVESSVVRLEI-KNPRPQVDYNEWDGLLRIVFVRKNRTISAGFK 46
           P P+V S++V++   +NP      N+ + + ++VF  + + I +  +
Sbjct: 186 PHPKVYSAIVKIIPLENPPSIALINKVEQITKLVFAGRRKMIKSSLR 232


 gi|4155935 (AE001555) DIMETHYLADENOSINE TRANSFERASE [Helicobacter pylori J99]
           Length = 271
           
 Score = 30.3 bits (66), Expect = 4.8
 Identities = 24/106 (22%), Positives = 38/106 (35%), Gaps = 19/106 (17%)

Query: 1   PPPQVESSVVRLEIKNPRP-----------------QVDYNEWDGLLRIVFVRKNRTISA 43
           PPP+V SSV  + IK P                   Q  +   +  L+  F    +T+S 
Sbjct: 167 PPPKVFSSVFEV-IKEPLKEKALASLAQAPFFEEALQKGFETLEDFLKACFSSPRKTLSN 225

Query: 44  GFKSTTVMDILEKNYKTFLAMNNEMVDDTKGSMHDVVKEKIDTVLK 89
             K +            FLA+ N+     + S      + + +VLK
Sbjct: 226 NLKKSVSYKEKLDKVLDFLALENQP-TSVRASEVKDYLKLLKSVLK 270


 sp|P45438|ERMK_BACLI RRNA ADENINE N-6-METHYLTRANSFERASE
           (MACROLIDE-LINCOSAMIDE-STREPTOGRAMIN B RESISTANCE
           PROTEIN) (ERYTHROMYCIN RESISTANCE PROTEIN)
           >gi|80107|pir||B42473 rRNA methylase ermK - Bacillus
           licheniformis >gi|143195 (M77505) erythromycin-inducible
           23S ribosomal RNA methylase [Bacillus licheniformis]
           Length = 287
           
 Score = 29.9 bits (65), Expect = 6.2
 Identities = 9/19 (47%), Positives = 14/19 (73%)

Query: 1   PPPQVESSVVRLEIKNPRP 19
           PPP+V+S++VR+  K   P
Sbjct: 182 PPPKVDSAMVRITRKKDAP 200


 pir||B48474 glycoprotein B - feline herpesvirus 1 (fragment)
           >gi|435794|bbs|138807 (S66371) glycoprotein B, gB
           [feline herpesvirus-1 FHV-1, C-27, Peptide Partial, 943
           aa] [Feline herpesvirus 1]
           Length = 943
           
 Score = 29.9 bits (65), Expect = 6.2
 Identities = 31/116 (26%), Positives = 45/116 (38%), Gaps = 12/116 (10%)

Query: 2   PPQVESSVVRLEIKN--PRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYK 59
           PP   S+VVRLE     P  ++  N  +G   I  + K       FK+    +I  KN  
Sbjct: 128 PPPSGSTVVRLEPPRACPDYKLGKNFTEG---IAVIFKENIAPYKFKA----NIYYKNII 180

Query: 60  TFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQ 115
                +      T     D V  K+  +   TDL D+R     + D++R  Y F  
Sbjct: 181 MTTVWSGSSYAVTTNRYTDRVPVKVQEI---TDLIDRRGMCLSKADYVRNNYQFTA 233


 pir||A56602 glycoprotein B homolog precursor - feline herpesvirus 1
           >gi|261095|bbs|120003 (S49775) glycoprotein B homolog
           [feline herpesvirus type 1 FHV-1, Peptide, 948 aa]
           [Feline herpesvirus 1]
           Length = 948
           
 Score = 29.9 bits (65), Expect = 6.2
 Identities = 31/116 (26%), Positives = 45/116 (38%), Gaps = 12/116 (10%)

Query: 2   PPQVESSVVRLEIKN--PRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYK 59
           PP   S+VVRLE     P  ++  N  +G   I  + K       FK+    +I  KN  
Sbjct: 129 PPPSGSTVVRLEPPRACPDYKLGKNFTEG---IAVIFKENIAPYKFKA----NIYYKNII 181

Query: 60  TFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQ 115
                +      T     D V  K+  +   TDL D+R     + D++R  Y F  
Sbjct: 182 MTTVWSGSSYAVTTNRYTDRVPVKVQEI---TDLIDRRGMCLSKADYVRNNYQFTA 234


 gi|471723 (U01887) gB homolog [Meleagrid herpesvirus type 1]
          Length = 864
          
 Score = 29.9 bits (65), Expect = 6.2
 Identities = 12/34 (35%), Positives = 18/34 (52%), Gaps = 1/34 (2%)

Query: 2  PPQVESSVVRLE-IKNPRPQVDYNEWDGLLRIVF 34
          PP V S++VRLE  +     +   EW   + I+F
Sbjct: 50 PPPVGSTIVRLEPPRKCPEPLKSTEWGEGIAILF 83


 prf||1911192A glycoprotein B [Feline herpesvirus 1]
           Length = 948
           
 Score = 29.9 bits (65), Expect = 6.2
 Identities = 31/116 (26%), Positives = 45/116 (38%), Gaps = 12/116 (10%)

Query: 2   PPQVESSVVRLEIKN--PRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYK 59
           PP   S+VVRLE     P  ++  N  +G   I  + K       FK+    +I  KN  
Sbjct: 129 PPPSGSTVVRLEPPRACPDYKLGKNFTEG---IAVIFKENIAPYKFKA----NIYYKNII 181

Query: 60  TFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQ 115
                +      T     D V  K+  +   TDL D+R     + D++R  Y F  
Sbjct: 182 MTTVWSGSSYAVTTNRYTDRVPVKVQEI---TDLIDRRGMCLSKADYVRNNYQFTA 234


 gi|2688507 (AE001160) dimethyladenosine transferase (ksgA) [Borrelia
           burgdorferi]
           Length = 281
           
 Score = 29.5 bits (64), Expect = 8.2
 Identities = 11/41 (26%), Positives = 25/41 (60%), Gaps = 1/41 (2%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTI 41
           P P+V+S+ ++L  K    + ++ E++ L+R VF  + + +
Sbjct: 196 PAPKVKSTTLKLIPKKNNIK-NFKEFNKLVRTVFSNRRKKL 235


 gnl|PID|e293782 (Y10669) lef-8 [Spodoptera littoralis nuclear polyhedrosis virus]
           Length = 909
           
 Score = 29.5 bits (64), Expect = 8.2
 Identities = 25/117 (21%), Positives = 46/117 (38%), Gaps = 13/117 (11%)

Query: 4   QVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKSTTVMDILEKNYKTFLA 63
           +V + VVR+   N R  +   E+D      FV   + + + F      ++     K FL 
Sbjct: 346 RVGNEVVRIINDNMRQDMLKQEYD------FV---KFVDSFFHG----EMTVAGKKFFLC 392

Query: 64  MNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQNDFLRLLYAFHQVGIHF 120
            +  + +    S+    K+ I++ L      D+     D +D   L+ AF+     F
Sbjct: 393 RDVRLPNVDYESIAKKFKDLIESNLIVRRDEDENGSGGDDDDDEWLMIAFNNRPTTF 449


 gi|471720 (U01886) gB homolog [Gallid herpesvirus 2]
           Length = 865
           
 Score = 29.5 bits (64), Expect = 8.2
 Identities = 26/108 (24%), Positives = 38/108 (35%), Gaps = 12/108 (11%)

Query: 2   PPQVESSVVRLE-IKNPRPQVDYNEWDGLLRIVFVRKNRTISAGFKST----TVMDILEK 56
           PP V ++VVRLE  +     +   EW   + I+F  K       FK T     V+     
Sbjct: 51  PPAVGTTVVRLEPPRKCPEPLKATEWGEGIAILF--KENINPYKFKVTLYYKNVIQTTTW 108

Query: 57  NYKTFLAMNNEMVDDTKGSMHDVVKEKIDTVLKETDLGDKRAGKCDQN 104
              T+  + N   D        V  E+I  V+        +A     N
Sbjct: 109 TGTTYRQITNRFTDR-----TPVSIEEITDVIDAKGRCSSKARYLRNN 151


 gi|1369852 (L46702) kinesin heavy chain-like protein [Solanum tuberosum]
           Length = 1265
           
 Score = 29.5 bits (64), Expect = 8.2
 Identities = 22/117 (18%), Positives = 45/117 (37%), Gaps = 8/117 (6%)

Query: 1   PPPQVESSVVRLEIKNPRPQVDYNEWDGLLRIVFVRKNRTISAGF--KSTTVMDILEK-N 57
           P P   SSV+  E+    P +D  + +G L+ +  + +     GF  K +    + EK  
Sbjct: 49  PTPNTLSSVLSPELAGAIPYIDRFQVEGFLKAMQKQIHSASKRGFFLKKSVGPQVREKFT 108

Query: 58  YKTFLAMNNEMVDDT----KGSMHDVVKEKIDTVLKETDLGD-KRAGKCDQNDFLRL 109
           ++  L    E +  +     G +     +   ++LK   +    R      ++ + L
Sbjct: 109 FEDMLCFQREPIPTSILKLNGDLISRAVKLFQSILKYMGIDSYDRVAPISLDERIEL 165


 gi|2314603 (AE000643) 16S rRNA (adenosine-N6,N6-)-dimethyltransferase (ksgA)
           [Helicobacter pylori 26695]
           Length = 271
           
 Score = 29.5 bits (64), Expect = 8.2
 Identities = 22/98 (22%), Positives = 34/98 (34%), Gaps = 18/98 (18%)

Query: 1   PPPQVESSVVRLEIKNPRP-----------------QVDYNEWDGLLRIVFVRKNRTISA 43
           PPP+V SSV  + IK P                   Q  +   +  L+  F    +T+S 
Sbjct: 167 PPPKVFSSVFEV-IKEPLKEKALASLAQAPFFEEALQKGFEMLEDFLKACFSSPRKTLSN 225

Query: 44  GFKSTTVMDILEKNYKTFLAMNNEMVDDTKGSMHDVVK 81
             K +            FLA+ N+        + D +K
Sbjct: 226 NLKKSVSYREKLDKVLDFLALENQPTSVRASEIKDYLK 263


CPU time:    14.07 user secs.	    0.57 sys. secs	   14.64 total secs.

  Database: Non-redundant GenBank CDS
  translations+PDB+SwissProt+SPupdate+PIR
    Posted date:  Jan 22, 1999 12:46 AM
  Number of letters in database: 107,844,596
  Number of sequences in database:  353,247
  
Lambda     K      H
   0.317    0.120    0.328 

Gapped
Lambda     K      H
   0.270   0.0421    0.230 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 17273284
Number of Sequences: 353247
Number of extensions: 481401
Number of successful extensions: 2224
Number of sequences better than 10.0: 43
Number of HSP's better than 10.0 without gapping: 27
Number of HSP's successfully gapped in prelim test: 16
Number of HSP's that attempted gapping in prelim test: 2172
Number of HSP's gapped (non-prelim): 45
length of query: 120
length of database: 107844596
effective HSP length: 51
effective length of query: 69
effective length of database: 89828999
effective search space: 6198200931
effective search space used: 6198200931
T: 11
A: 40
X1: 16 ( 7.3 bits)
X2: 38 (14.8 bits)
X3: 64 (24.9 bits)
S1: 41 (21.8 bits)
S2: 64 (29.5 bits)