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: genome.fa 
           697 sequences; 71,881,135 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_36                                                           152   1e-36

>scaffold_36
          Length = 507435

 Score =  152 bits (384), Expect = 1e-36
 Identities = 85/199 (42%), Positives = 115/199 (57%), Gaps = 1/199 (0%)
 Frame = -3

Query: 222    IIGGVATVVCQPNEFIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDN-PERWNKVKMVVS 280
              I GGV       ++FIMP+ A  GD LVLTKPLG Q+AVN +QW+   P +W+ +K + +
Sbjct: 409540 IRGGVGIASNHKDDFIMPNGAEAGDYLVLTKPLGVQLAVNTYQWMTTIPSKWDAIKGITN 409361

Query: 281    REEVELAYQEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQRNEVSFV 340
                E+  A+  A  +M TLN  AA L    N  A+TD+TGFGILGH++ LA+  R +V FV
Sbjct: 409360 ESEIN*AFDHATKSMCTLNAAAASLFKKHNIKASTDVTGFGILGHARYLARA*RRKV*FV 409181

Query: 341    IHNLPIIAKMAAVSKASGRFGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQ 400
              I  LPII  +  + K +  F  L+G +AETSGGLLI           S +K  K+  G  
Sbjct: 409180 IDKLPIIQNLYKLDKKARDFRFLEGLAAETSGGLLIATKDPD-----SFLKDYKHLIGEW 409016

Query: 401    AWIVGIVEKGNRTARIIDK 419
               W++G VEK      +I K
Sbjct: 409015 GWVIGRVEKSQERDAVISK 408959



 Score = 78.6 bits (192), Expect(2) = 2e-17
 Identities = 56/175 (32%), Positives = 76/175 (43%), Gaps = 1/175 (0%)
 Frame = -2

Query: 28     LGRSFSNYRPFEPQALGLSPSWRLTGFSGMKGUGCKVPQEALLKLLAGLTRPDVRPPLGR 87
              LG  +     + P+  GL   + +T F+ +KG GCKVP   LLK L  +   D+      
Sbjct: 410606 LGTKYPKLFNW*PENYGLPKGFAMTNFTELKGCGCKVP*T*LLKYLG*IGEGDI------ 410445

Query: 88     GLVGGQEEASQEAGLPAGAGPSPTFPALGIGMDSCVIPLRHG-GLSLVQTTDFFYPLVED 146
                  G+E                         D  V P++       + T D+FYPLV +
Sbjct: 410444 ----GKETP-----------------------DVSVYPVKRT*DFHTISTIDYFYPLVAN 410346

Query: 147    PYMMGRIACANVLSDLYAMGITECDNXXXXXXXXXXXXXXXREKVTPLMVKGFRD 201
              PY  GRIAC NVLSDLYAMG+T+ D                +   T LM++GF D
Sbjct: 410345 PYA*GRIACCNVLSDLYAMGVTQVDTMLMVLGVCTKMTDEEK*ISTSLMLEGFND 410181



 Score = 30.4 bits (67), Expect(2) = 2e-17
 Identities = 18/71 (25%), Positives = 37/71 (52%), Gaps = 1/71 (1%)
 Frame = -3

Query: 222    IIGGVATVVCQPNEFIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDN-PERWNKVKMVVS 280
              I GGV       + FI+P+ +     L LTKP   Q++ N +Q++   P ++  +K + +
Sbjct: 410095 IRGGVGIASNHKDYFIIPNFSYSFYYLFLTKPFFFQLSFNTYQFITTIPSKFYSIKFITN 409916

Query: 281    REEVELAYQEA 291
                 + 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.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: 22,635,223
Number of Sequences: 697
Number of extensions: 238525
Number of successful extensions: 513
Number of sequences better than 10.0: 1
Number of HSP's better than 10.0 without gapping: 2
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 506
Number of HSP's gapped (non-prelim): 6
length of query: 448
length of database: 23,960,378
effective HSP length: 109
effective length of query: 339
effective length of database: 23,884,405
effective search space: 8096813295
effective search space used: 8096813295
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: 65 (29.6 bits)