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= SPP00000108_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Mus musculus # Complete
         (452 letters)

Database: genome.fa 
           697 sequences; 71,881,135 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_36                                                           142   9e-34

>scaffold_36
          Length = 507435

 Score =  142 bits (359), Expect = 9e-34
 Identities = 82/199 (41%), Positives = 109/199 (54%), Gaps = 1/199 (0%)
 Frame = -3

Query: 225    IIGGVATVVCQQNEFIMPDSAVVGDVLVLTKPLGTQVAANAHQWLDN-PEKWNKIKMVVS 283
              I GGV      +++FIMP+ A  GD LVLTKPLG Q+A N +QW+   P KW+ IK + +
Sbjct: 409540 IRGGVGIASNHKDDFIMPNGAEAGDYLVLTKPLGVQLAVNTYQWMTTIPSKWDAIKGITN 409361

Query: 284    REEVELAYQEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQKNEVSFV 343
                E+  A+  A  +M TLN  AA L    N  A+TD+TGFGILGH++ LA+  + +V FV
Sbjct: 409360 ESEIN*AFDHATKSMCTLNAAAASLFKKHNIKASTDVTGFGILGHARYLARA*RRKV*FV 409181

Query: 344    IHNLPXXXXXXXXXXXXGRFGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQ 403
              I  LP              F  L+G +AETSGGLLI           S +K  K+  G  
Sbjct: 409180 IDKLPIIQNLYKLDKKARDFRFLEGLAAETSGGLLIATKDPD-----SFLKDYKHLIGEW 409016

Query: 404    AWIVGIVEKGNRTARIIDK 422
               W++G VEK      +I K
Sbjct: 409015 GWVIGRVEKSQERDAVISK 408959



 Score = 78.6 bits (192), Expect(2) = 5e-18
 Identities = 53/163 (32%), Positives = 74/163 (45%), Gaps = 1/163 (0%)
 Frame = -2

Query: 43     PQTLGFSPSWRLTSFSGMKGUGCKVPQETLLKLLEGLTRPALQPPLTSXXXXXXXXXXXX 102
              P+  G    + +T+F+ +KG GCKVP   LLK L  +    +                  
Sbjct: 410570 PENYGLPKGFAMTNFTELKGCGCKVP*T*LLKYLG*IGEGDI------------------ 410445

Query: 103    XXXSTRPGPGSAFPSLSIGMDSCVIPLRHG-GLSLVQTTDFFYPLVEDPYMMGRIACANV 161
                       G   P +S      V P++       + T D+FYPLV +PY  GRIAC NV
Sbjct: 410444 ---------GKETPDVS------VYPVKRT*DFHTISTIDYFYPLVANPYA*GRIACCNV 410310

Query: 162    LSDLYAMGITECDNXXXXXXXXXXXXEKEREKVTPLMIKGFRD 204
              LSDLYAMG+T+ D             ++E+   T LM++GF D
Sbjct: 410309 LSDLYAMGVTQVDTMLMVLGVCTKMTDEEK*ISTSLMLEGFND 410181



 Score = 32.3 bits (72), Expect(2) = 5e-18
 Identities = 20/71 (28%), Positives = 38/71 (53%), Gaps = 1/71 (1%)
 Frame = -3

Query: 225    IIGGVATVVCQQNEFIMPDSAVVGDVLVLTKPLGTQVAANAHQWLDN-PEKWNKIKMVVS 283
              I GGV      ++ FI+P+ +     L LTKP   Q++ N +Q++   P K+  IK + +
Sbjct: 410095 IRGGVGIASNHKDYFIIPNFSYSFYYLFLTKPFFFQLSFNTYQFITTIPSKFYSIKFITN 409916

Query: 284    REEVELAYQEA 294
                 + L++ ++
Sbjct: 409915 YI*IYLSFYQS 409883


  Database: genome.fa
    Posted date:  Feb 10, 2010 10:01 AM
  Number of letters in database: 71,881,135
  Number of sequences in database:  697
  
Lambda     K      H
   0.319    0.136    0.413 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 22,710,678
Number of Sequences: 697
Number of extensions: 235179
Number of successful extensions: 476
Number of sequences better than 10.0: 1
Number of HSP's better than 10.0 without gapping: 3
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 465
Number of HSP's gapped (non-prelim): 10
length of query: 452
length of database: 23,960,378
effective HSP length: 109
effective length of query: 343
effective length of database: 23,884,405
effective search space: 8192350915
effective search space used: 8192350915
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.7 bits)
S2: 65 (29.6 bits)