TBLASTN 2.2.13 [Nov-27-2005]


Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs",  Nucleic Acids Res. 25:3389-3402.

Query= SPP00000033_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Homo sapiens # Complete
         (448 letters)

Database: C.merolae/genome.fa 
           20 sequences; 16,546,747 total letters

Searching....................done

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

dbj|AP006498.2| Cyanidioschyzon merolae DNA, chromosome 16, comp...   153   2e-37
dbj|AP006500.2| Cyanidioschyzon merolae DNA, chromosome 18, comp...    28   9.7  

>dbj|AP006498.2| Cyanidioschyzon merolae DNA, chromosome 16, complete genome, complete
              sequence
          Length = 908485

 Score =  153 bits (386), Expect = 2e-37
 Identities = 114/328 (34%), Positives = 155/328 (47%), Gaps = 27/328 (8%)
 Frame = -2

Query: 49     WRLTGFSGMKGUGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGP 108
              +RLT +S +KG GCKVPQ  LL+ +   T+  +R P+        +  S+  GL + +  
Sbjct: 479109 FRLTRYSTLKGCGCKVPQSELLEWMTA-TKNLLRSPM--------KSFSRPGGLMSTSNT 478957

Query: 109    SPTFPA-----LGIGMDSCVIPLRHGG--------------------LSLVQTTDFFYPL 143
                + P      LG+GMD+ V+ L                        L L+ TTDFF+PL
Sbjct: 478956 ESSTPGDVSVELGVGMDASVVRLTQLSSRQQSTEKDAENAPSFGTQPLYLISTTDFFFPL 478777

Query: 144    VEDPYMMGRIACANVLSDLYAMGITECDNXXXXXXXXXXXXXXX--REKVTPLMVKGFRD 201
              VEDPY+MGRIA ANVLSD YA+GI    +                 R  VT  M+ GF D
Sbjct: 478776 VEDPYLMGRIAAANVLSDAYALGIRTPPSTMLMLVCAATEIEPPALRAGVTQTMMAGFID 478597

Query: 202    XXXXXXXXXXXXXXXXNPWIIIGGVATVVCQPNEFIMPDSAVVGDVLVLTKPLGTQVAVN 261
                              NPW +IGGVA VV    EF+MPD+A  GD +VLTKPLGTQVA+N
Sbjct: 478596 ACHEAGTQVTGGQTTLNPWPLIGGVAEVVLPRGEFLMPDAAQPGDTIVLTKPLGTQVAIN 478417

Query: 262    AHQWLDNPERWNKVKMVVSREEVELAYQEAMFNMATLNRTAAGLMHTFNAHAATDITGFG 321
               +QWL   E  + + ++   E +E   Q+      T        +       A DI G G
Sbjct: 478416 LYQWLGQMEDAHYLHLL---ESLETTQQQEYGPEDTRRSVHGQFVPPSEETHARDIPGPG 478246

Query: 322    ILGHSQNLAKQQRNEVSFVIHNLPIIAK 349
                 H+ + A+  R   S ++  LP  A+
Sbjct: 478245 -PAHTTHDARLSRALSSSILEQLPPAAR 478165



 Score =  112 bits (280), Expect = 3e-25
 Identities = 72/166 (43%), Positives = 97/166 (58%), Gaps = 12/166 (7%)
 Frame = -2

Query: 272    WNKVKMVVSREEVELAYQEAMFNMATLNRTAAGLM-HTFNAHAATDITGFGILGHSQNLA 330
              ++K++ V++  ++  AY  A  +MA LN  AA L+  +  AHA TDITGFGILGH+QNLA
Sbjct: 478128 YSKIEQVLTPADICRAYDVACLSMAHLNARAAQLLLESGAAHAVTDITGFGILGHAQNLA 477949

Query: 331    KQQRNEVSFVIHNLPIIAKMAAVSKA---SGRFGLLQGTSAETSGGLLICLPREQ----- 382
              + Q+  V   I  LPII  M +  +A      F LL+G SAETSGGLLIC          
Sbjct: 477948 RNQKAPVDMEIDVLPIIRGMTSADRALQGGKMFRLLEGFSAETSGGLLICFGASSTERSS 477769

Query: 383    ---AARFCSEIKSSKYGEGHQAWIVGIVEKGNRTARIIDKPRVIEV 425
                 A +FC   K +   +   AWIVG V  G+R ARI+D  +++EV
Sbjct: 477768 QVVARQFCE--KLATVDDHWPAWIVGRVVPGSREARILDHVKILEV 477637


>dbj|AP006500.2| Cyanidioschyzon merolae DNA, chromosome 18, complete genome, complete
              sequence
          Length = 1253087

 Score = 27.7 bits (60), Expect = 9.7
 Identities = 17/42 (40%), Positives = 23/42 (54%)
 Frame = +3

Query: 242    AVVGDVLVLTKPLGTQVAVNAHQWLDNPERWNKVKMVVSREE 283
              A++    +LT  +G  + VNAH WLD   RW+  K    REE
Sbjct: 982623 ALLSPCCILTGAVG--LIVNAHSWLD--ARWSDWKDSYLREE 982736


  Database: C.merolae/genome.fa
    Posted date:  Feb 9, 2011  9:02 AM
  Number of letters in database: 16,546,747
  Number of sequences in database:  20
  
Lambda     K      H
   0.320    0.136    0.412 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 7,664,744
Number of Sequences: 20
Number of extensions: 113487
Number of successful extensions: 593
Number of sequences better than 10.0: 2
Number of HSP's better than 10.0 without gapping: 1
Number of HSP's successfully gapped in prelim test: 1
Number of HSP's that attempted gapping in prelim test: 577
Number of HSP's gapped (non-prelim): 33
length of query: 448
length of database: 5,515,582
effective HSP length: 99
effective length of query: 349
effective length of database: 5,513,602
effective search space: 1924247098
effective search space used: 1924247098
frameshift window, decay const: 40,  0.1
T: 13
A: 40
X1: 16 ( 7.4 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.8 bits)
S2: 60 (27.7 bits)