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= SPP00000056_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Anopheles gambiae # Complete
         (369 letters)

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

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_36                                                           129   1e-29

>scaffold_36
          Length = 507435

 Score =  129 bits (323), Expect = 1e-29
 Identities = 81/203 (39%), Positives = 109/203 (53%), Gaps = 1/203 (0%)
 Frame = -3

Query: 167    GGAASAVCHRSELIMPYNAQPGDALVLTKPLGTQLATNAYIWMGEQQSDSWAQLRERFTV 226
              GG   A  H+ + IMP  A+ GD LVLTKPLG QLA N Y WM    S  W  ++     
Sbjct: 409534 GGVGIASNHKDDFIMPNGAEAGDYLVLTKPLGVQLAVNTYQWMTTIPS-KWDAIKGITNE 409358

Query: 227    ADIEQTYRIALESMSRLNKTGAELMKKYGAHAATDVTGFGLYGHAENLASHQTADVDFHL 286
              ++I   +  A +SM  LN   A L KK+   A+TDVTGFG+ GHA  LA      V F +
Sbjct: 409357 SEIN*AFDHATKSMCTLNAAAASLFKKHNIKASTDVTGFGILGHARYLARA*RRKV*FVI 409178

Query: 287    DTLPIIKNVREIAETLGRGAKLLAGKAVETSGGLLICLPREAAAGFCEEYRRCTKHEAWI 346
              D LPII+N+ ++ +   R  + L G A ETSGGLLI    +    F ++Y+       W+
Sbjct: 409177 DKLPIIQNLYKL-DKKARDFRFLEGLAAETSGGLLIA--TKDPDSFLKDYKHLIGEWGWV 409007

Query: 347    VGHVEKG-ARGVKMNPNLNILSV 368
              +G VEK   R   ++ N+ I  V
Sbjct: 409006 IGRVEKSQERDAVISKNVKIQEV 408938



 Score = 94.4 bits (233), Expect(2) = 7e-21
 Identities = 55/162 (33%), Positives = 91/162 (56%), Gaps = 2/162 (1%)
 Frame = -2

Query: 4      PESYGLSHDFRLTKFSTLRGUGSKVPQDVLNRLLAGVYGEQLGDKDGKGPNKEEGVGIGL 63
              PE+YGL   F +T F+ L+G G KVP   L + L G  GE  GD   + P          
Sbjct: 410570 PENYGLPKGFAMTNFTELKGCGCKVP*T*LLKYL-G*IGE--GDIGKETP---------- 410430

Query: 64     DSSVIALKH--DLFLVQSVDFFYPLIDDPFMLGKIALANVVSDVFAVGATEIDQIKLIVT 121
              D SV  +K   D   + ++D+FYPL+ +P+  G+IA  NV+SD++A+G T++D + +++ 
Sbjct: 410429 DVSVYPVKRT*DFHTISTIDYFYPLVANPYA*GRIACCNVLSDLYAMGVTQVDTMLMVLG 410250

Query: 122    APTEFTEQEREVVVPMVMQGFLEAAKACNAPVQIGSIAENPW 163
                T+ T++E+ +   ++++GF + AK     V  G    +PW
Sbjct: 410249 VCTKMTDEEK*ISTSLMLEGFNDCAKEAKTSVTGGQSV*SPW 410124



 Score = 25.8 bits (55), Expect(2) = 7e-21
 Identities = 15/40 (37%), Positives = 19/40 (47%)
 Frame = -3

Query: 167    GGAASAVCHRSELIMPYNAQPGDALVLTKPLGTQLATNAY 206
              GG   A  H+   I+P  +     L LTKP   QL+ N Y
Sbjct: 410089 GGVGIASNHKDYFIIPNFSYSFYYLFLTKPFFFQLSFNTY 409970


  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.399 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 21,122,601
Number of Sequences: 697
Number of extensions: 220269
Number of successful extensions: 523
Number of sequences better than 10.0: 1
Number of HSP's better than 10.0 without gapping: 11
Number of HSP's successfully gapped in prelim test: 3
Number of HSP's that attempted gapping in prelim test: 500
Number of HSP's gapped (non-prelim): 36
length of query: 369
length of database: 23,960,378
effective HSP length: 107
effective length of query: 262
effective length of database: 23,885,799
effective search space: 6258079338
effective search space used: 6258079338
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: 64 (29.3 bits)