RPS-BLAST 2.2.6 [Apr-09-2003]

Database: All 
           18,039 sequences; 5,506,404 total letters



Query= 2166 (315 letters)

Distribution of 8 Blast Hits on the Query Sequence




                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gnl|CDD|10620  COG0752, GlyQ, Glycyl-tRNA synthetase, alpha ...   568   e-163 
gnl|CDD|6639  pfam02091, tRNA_synt_2e, Glycyl-tRNA synthetas...   526   e-150 
gnl|CDD|9889  COG0013, AlaS, Alanyl-tRNA synthetase [Transla...    32   0.20  
gnl|CDD|14302  COG5178, PRP8, U5 snRNP spliceosome subunit [...    29   1.2   
gnl|CDD|20870  KOG3085, KOG3085, Predicted hydrolase (HAD su...    29   1.3   
gnl|CDD|25824  pfam01808, AICARFT_IMPCHas, AICARFT/IMPCHase ...    27   5.3   
gnl|CDD|19581  KOG1795, KOG1795, U5 snRNP spliceosome subuni...    27   6.2   
gnl|CDD|10013  COG0138, PurH, AICAR transformylase/IMP cyclo...    27   8.4   

>gnl|CDD|10620 COG0752, GlyQ, Glycyl-tRNA synthetase, alpha subunit [Translation, ribosomal structure and biogenesis] Length = 298 Score = 568 bits (1465), Expect = e-163 Identities = 235/298 (78%), Positives = 261/298 (87%)
Query:  5    TPAVRTFQDLILALQNYWAEQGCVVLQPYDMEVGAGTFHTATFLRAIGPETWNAAYVQPS  64
Sbjct:  1    MNKKLTFQGLILTLQNYWAEQGCTILQPYDMEVGAGTFHPATFLRALGPEPWNAAYVQPS  60

Query:  65   RRPTDGRYGENPNRLQHYYQFQVVLKPNPENFQELYLGSLKAIGIDPLVHDIRFVEDNWE  124
Sbjct:  61   RRPTDGRYGENPNRLQHYYQFQVIIKPSPDNIQELYLGSLEALGIDPLEHDIRFVEDNWE  120

Query:  125  SPTLGAWGLGWEIWLNGMEVTQFTYFQQVGGIECYPVTGEITYGLERLAMYLQGVDSVYD  184
Sbjct:  121  NPTLGAWGLGWEVWLDGMEVTQFTYFQQVGGLECKPVSGEITYGLERLAMYIQGVDNVYD  180

Query:  185  LVWTDGPFGKVTYGDVFHQNEVEQSTFNFEHANVPKLFELFDFYESEANRLIALELPLPT  244
Sbjct:  181  LEWNDGPGGKVTYGDVFLQNEVEQSKYNFEYADVDMLFRHFDDYEKEAKRLLELGLVLPA  240

Query:  245  YEMVLKASHTFNLLDARRAISVTERQRYILRVRTLARAVAQSYLQARARLGFPMATPE  302
Sbjct:  241  YDYVLKASHTFNLLDARGAISVTERQRYILRIRNLARAVAEAYLESREALGFPLLKKE  298


>gnl|CDD|6639 pfam02091, tRNA_synt_2e, Glycyl-tRNA synthetase alpha subunit Length = 285 Score = 526 bits (1357), Expect = e-150 Identities = 232/288 (80%), Positives = 257/288 (89%), Gaps = 3/288 (1%)
Query:  11   FQDLILALQNYWAEQGCVVLQPYDMEVGAGTFHTATFLRAIGPETWNAAYVQPSRRPTDG  70
Sbjct:  1    FQQMILTLQEYWASQGCVVMQPYDMEVGAGTFHPATFLRSLGPEPWNVAYVEPSRRPTDG  60

Query:  71   RYGENPNRLQHYYQFQVVLKPNPENFQELYLGSLKAIGIDPLVHDIRFVEDNWESPTLGA  130
Sbjct:  61   RYGENPNRLQHYYQFQVVLKPSPDNIQELYLGSLKALGIDPLDHDIRFVEDNWESPTLGA  120

Query:  131  WGLGWEIWLNGMEVTQFTYFQQVGGIECYPVTGEITYGLERLAMYLQGVDSVYDLVWTDG  190
Sbjct:  121  WGLGWEVWLDGMEITQFTYFQQVGGLECKPVSGEITYGLERLAMYLQGVDNVYDLVWADG  180

Query:  191  PFGKVTYGDVFHQNEVEQSTFNFEHANVPKLFELFDFYESEANRLIALELPLPTYEMVLK  250
Sbjct:  181  P---VTYGDIFLQNEVEQSTYNFETANVDRLFKLFDFYEKEALQLLENGLPLPAYDFVLK  237

Query:  251  ASHTFNLLDARRAISVTERQRYILRVRTLARAVAQSYLQARARLGFPM  298
Sbjct:  238  ASHAFNLLDARGAISVTERQRYILRIRTLARAVAELYYAQREKLGFPL  285


>gnl|CDD|9889 COG0013, AlaS, Alanyl-tRNA synthetase [Translation, ribosomal structure and biogenesis] Length = 879 Score = 31.8 bits (72), Expect = 0.20 Identities = 13/23 (56%), Positives = 15/23 (65%), Gaps = 3/23 (13%)
Query:  162  TGEITYGLERLAMYLQGVDSVYD  184
Sbjct:  227  TG---MGLERIAAVLQGVPTNYD  246


>gnl|CDD|14302 COG5178, PRP8, U5 snRNP spliceosome subunit [RNA processing and modification] Length = 2365 Score = 29.3 bits (65), Expect = 1.2 Identities = 19/71 (26%), Positives = 26/71 (36%), Gaps = 2/71 (2%)
Query:  115   DIRFVEDNWESPTLGAWGLGWEIWLNGMEVTQFTYFQQVGGIECYPVTGEITYGLERLAM  174
Sbjct:  643   YYAFNEGVGKGPGCGFWGPQWRVWLFFLRGHIPLLERYIGNLVTRQFEGRSDYNPKPLTK  702

Query:  175   YLQGVDSVYDL  185
Sbjct:  703   --QRSDSGYDL  711


>gnl|CDD|20870 KOG3085, KOG3085, Predicted hydrolase (HAD superfamily) [General function prediction only] Length = 237 Score = 29.1 bits (65), Expect = 1.3 Identities = 19/68 (27%), Positives = 29/68 (42%), Gaps = 6/68 (8%)
Query:  90   KPNPENFQELYLGSLKAIGIDPLVHDIRFVEDNWESPTLGAWGLGWEIWLNGMEVTQFTY  149
Sbjct:  168  KPDPRIFQLA----LERLGVKP--EECVHIGDLLENDYEGARNLGWHAILVDNSITALKE  221

Query:  150  FQQVGGIE  157
Sbjct:  222  LEYKLGID  229


>gnl|CDD|25824 pfam01808, AICARFT_IMPCHas, AICARFT/IMPCHase bienzyme. This is a family of bifunctional enzymes catalysing the last two steps in de novo purine biosynthesis. The bifunctional enzyme is found in both prokaryotes and eukaryotes. The second last step is catalysed by 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase EC:2.1.2.3 (AICARFT), this enzyme catalyses the formylation of AICAR with 10-formyl-tetrahydrofolate to yield FAICAR and tetrahydrofolate. The last step is catalysed by IMP (Inosine monophosphate) cyclohydrolase EC:3.5.4.10 (IMPCHase), cyclizing FAICAR (5-formylaminoimidazole-4-carboxamide ribonucleotide) to IMP Length = 315 Score = 27.2 bits (60), Expect = 5.3 Identities = 16/65 (24%), Positives = 27/65 (41%), Gaps = 8/65 (12%)
Query:  43   HTATFLRAI--------GPETWNAAYVQPSRRPTDGRYGENPNRLQHYYQFQVVLKPNPE  94
Sbjct:  45   HTAAYDAAIADYFRKQYSPKHPERTFTLPGVRKLSLRYGENPHQQAAFYVDPTLKEGSIA  104

Query:  95   NFQEL  99
Sbjct:  105  TAEQL  109


>gnl|CDD|19581 KOG1795, KOG1795, U5 snRNP spliceosome subunit [RNA processing and modification] Length = 2321 Score = 26.9 bits (59), Expect = 6.2 Identities = 9/26 (34%), Positives = 11/26 (42%)
Query:  115   DIRFVEDNWESPTLGAWGLGWEIWLN  140
Sbjct:  609   YYRFNTGVGKGPGCGFWAPGWRVWLF  634


>gnl|CDD|10013 COG0138, PurH, AICAR transformylase/IMP cyclohydrolase PurH (only IMP cyclohydrolase domain in Aful) [Nucleotide transport and metabolism] Length = 515 Score = 26.7 bits (59), Expect = 8.4 Identities = 9/41 (21%), Positives = 17/41 (41%)
Query:  45   ATFLRAIGPETWNAAYVQPSRRPTDGRYGENPNRLQHYYQF  85
Sbjct:  187  ANYFSELLGGEFPEQLNLSFTKKQDLRYGENPHQSAAFYRD  227

Lambda     K      H
   0.321    0.138    0.423 

Gapped
Lambda     K      H
   0.267   0.0539    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 1100, Extension: 100
Number of Hits to DB: 5,577,383
Number of Sequences: 0
Number of extensions: 565671
Number of successful extensions: 211
Number of sequences better than 10.0: 1
Number of HSP's better than 10.0 without gapping: 1
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 0
Number of HSP's gapped (non-prelim): 210
length of query: 0
length of database: 5,506,404
effective HSP length: 94
effective length of query: -94
effective length of database: 3,810,738
effective search space: -358209372
effective search space used: 842173098
T: 11
A: 40
X1: 1600 (741.1 bits)
X2: 3800 (1463.8 bits)
X3: 6400 (2465.3 bits)
S1: 4100 (1902.0 bits)
S2: 58 (26.6 bits)