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= SPP00000048_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Drosophila melanogaster # Complete
         (370 letters)

Database: C.fasciculata/genome.fa 
           1089 sequences; 37,048,325 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

Cf_Contig765  |  | 1 to 162644                                        223   2e-58
Cf_Contig741  |  | 1 to 201246                                         32   1.1  
Cf_Contig972  |  | 1 to 178207                                         31   1.9  

>Cf_Contig765  |  | 1 to 162644
          Length = 162644

 Score =  223 bits (568), Expect = 2e-58
 Identities = 138/362 (38%), Positives = 202/362 (55%), Gaps = 12/362 (3%)
 Frame = +3

Query: 2      FQPEKHGLEPDFQLTKFTTHTGUSCKIPQKVLEKYLRGTEIENKNNDGYLIGSGMDCAVI 61
              F P   GL  +FQLT +T   G SCK+PQ  L   L+         D      GMDC+++
Sbjct: 101625 FDPVSLGLPAEFQLTDYTRLKGCSCKLPQPKLLALLQAVATSPGQKD-----VGMDCSIV 101789

Query: 62     PLK----RHKDYLLIQTVDFFYPMVNDPELLGRIALANVLSDVYAVGVTQFDTVEMIVST 117
              PL     + +   L+ T DFF+P V DP L G I  ANVLSD+Y++G+   DT+ M+++ 
Sbjct: 101790 PLHHTNGKGEALFLVSTTDFFFPSVADPFLQGHIGAANVLSDLYSMGIAHCDTMLMLLAA 101969

Query: 118    STSFSEKERDVVIGLVMKGFQNSLKANGYRNTPLIIRQLKINPWCIIGGIATSVCRSEEI 177
              ST   E ER V    +M+GF    +      T +   Q  +NPW +IGG+A +V    E+
Sbjct: 101970 STDMDEDERLVTTREIMRGFAERARLAA---TSVTGGQTVMNPWPLIGGVAMAVVSEAEM 102140

Query: 178    ILPSN-AQPGDVLVLTKPLGGQMAMDAHLWQLNQTEKYKKLL-SECSDADIKETFEIAVK 235
              + P+   + GDVLVLTKPLG Q+AM+   W +  +  Y+ ++    S  +I+E +  A  
Sbjct: 102141 VRPAGLLRAGDVLVLTKPLGCQIAMNLKQWLMRPSPLYEDVIEGHISRDEIEELYNKATD 102320

Query: 236    SMTYLNKNAALLMHKYQAHCATDITGFGLLGHANNLAQFQ-----KEKVLFQINKLPIIK 290
              SM  LN+  A LM K+ AH ATD+TGFG+LGHANN    Q      + +   +++LP++K
Sbjct: 102321 SMRRLNREGARLMLKHGAHGATDVTGFGILGHANNFGAAQAVGSAPQSLCLVLDRLPMLK 102500

Query: 291    NVLKFSTLVGQSTKFRSGRSVETSGGLLICLPA-DAADKFCRDFEEATNGEQKSFQIGHV 349
                +  S L+    +   G S ETSGGLL+ L + +AA+ F ++ EEA +G   S+ +G V
Sbjct: 102501 TAVAASRLMQDKYRLFEGYSAETSGGLLVALASREAAEAFQKELEEA-DGGCPSWIVGRV 102677

Query: 350    TA 351
               A
Sbjct: 102678 EA 102683


>Cf_Contig741  |  | 1 to 201246
          Length = 201246

 Score = 31.6 bits (70), Expect = 1.1
 Identities = 19/76 (25%), Positives = 34/76 (44%), Gaps = 9/76 (11%)
 Frame = +1

Query: 296    STLVGQSTKFRSGRSVET---------SGGLLICLPADAADKFCRDFEEATNGEQKSFQI 346
              + + G+++  R G+ ++T         S   L C P    DK+C     ++ GEQ     
Sbjct: 121363 TAIAGRTSAMRGGQKIDTKT*QGL*SGSVAFLACAPKALLDKYCFCGVRSSEGEQTRMVP 121542

Query: 347    GHVTAANESDAVLCED 362
                ++ +  S +VLC D
Sbjct: 121543 SVLSFSLLSSSVLCHD 121590


>Cf_Contig972  |  | 1 to 178207
          Length = 178207

 Score = 30.8 bits (68), Expect = 1.9
 Identities = 12/36 (33%), Positives = 20/36 (55%)
 Frame = +3

Query: 160    PWCIIGGIATSVCRSEEIILPSNAQPGDVLVLTKPL 195
              PWC++GG+   VC+   +  P  +    +L+ TK L
Sbjct: 166710 PWCVVGGLVECVCKGFVLKAPP*SSSEGLLLSTKNL 166817


  Database: C.fasciculata/genome.fa
    Posted date:  Nov 21, 2011  7:46 PM
  Number of letters in database: 37,048,325
  Number of sequences in database:  1089
  
Lambda     K      H
   0.319    0.135    0.394 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 13,919,428
Number of Sequences: 1089
Number of extensions: 186272
Number of successful extensions: 782
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 40
Number of HSP's successfully gapped in prelim test: 28
Number of HSP's that attempted gapping in prelim test: 709
Number of HSP's gapped (non-prelim): 169
length of query: 370
length of database: 12,349,441
effective HSP length: 103
effective length of query: 267
effective length of database: 12,237,274
effective search space: 3267352158
effective search space used: 3267352158
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: 62 (28.5 bits)