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= SPP00000032_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Homo sapiens # Complete
         (392 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             82   1e-15
scaffold_15                                                            32   1.1  
scaffold_12                                                            31   2.4  

>scaffold_6
          Length = 2755913

 Score = 82.0 bits (201), Expect = 1e-15
 Identities = 76/288 (26%), Positives = 121/288 (42%), Gaps = 28/288 (9%)
 Frame = -1

Query: 82      LVQTTDYIYPIVDDPYMMGRIACANVLSDLYAMGVTECDNMLMLLGVSNKMTDRERDKVM 141
               +V T D+    V DPY+ G+IA  + LSD  AMG      + ++  V    ++R  +  +
Sbjct: 1819187 VVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALV--VVPYGSNRCVEDTL 1819014

Query: 142     PLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNE-FIMPDNAVPGDVLVLT 200
                 ++ G + A + AG ++ GG T     + LG     V      F+  +  V GD L+LT
Sbjct: 1819013 FQLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILT 1818834

Query: 201     KPLGTQVAVAVHQWLDIPEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTAAGLMHTFNA 260
               KPLGT    A             ++ +    +V    +E + +M + NR AA  +     
Sbjct: 1818833 KPLGTGAVFA-----------GDMRALSRAAEV----RECLASMTQSNRAAATELKEHGC 1818699

Query: 261     HAATDITGFGILGHAQNLAKQQRN-EVSFVIHNLPVLAKMAAVSKACGNMFGLM------ 313
                + TD+TGFG+LGH   + +      V   +   PV A   A   A   +F  +      
Sbjct: 1818698 TSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFA--GAKRLARDGVFSSLQDENLR 1818525

Query: 314     -------HGTC-------------PETSGGLLICLPREQAARFCAEIK 341
                      HG               P+T+GGLL  +P ++A    A I+
Sbjct: 1818524 AKRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIR 1818381



 Score = 29.3 bits (64), Expect = 9.3
 Identities = 17/58 (29%), Positives = 26/58 (44%)
 Frame = +1

Query: 34      PQDVLQKLLESLQENHFQEDEQFLGAVMPRLGIGMDTCVIPLRHGGLSLVQTTDYIYP 91
               PQ+VL++  +   E   +++E  LG   P L +     V+P R G   L    D   P
Sbjct: 2275972 PQEVLERQTQDEDERPDEQEEVALGDRAPGLALAPGVAVLPHRFGEGRLAHEGDDFLP 2276145


>scaffold_15
          Length = 987846

 Score = 32.3 bits (72), Expect = 1.1
 Identities = 19/65 (29%), Positives = 32/65 (49%), Gaps = 5/65 (7%)
 Frame = +2

Query: 56     FLGAVMPRLGIGMDTC-----VIPLRHGGLSLVQTTDYIYPIVDDPYMMGRIACANVLSD 110
              FL A++P  G+G+ T      ++  R G L++ +    +  I +    +GR+AC  V   
Sbjct: 899645 FLAAILPGAGMGLMTINNLGQIVSARGGNLAMQEACTTLVSIAN---CLGRLACGRVADS 899815

Query: 111    LYAMG 115
              L A G
Sbjct: 899816 LVARG 899830


>scaffold_12
          Length = 1412611

 Score = 31.2 bits (69), Expect = 2.4
 Identities = 14/34 (41%), Positives = 20/34 (58%)
 Frame = +1

Query: 68      MDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMMGR 101
               +D  V P+RHGG+ +V   D + P+VD     GR
Sbjct: 1329898 LDDEVAPVRHGGVVVVLADDEVAPVVDGRRRRGR 1329999


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.320    0.135    0.405 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 25,202,564
Number of Sequences: 1185
Number of extensions: 387552
Number of successful extensions: 1610
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 48
Number of HSP's successfully gapped in prelim test: 29
Number of HSP's that attempted gapping in prelim test: 1149
Number of HSP's gapped (non-prelim): 851
length of query: 392
length of database: 18,886,859
effective HSP length: 106
effective length of query: 286
effective length of database: 18,761,249
effective search space: 5365717214
effective search space used: 5365717214
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: 64 (29.3 bits)
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 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             60   4e-09
scaffold_56                                                            33   1.00 
scaffold_8                                                             33   1.00 
scaffold_22                                                            32   1.3  
scaffold_4                                                             32   1.7  
scaffold_69                                                            31   3.8  
scaffold_24                                                            31   3.8  
scaffold_10                                                            31   3.8  
scaffold_37                                                            30   5.0  
scaffold_11                                                            30   6.5  
scaffold_73                                                            30   8.5  
scaffold_5                                                             30   8.5  
scaffold_1                                                             30   8.5  

>scaffold_6
          Length = 2755913

 Score = 60.5 bits (145), Expect = 4e-09
 Identities = 83/354 (23%), Positives = 128/354 (36%), Gaps = 33/354 (9%)
 Frame = -1

Query: 108     PSPTFPALGIGM--DSCVI---PLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDL 162
               P P    +G+    D  V+   P +     +V T DFF   V DPY+ G+IA  + LSD 
Sbjct: 1819277 PPPCDVVVGVDAPDDGAVVAYGPEKDQRTQVVHTVDFFRSFVSDPYVFGQIAANHALSDC 1819098

Query: 163     YAMGITECDNXXXXXXXXXXXXXXXREKVTPLMVKGFRDXXXXXXXXXXXXXXXXNPWII 222
                AMG                      E     ++ G R                    + 
Sbjct: 1819097 DAMGALPASALALVVVPYGSNRCV--EDTLFQLMSGARVALDAAGCTLVGGHTCEGAELG 1818924

Query: 223     IGGVATVVCQPNE-FIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDNPERWNKVKMVVSR 281
               +G     V      F+  +  V GD L+LTKPLGT  AV A            ++ +   
Sbjct: 1818923 LGFAVNGVPDGGAAFLSKEGMVAGDALILTKPLGTG-AVFA----------GDMRALSRA 1818777

Query: 282     EEVELAYQEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQRN-EVSFV 340
                EV    +E + +M   NR AA  +      + TD+TGFG+LGH   + +      V   
Sbjct: 1818776 AEV----RECLASMTQSNRAAATELKEHGCTSCTDVTGFGLLGHLAEMCRASGGARVELS 1818609

Query: 341     IHNLPIIAKMAAVSKASGRFGLLQGTS--------------------------AETSGGL 374
               +   P+ A    +++  G F  LQ  +                           +T+GGL
Sbjct: 1818608 LAATPVFAGAKRLAR-DGVFSSLQDENLRAKRAVENHGEAVELDAAAYALLFDPQTAGGL 1818432

Query: 375     LICLPREQAARFCSEIKSSKYGEGHQAWIVGIVEKGNRTARIIDKPRVIEVLPR 428
               L  +P ++A    + I+            V  V+  +  A++    R + VL R
Sbjct: 1818431 LATVPADRADAALAAIRRVGCARAAVIGTVAAVDAAHFPAKVAVV*RELGVLRR 1818270



 Score = 31.2 bits (69), Expect = 2.9
 Identities = 20/53 (37%), Positives = 26/53 (49%)
 Frame = -1

Query: 78      RPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSPTFPALGIGMDSCVIPLRHGG 130
               RP+   P+G    G    A   AGLP+GAGP P   A G+  +  +    HGG
Sbjct: 1795520 RPEASAPVGTCQQGRPPRAP--AGLPSGAGP-PAHTATGLSFEDVI--FTHGG 1795377



 Score = 29.6 bits (65), Expect = 8.5
 Identities = 21/58 (36%), Positives = 27/58 (46%), Gaps = 3/58 (5%)
 Frame = +1

Query: 56      GMKGXGCKVPQEALLKLLAGLTRPDVR--PPLGRGLVGGQEEA-SQEAGLPAGAGPSP 110
               G +G G   P +A    + G  R   R  P  GRG   GQ    ++  G PAG GP+P
Sbjct: 2389063 GRRGAGRARPGDAARARVRGEPRRRARRRPAPGRGRARGQRRRRARRRGHPAGHGPAP 2389236


>scaffold_56
          Length = 168398

 Score = 32.7 bits (73), Expect = 1.00
 Identities = 21/53 (39%), Positives = 23/53 (43%)
 Frame = -3

Query: 58    KGXGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSP 110
             KG G   P+       AGL R D R P GRGL   +         P GAGP P
Sbjct: 24621 KGRGLPAPEREPAAEPAGLRRGDAREP-GRGLPRRRGPQEGRGPRPGGAGPLP 24466


>scaffold_8
          Length = 1961717

 Score = 32.7 bits (73), Expect = 1.00
 Identities = 24/60 (40%), Positives = 27/60 (45%)
 Frame = +2

Query: 51      LTGFSGMKGXGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSP 110
               L G  G +  G +V   A      G  RPD  PP  R    GQEE     GLPAG  P+P
Sbjct: 1561985 LPGEDGPRRQGRRVGPAAAAP---GDLRPDALPPQVRRQADGQEEVQ---GLPAGRRPAP 1562146


>scaffold_22
          Length = 725242

 Score = 32.3 bits (72), Expect = 1.3
 Identities = 22/45 (48%), Positives = 24/45 (53%), Gaps = 2/45 (4%)
 Frame = +3

Query: 72    LLAGLTRPDVRPPLGR--GLVGGQEEASQEAGLPAGAGPSPTFPA 114
             LLA   RP  RPP  R  GLV G    ++  GL  GAGP P  PA
Sbjct: 19860 LLAAAARPAPRPPRARRHGLVRG----ARGDGLRRGAGPRPGAPA 19982


>scaffold_4
          Length = 3024044

 Score = 32.0 bits (71), Expect = 1.7
 Identities = 26/82 (31%), Positives = 35/82 (42%), Gaps = 7/82 (8%)
 Frame = -2

Query: 239     PDSAVVGDVLVLTKPLGTQVAVNAHQWLDN-----PERWNKVKMV--VSREEVELAYQEA 291
               P  AVVGD  V T+ LGT+    +  W         +RW + K V    R +VEL     
Sbjct: 2645338 PPPAVVGDARVYTRVLGTKHVPQSPDWNCGLACWICDRWRETKFVWAPGRSDVELDGGAK 2645159

Query: 292     MFNMATLNRTAAGLMHTFNAHA 313
                 + A  NR A  +   F+  A
Sbjct: 2645158 ARDAAASNRLAVQVRCAFDGWA 2645093



 Score = 32.0 bits (71), Expect = 1.7
 Identities = 31/104 (29%), Positives = 44/104 (42%), Gaps = 9/104 (8%)
 Frame = +3

Query: 35      YRPFEPQ-ALGLSPSWRLTGFSGMKGXGCKVPQEA---LLKLLAGLTR-----PDVRPPL 85
               +RP  P  +  +    + T F G+KG G ++P  A   L +    + R     P V P +
Sbjct: 2239392 HRPLRPGGSTRVQRKGKRT*FVGLKGEGPELPGVAGVALARRPGDVGRHVEDVPAVAPRV 2239571

Query: 86      GRGLVGGQEEASQEAGLPAGAGPSPTFPALGIGMDSCVIPLRHG 129
               G  L+G  E   +   L  GAG     PA   G+D    P  HG
Sbjct: 2239572 GGLLLGDDEAQLRRVDLVGGAGRGRLGPAEDGGVDGVEGPRGHG 2239703



 Score = 31.6 bits (70), Expect = 2.2
 Identities = 14/27 (51%), Positives = 16/27 (59%)
 Frame = +2

Query: 103    PAGAGPSPTFPALGIGMDSCVIPLRHG 129
              P+G GP+P  PA G G  SC  P R G
Sbjct: 114905 PSGPGPAPRGPAGGPGTPSCTPPRRRG 114985


>scaffold_69
          Length = 78151

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 16/31 (51%), Positives = 17/31 (54%)
 Frame = -3

Query: 75    GLTRPDVRPPLGRGLVGGQEEASQEAGLPAG 105
             G  RP  R   GRG VG +EEA QE G   G
Sbjct: 60266 GARRPGGRESRGRGAVGVEEEAQQEEGRRQG 60174


>scaffold_24
          Length = 564861

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 24/67 (35%), Positives = 31/67 (46%), Gaps = 8/67 (11%)
 Frame = -1

Query: 51     LTGFSGMKGXGCKVPQEALLKL-LAGLTRPDVRPPLGRGLVGGQEEASQEAG-------L 102
              ++G SG +G G ++ QEA     + G T       LGRG VGG    S   G       L
Sbjct: 490980 MSGSSGARGQGGELEQEAGRGAQVGGHTEISAGESLGRGAVGGMTSLSAGVGRDGGGLIL 490801

Query: 103    PAGAGPS 109
               AG+G S
Sbjct: 490800 SAGSGKS 490780



 Score = 30.8 bits (68), Expect = 3.8
 Identities = 24/67 (35%), Positives = 31/67 (46%), Gaps = 8/67 (11%)
 Frame = -1

Query: 51     LTGFSGMKGXGCKVPQEALLKL-LAGLTRPDVRPPLGRGLVGGQEEASQEAG-------L 102
              ++G SG +G G ++ QEA     + G T       LGRG VGG    S   G       L
Sbjct: 483825 MSGSSGARGQGGELEQEAGRGAQVGGHTEISAGESLGRGAVGGMTSLSAGVGRDGGGLIL 483646

Query: 103    PAGAGPS 109
               AG+G S
Sbjct: 483645 SAGSGKS 483625


>scaffold_10
          Length = 1654512

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 27/79 (34%), Positives = 35/79 (44%), Gaps = 8/79 (10%)
 Frame = -1

Query: 61      GCKVPQ-EALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSPTFPALG--I 117
               G + PQ  A L+  A    P   PP  RGL   +EEA++E G  AG         LG  +
Sbjct: 1418976 GLQGPQAHAPLRARAEPEGPRPGPPQARGLRAAEEEAAKEEGEEAGVSQEEAGAPLGSRL 1418797

Query: 118     GMDSCVIP-----LRHGGL 131
                M + V+        HGGL
Sbjct: 1418796 PMGAAVVEQLVDVADHGGL 1418740


>scaffold_37
          Length = 376641

 Score = 30.4 bits (67), Expect = 5.0
 Identities = 21/68 (30%), Positives = 29/68 (42%), Gaps = 3/68 (4%)
 Frame = -1

Query: 46   SPSWRLTGFSGMKGXGCKVPQEALLKLLA---GLTRPDVRPPLGRGLVGGQEEASQEAGL 102
            SP WR +          ++P+ +  +  A   G TR   R P   G   G   A  EAG+
Sbjct: 7356 SPQWRQSSTR-------RLPRGSPTRAAAAAHGPTRAAQRAPHAAGATPGPRRARPEAGM 7198

Query: 103  PAGAGPSP 110
            P+GA   P
Sbjct: 7197 PSGAASGP 7174


>scaffold_11
          Length = 1405779

 Score = 30.0 bits (66), Expect = 6.5
 Identities = 22/53 (41%), Positives = 24/53 (45%)
 Frame = +3

Query: 78     RPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSPTFPALGIGMDSCVIPLRHGG 130
              RP    P GRG  GG+  A   AG  A AGP P   A   G +    P  HGG
Sbjct: 785214 RPARDAPAGRGAPGGRGAARGRAG-RAAAGPGPRGRARAAGANR---PPFHGG 785360


>scaffold_73
          Length = 105214

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 18/34 (52%), Positives = 19/34 (55%), Gaps = 3/34 (8%)
 Frame = -3

Query: 78    RPDVRPPLGRG---LVGGQEEASQEAGLPAGAGP 108
             RP   PP GRG     GGQ    + AG PAGAGP
Sbjct: 91226 RPGASPPGGRGSPDAAGGQP--GRGAGNPAGAGP 91131


>scaffold_5
          Length = 3000003

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 23/65 (35%), Positives = 27/65 (41%), Gaps = 3/65 (4%)
 Frame = +2

Query: 49      WRLTGFSGMKGXGCKVPQEALLKLLAGLTRPDVRPPLGRG---LVGGQEEASQEAGLPAG 105
               W LTG +            A  K   G  RP  R   GRG    +G +  AS+  G PAG
Sbjct: 1324232 WELTGLAAAGSGRG*GAGRAHKKGSRGFIRPG*R---GRGPCRRMGSRRRASRPRGAPAG 1324402

Query: 106     AGPSP 110
                GP P
Sbjct: 1324403 PGPGP 1324417


>scaffold_1
          Length = 3520325

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 16/44 (36%), Positives = 21/44 (47%)
 Frame = -3

Query: 75      GLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGPSPTFPALGIG 118
               GL R ++ P LGR L+   +E   E  L  G G       LG+G
Sbjct: 2765031 GLQRRELGPDLGRDLLVEDDERGVERALAPGPGRERVLLVLGLG 2764900


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
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: 27,654,991
Number of Sequences: 1185
Number of extensions: 479513
Number of successful extensions: 4839
Number of sequences better than 10.0: 14
Number of HSP's better than 10.0 without gapping: 154
Number of HSP's successfully gapped in prelim test: 69
Number of HSP's that attempted gapping in prelim test: 1223
Number of HSP's gapped (non-prelim): 4173
length of query: 448
length of database: 18,886,859
effective HSP length: 107
effective length of query: 341
effective length of database: 18,760,064
effective search space: 6397181824
effective search space used: 6397181824
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: 64 (29.3 bits)
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= SPP00000072_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Pan troglodytes # Complete
         (391 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             82   1e-15
scaffold_15                                                            32   1.1  
scaffold_12                                                            31   2.4  

>scaffold_6
          Length = 2755913

 Score = 82.0 bits (201), Expect = 1e-15
 Identities = 76/288 (26%), Positives = 121/288 (42%), Gaps = 28/288 (9%)
 Frame = -1

Query: 81      LVQTTDYIYPIVDDPYMMGRIACANVLSDLYAMGVTECDNMLMLLGVSNKMTDRERDKVM 140
               +V T D+    V DPY+ G+IA  + LSD  AMG      + ++  V    ++R  +  +
Sbjct: 1819187 VVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALV--VVPYGSNRCVEDTL 1819014

Query: 141     PLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNE-FIMPDNAVPGDVLVLT 199
                 ++ G + A + AG ++ GG T     + LG     V      F+  +  V GD L+LT
Sbjct: 1819013 FQLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILT 1818834

Query: 200     KPLGTQVAVAVHQWLDIPEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTAAGLMHTFNA 259
               KPLGT    A             ++ +    +V    +E + +M + NR AA  +     
Sbjct: 1818833 KPLGTGAVFA-----------GDMRALSRAAEV----RECLASMTQSNRAAATELKEHGC 1818699

Query: 260     HAATDITGFGILGHAQNLAKQQRN-EVSFVIHNLPVLAKMAAVSKACGNMFGLM------ 312
                + TD+TGFG+LGH   + +      V   +   PV A   A   A   +F  +      
Sbjct: 1818698 TSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFA--GAKRLARDGVFSSLQDENLR 1818525

Query: 313     -------HGTC-------------PETSGGLLICLPREQAARFCAEIK 340
                      HG               P+T+GGLL  +P ++A    A I+
Sbjct: 1818524 AKRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIR 1818381



 Score = 29.3 bits (64), Expect = 9.2
 Identities = 17/58 (29%), Positives = 26/58 (44%)
 Frame = +1

Query: 33      PQDVLQKLLESLQENHFQEDEQFLGAVMPRLGIGMDTCVIPLRHGGLSLVQTTDYIYP 90
               PQ+VL++  +   E   +++E  LG   P L +     V+P R G   L    D   P
Sbjct: 2275972 PQEVLERQTQDEDERPDEQEEVALGDRAPGLALAPGVAVLPHRFGEGRLAHEGDDFLP 2276145


>scaffold_15
          Length = 987846

 Score = 32.3 bits (72), Expect = 1.1
 Identities = 19/65 (29%), Positives = 32/65 (49%), Gaps = 5/65 (7%)
 Frame = +2

Query: 55     FLGAVMPRLGIGMDTC-----VIPLRHGGLSLVQTTDYIYPIVDDPYMMGRIACANVLSD 109
              FL A++P  G+G+ T      ++  R G L++ +    +  I +    +GR+AC  V   
Sbjct: 899645 FLAAILPGAGMGLMTINNLGQIVSARGGNLAMQEACTTLVSIAN---CLGRLACGRVADS 899815

Query: 110    LYAMG 114
              L A G
Sbjct: 899816 LVARG 899830


>scaffold_12
          Length = 1412611

 Score = 31.2 bits (69), Expect = 2.4
 Identities = 14/34 (41%), Positives = 20/34 (58%)
 Frame = +1

Query: 67      MDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMMGR 100
               +D  V P+RHGG+ +V   D + P+VD     GR
Sbjct: 1329898 LDDEVAPVRHGGVVVVLADDEVAPVVDGRRRRGR 1329999


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.320    0.135    0.404 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 25,155,390
Number of Sequences: 1185
Number of extensions: 387096
Number of successful extensions: 1603
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 48
Number of HSP's successfully gapped in prelim test: 27
Number of HSP's that attempted gapping in prelim test: 1150
Number of HSP's gapped (non-prelim): 841
length of query: 391
length of database: 18,886,859
effective HSP length: 106
effective length of query: 285
effective length of database: 18,761,249
effective search space: 5346955965
effective search space used: 5346955965
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: 64 (29.3 bits)
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= SPP00000073_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Pan troglodytes # Complete
         (448 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             62   2e-09
scaffold_8                                                             33   1.00 
scaffold_4                                                             32   1.7  
scaffold_69                                                            31   3.8  
scaffold_10                                                            31   3.8  
scaffold_24                                                            30   6.5  
scaffold_37                                                            30   8.5  
scaffold_22                                                            30   8.5  
scaffold_13                                                            30   8.5  
scaffold_9                                                             30   8.5  

>scaffold_6
          Length = 2755913

 Score = 62.0 bits (149), Expect = 2e-09
 Identities = 92/398 (23%), Positives = 138/398 (34%), Gaps = 28/398 (7%)
 Frame = -1

Query: 59      GXGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGADPSPTFPALGIG 118
               G G KV   AL + LA L +                  + + G P   D      A   G
Sbjct: 1819361 GCGAKVGASALTRALANLPK----------------SLTFDNGTPPPCDVVVGVDAPDDG 1819230

Query: 119     MDSCVIPLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGITECDNXXXXXX 178
                     P +     +V T DFF   V DPY+ G+IA  + LSD  AMG            
Sbjct: 1819229 AVVAYGPEKDQRTQVVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALVVV 1819050

Query: 179     XXXXXXXXXREKVTPLMVKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVATVVCQPNE-FI 237
                         E     ++ G R                    + +G     V      F+
Sbjct: 1819049 PYGSNRCV--EDTLFQLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFL 1818876

Query: 238     MPDSAVVGDVLVLTKPLGTQVAVNAHQWLDNPERWNKVKMVVSREEVELAYQEAMFNMAT 297
                 +  V GD L+LTKPLGT  AV A            ++ +    EV    +E + +M  
Sbjct: 1818875 SKEGMVAGDALILTKPLGTG-AVFA----------GDMRALSRAAEV----RECLASMTQ 1818741

Query: 298     LNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQRN-EVSFVIHNLPIIAKMAAVSKA 356
                NR AA  +      + TD+TGFG+LGH   + +      V   +   P+ A    +++ 
Sbjct: 1818740 SNRAAATELKEHGCTSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFAGAKRLAR- 1818564

Query: 357     SGRFGLLQGTS--------------------------AETSGGLLICLPREQAARFCSEI 390
                G F  LQ  +                           +T+GGLL  +P ++A    + I
Sbjct: 1818563 DGVFSSLQDENLRAKRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAI 1818384

Query: 391     KSSKYGEGHQAWIVGIVEKGNRTARIIDKPRVIEVLPR 428
               +            V  V+  +  A++    R + VL R
Sbjct: 1818383 RRVGCARAAVIGTVAAVDAAHFPAKVAVV*RELGVLRR 1818270


>scaffold_8
          Length = 1961717

 Score = 32.7 bits (73), Expect = 1.00
 Identities = 24/60 (40%), Positives = 27/60 (45%)
 Frame = +2

Query: 51      LTGFSGMKGXGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGADPSP 110
               L G  G +  G +V   A      G  RPD  PP  R    GQEE     GLPAG  P+P
Sbjct: 1561985 LPGEDGPRRQGRRVGPAAAAP---GDLRPDALPPQVRRQADGQEEVQ---GLPAGRRPAP 1562146


>scaffold_4
          Length = 3024044

 Score = 32.0 bits (71), Expect = 1.7
 Identities = 26/82 (31%), Positives = 35/82 (42%), Gaps = 7/82 (8%)
 Frame = -2

Query: 239     PDSAVVGDVLVLTKPLGTQVAVNAHQWLDN-----PERWNKVKMV--VSREEVELAYQEA 291
               P  AVVGD  V T+ LGT+    +  W         +RW + K V    R +VEL     
Sbjct: 2645338 PPPAVVGDARVYTRVLGTKHVPQSPDWNCGLACWICDRWRETKFVWAPGRSDVELDGGAK 2645159

Query: 292     MFNMATLNRTAAGLMHTFNAHA 313
                 + A  NR A  +   F+  A
Sbjct: 2645158 ARDAAASNRLAVQVRCAFDGWA 2645093


>scaffold_69
          Length = 78151

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 16/31 (51%), Positives = 17/31 (54%)
 Frame = -3

Query: 75    GLTRPDVRPPLGRGLVGGQEEASQEAGLPAG 105
             G  RP  R   GRG VG +EEA QE G   G
Sbjct: 60266 GARRPGGRESRGRGAVGVEEEAQQEEGRRQG 60174


>scaffold_10
          Length = 1654512

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 27/79 (34%), Positives = 35/79 (44%), Gaps = 8/79 (10%)
 Frame = -1

Query: 61      GCKVPQ-EALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGADPSPTFPALG--I 117
               G + PQ  A L+  A    P   PP  RGL   +EEA++E G  AG         LG  +
Sbjct: 1418976 GLQGPQAHAPLRARAEPEGPRPGPPQARGLRAAEEEAAKEEGEEAGVSQEEAGAPLGSRL 1418797

Query: 118     GMDSCVIP-----LRHGGL 131
                M + V+        HGGL
Sbjct: 1418796 PMGAAVVEQLVDVADHGGL 1418740


>scaffold_24
          Length = 564861

 Score = 30.0 bits (66), Expect = 6.5
 Identities = 20/56 (35%), Positives = 26/56 (46%), Gaps = 1/56 (1%)
 Frame = -1

Query: 51     LTGFSGMKGXGCKVPQEALLKL-LAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAG 105
              ++G SG +G G ++ QEA     + G T       LGRG VGG    S   G   G
Sbjct: 490980 MSGSSGARGQGGELEQEAGRGAQVGGHTEISAGESLGRGAVGGMTSLSAGVGRDGG 490813



 Score = 30.0 bits (66), Expect = 6.5
 Identities = 20/56 (35%), Positives = 26/56 (46%), Gaps = 1/56 (1%)
 Frame = -1

Query: 51     LTGFSGMKGXGCKVPQEALLKL-LAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAG 105
              ++G SG +G G ++ QEA     + G T       LGRG VGG    S   G   G
Sbjct: 483825 MSGSSGARGQGGELEQEAGRGAQVGGHTEISAGESLGRGAVGGMTSLSAGVGRDGG 483658


>scaffold_37
          Length = 376641

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 21/68 (30%), Positives = 29/68 (42%), Gaps = 3/68 (4%)
 Frame = -1

Query: 46   SPSWRLTGFSGMKGXGCKVPQEALLKLLA---GLTRPDVRPPLGRGLVGGQEEASQEAGL 102
            SP WR +          ++P+ +  +  A   G TR   R P   G   G   A  EAG+
Sbjct: 7356 SPQWRQSSTR-------RLPRGSPTRAAAAAHGPTRAAQRAPHAAGATPGPRRARPEAGM 7198

Query: 103  PAGADPSP 110
            P+GA   P
Sbjct: 7197 PSGAASGP 7174


>scaffold_22
          Length = 725242

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 21/45 (46%), Positives = 23/45 (51%), Gaps = 2/45 (4%)
 Frame = +3

Query: 72    LLAGLTRPDVRPPLGR--GLVGGQEEASQEAGLPAGADPSPTFPA 114
             LLA   RP  RPP  R  GLV G    ++  GL  GA P P  PA
Sbjct: 19860 LLAAAARPAPRPPRARRHGLVRG----ARGDGLRRGAGPRPGAPA 19982


>scaffold_13
          Length = 1338096

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 17/37 (45%), Positives = 18/37 (48%)
 Frame = -3

Query: 78     RPDVRPPLGRGLVGGQEEASQEAGLPAGADPSPTFPA 114
              RP VRPP           A +  GLPA ADP P  PA
Sbjct: 530476 RPRVRPP-----------ADEPPGLPADADPRPRRPA 530399


>scaffold_9
          Length = 1602328

 Score = 29.6 bits (65), Expect = 8.5
 Identities = 15/44 (34%), Positives = 22/44 (50%)
 Frame = +2

Query: 72      LLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGADPSPTFPAL 115
               L +G +RP   PP  R L   +   S+ A   +GA   PT P++
Sbjct: 1224089 LKSGTSRPFPAPPGARNLRASESAPSKAASSSSGAQARPTAPSM 1224220


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.320    0.136    0.411 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 26,946,217
Number of Sequences: 1185
Number of extensions: 453703
Number of successful extensions: 3936
Number of sequences better than 10.0: 11
Number of HSP's better than 10.0 without gapping: 129
Number of HSP's successfully gapped in prelim test: 60
Number of HSP's that attempted gapping in prelim test: 1283
Number of HSP's gapped (non-prelim): 3227
length of query: 448
length of database: 18,886,859
effective HSP length: 107
effective length of query: 341
effective length of database: 18,760,064
effective search space: 6397181824
effective search space used: 6397181824
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: 64 (29.3 bits)
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= SPP00000107_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Mus musculus # Complete
         (392 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             82   1e-15
scaffold_15                                                            32   1.1  
scaffold_12                                                            31   2.4  

>scaffold_6
          Length = 2755913

 Score = 82.0 bits (201), Expect = 1e-15
 Identities = 76/288 (26%), Positives = 121/288 (42%), Gaps = 28/288 (9%)
 Frame = -1

Query: 82      LVQTTDYIYPIVDDPYMMGRIACANVLSDLYAMGVTECDNMLMLLGVSNKMTDRERDKVI 141
               +V T D+    V DPY+ G+IA  + LSD  AMG      + ++  V    ++R  +  +
Sbjct: 1819187 VVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALV--VVPYGSNRCVEDTL 1819014

Query: 142     PLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNE-FIMPDNAVPGDVLVLT 200
                 ++ G + A + AG ++ GG T     + LG     V      F+  +  V GD L+LT
Sbjct: 1819013 FQLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILT 1818834

Query: 201     KPLGTQVAVAVHQWLDIPEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTAAGLMHTFNA 260
               KPLGT    A             ++ +    +V    +E + +M + NR AA  +     
Sbjct: 1818833 KPLGTGAVFA-----------GDMRALSRAAEV----RECLASMTQSNRAAATELKEHGC 1818699

Query: 261     HAATDITGFGILGHAQNLAKQQRN-EVSFVIHNLPVLAKMAAVSKACGNMFGLM------ 313
                + TD+TGFG+LGH   + +      V   +   PV A   A   A   +F  +      
Sbjct: 1818698 TSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFA--GAKRLARDGVFSSLQDENLR 1818525

Query: 314     -------HGTC-------------PETSGGLLICLPREQAARFCAEIK 341
                      HG               P+T+GGLL  +P ++A    A I+
Sbjct: 1818524 AKRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIR 1818381



 Score = 29.3 bits (64), Expect = 9.3
 Identities = 17/58 (29%), Positives = 26/58 (44%)
 Frame = +1

Query: 34      PQDVLQKLLESLQENHFQEDEQFLGAVMPRLGIGMDTCVIPLRHGGLSLVQTTDYIYP 91
               PQ+VL++  +   E   +++E  LG   P L +     V+P R G   L    D   P
Sbjct: 2275972 PQEVLERQTQDEDERPDEQEEVALGDRAPGLALAPGVAVLPHRFGEGRLAHEGDDFLP 2276145


>scaffold_15
          Length = 987846

 Score = 32.3 bits (72), Expect = 1.1
 Identities = 19/65 (29%), Positives = 32/65 (49%), Gaps = 5/65 (7%)
 Frame = +2

Query: 56     FLGAVMPRLGIGMDTC-----VIPLRHGGLSLVQTTDYIYPIVDDPYMMGRIACANVLSD 110
              FL A++P  G+G+ T      ++  R G L++ +    +  I +    +GR+AC  V   
Sbjct: 899645 FLAAILPGAGMGLMTINNLGQIVSARGGNLAMQEACTTLVSIAN---CLGRLACGRVADS 899815

Query: 111    LYAMG 115
              L A G
Sbjct: 899816 LVARG 899830


>scaffold_12
          Length = 1412611

 Score = 31.2 bits (69), Expect = 2.4
 Identities = 14/34 (41%), Positives = 20/34 (58%)
 Frame = +1

Query: 68      MDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMMGR 101
               +D  V P+RHGG+ +V   D + P+VD     GR
Sbjct: 1329898 LDDEVAPVRHGGVVVVLADDEVAPVVDGRRRRGR 1329999


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.320    0.135    0.405 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 25,231,875
Number of Sequences: 1185
Number of extensions: 388946
Number of successful extensions: 1615
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 48
Number of HSP's successfully gapped in prelim test: 30
Number of HSP's that attempted gapping in prelim test: 1155
Number of HSP's gapped (non-prelim): 849
length of query: 392
length of database: 18,886,859
effective HSP length: 106
effective length of query: 286
effective length of database: 18,761,249
effective search space: 5365717214
effective search space used: 5365717214
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: 64 (29.3 bits)
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 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             57   4e-08
scaffold_4                                                             31   3.8  

>scaffold_6
          Length = 2755913

 Score = 57.4 bits (137), Expect = 4e-08
 Identities = 52/200 (26%), Positives = 76/200 (38%), Gaps = 1/200 (0%)
 Frame = -1

Query: 136     LVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGITECDNXXXXXXXXXXXXEKEREKVT 195
               +V T DFF   V DPY+ G+IA  + LSD  AMG                   +  E   
Sbjct: 1819187 VVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPAS--ALALVVVPYGSNRCVEDTL 1819014

Query: 196     PLMIKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVATVVCQQNE-FIMPDSAVVGDVLVLT 254
                 ++ G R                    + +G     V      F+  +  V GD L+LT
Sbjct: 1819013 FQLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILT 1818834

Query: 255     KPLGTQVAANAHQWLDNPEKWNKIKMVVSREEVELAYQEAMFNMATLNRTAAGLMHTFNA 314
               KPLGT                  ++ +    EV    +E + +M   NR AA  +     
Sbjct: 1818833 KPLGTGAVFAGD-----------MRALSRAAEV----RECLASMTQSNRAAATELKEHGC 1818699

Query: 315     HAATDITGFGILGHSQNLAK 334
                + TD+TGFG+LGH   + +
Sbjct: 1818698 TSCTDVTGFGLLGHLAEMCR 1818639


>scaffold_4
          Length = 3024044

 Score = 30.8 bits (68), Expect = 3.8
 Identities = 25/82 (30%), Positives = 35/82 (42%), Gaps = 7/82 (8%)
 Frame = -2

Query: 242     PDSAVVGDVLVLTKPLGTQVAANAHQWLDN-----PEKWNKIKMV--VSREEVELAYQEA 294
               P  AVVGD  V T+ LGT+    +  W         ++W + K V    R +VEL     
Sbjct: 2645338 PPPAVVGDARVYTRVLGTKHVPQSPDWNCGLACWICDRWRETKFVWAPGRSDVELDGGAK 2645159

Query: 295     MFNMATLNRTAAGLMHTFNAHA 316
                 + A  NR A  +   F+  A
Sbjct: 2645158 ARDAAASNRLAVQVRCAFDGWA 2645093


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
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: 25,556,942
Number of Sequences: 1185
Number of extensions: 363395
Number of successful extensions: 1717
Number of sequences better than 10.0: 2
Number of HSP's better than 10.0 without gapping: 46
Number of HSP's successfully gapped in prelim test: 25
Number of HSP's that attempted gapping in prelim test: 1092
Number of HSP's gapped (non-prelim): 928
length of query: 452
length of database: 18,886,859
effective HSP length: 108
effective length of query: 344
effective length of database: 18,758,879
effective search space: 6453054376
effective search space used: 6453054376
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)
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= SPP00000047_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Drosophila melanogaster # Complete
         (398 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             88   2e-17
scaffold_9                                                             33   0.50 
scaffold_4                                                             33   0.50 
scaffold_25                                                            32   1.5  
scaffold_20                                                            32   1.9  
scaffold_12                                                            31   3.3  
scaffold_1                                                             31   3.3  
scaffold_15                                                            30   7.2  
scaffold_70                                                            29   9.5  
scaffold_56                                                            29   9.5  
scaffold_5                                                             29   9.5  

>scaffold_6
          Length = 2755913

 Score = 88.2 bits (217), Expect = 2e-17
 Identities = 97/360 (26%), Positives = 153/360 (42%), Gaps = 31/360 (8%)
 Frame = -1

Query: 48      GRGCKVPQDVLSKLVSALQQDYSAQDQEPQFLNVAIPRIGIGLDCSVI---PLRHGGLCL 104
               G G KV    L++ ++ L +  +  +  P   +V +  +    D +V+   P +     +
Sbjct: 1819361 GCGAKVGASALTRALANLPKSLTFDNGTPPPCDVVVG-VDAPDDGAVVAYGPEKDQRTQV 1819185

Query: 105     VQTTDFFYPIVDDPYMMGKIACANVLSDLYAMGVTDCDNMLMLLAVSTKMTEKERDVVIP 164
               V T DFF   V DPY+ G+IA  + LSD  AMG     + L L+ V         D +  
Sbjct: 1819184 VHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPA-SALALVVVPYGSNRCVEDTLFQ 1819008

Query: 165     LIMRGFKDSALEAGTTVTGGQSVVNPWCTIGGVASTICQPN-EYIVPDNAVVGDVLVLTK 223
               L M G + +   AG T+ GG +       +G   + +      ++  +  V GD L+LTK
Sbjct: 1819007 L-MSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILTK 1818831

Query: 224     PLGTQVAVNAHQWIDQPERWNRIKLVVSEEDVRKAYHRAMNSMARLNRVAARLMHKYNAH 283
               PLGT  AV A            ++ +    +VR+     + SM + NR AA  + ++   
Sbjct: 1818830 PLGTG-AVFA----------GDMRALSRAAEVRE----CLASMTQSNRAAATELKEHGCT 1818696

Query: 284     GATDITGFGLLGH-AQTLAAHQKKDVSFVIHNLPVIAKMAAVAKACGNMFQLLQGHS--- 339
                 TD+TGFGLLGH A+   A     V   +   PV A    +A+    +F  LQ  +   
Sbjct: 1818695 SCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFAGAKRLAR--DGVFSSLQDENLRA 1818522

Query: 340     -----------------------AETSGGLLICLPREQAAAYCKDIEKQEGYQAWIIGIV 376
                                       +T+GGLL  +P ++A A    I +    +A +IG V
Sbjct: 1818521 KRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIRRVGCARAAVIGTV 1818342


>scaffold_9
          Length = 1602328

 Score = 33.5 bits (75), Expect = 0.50
 Identities = 24/64 (37%), Positives = 32/64 (50%), Gaps = 1/64 (1%)
 Frame = +1

Query: 260    HRAMNSMARLNRVAARLMHKYNAHGATDITGFGLLGHAQTLAAHQKKDVSFVIHN-LPVI 318
              HR    + RL R  AR      A GA  +TG G LGHA+    H +  V  VI + +P  
Sbjct: 652615 HRRQQLLPRLGRQGARDFCG-EAQGAAGVTGHGALGHARDGRGHGRGRVGHVIFSRIPRG 652791

Query: 319    AKMA 322
              A++A
Sbjct: 652792 AELA 652803


>scaffold_4
          Length = 3024044

 Score = 33.5 bits (75), Expect = 0.50
 Identities = 23/64 (35%), Positives = 32/64 (50%)
 Frame = -3

Query: 1     MSYAADVLNSAHLELHGGGDAELRRPFDPTAHDLDASFRLTRFADLKGRGCKVPQDVLSK 60
             +S   D  +SA L+L G G+ + ++  D   H  DAS RL R     G  C+VP D    
Sbjct: 93684 VSRPRDQNDSAALDLEGVGEMDWKQVKD-LCHGHDASVRLGRVPQPVGSRCRVPFDAREH 93508

Query: 61    LVSA 64
              V+A
Sbjct: 93507 AVAA 93496


>scaffold_25
          Length = 552203

 Score = 32.0 bits (71), Expect = 1.5
 Identities = 22/77 (28%), Positives = 35/77 (45%), Gaps = 3/77 (3%)
 Frame = -2

Query: 223    KPLGTQVAVNAHQWIDQPE-RWNRIKLVVSEEDVRKAYHRAMNSMARLNRVAARLMHKYN 281
              +P G + A   H + D+   R  R+      E  R    R ++++ R  R   RL  + +
Sbjct: 309103 RPPGRRRAAQGHDFPDEAAPRGLRLLGAPQREQRRNVLRRPVDALPRFGRRPLRLRPRAS 308924

Query: 282    AH--GATDITGFGLLGH 296
              A   G  D+ GFG +GH
Sbjct: 308923 APRAGQDDVRGFGAVGH 308873


>scaffold_20
          Length = 721792

 Score = 31.6 bits (70), Expect = 1.9
 Identities = 19/55 (34%), Positives = 26/55 (47%)
 Frame = -1

Query: 174    ALEAGTTVTGGQSVVNPWCTIGGVASTICQPNEYIVPDNAVVGDVLVLTKPLGTQ 228
              A+ A   V  G  V       G V +T+ QP   ++P  AVVG V+ +  P G Q
Sbjct: 267961 AMSAEMPVVQGAVVAPGQVQPGIVQATVLQPVAVVIPQGAVVGTVVAVQGPGGPQ 267797


>scaffold_12
          Length = 1412611

 Score = 30.8 bits (68), Expect = 3.3
 Identities = 20/58 (34%), Positives = 28/58 (48%)
 Frame = +1

Query: 66      QQDYSAQDQEPQFLNVAIPRIGIGLDCSVIPLRHGGLCLVQTTDFFYPIVDDPYMMGK 123
               Q++  A+D  P  L       G  LD  V P+RHGG+ +V   D   P+VD     G+
Sbjct: 1329844 QEEVLAEDFGPGGLR------GPVLDDEVAPVRHGGVVVVLADDEVAPVVDGRRRRGR 1329999


>scaffold_1
          Length = 3520325

 Score = 30.8 bits (68), Expect = 3.3
 Identities = 20/68 (29%), Positives = 32/68 (47%)
 Frame = +1

Query: 300    LAAHQKKDVSFVIHNLPVIAKMAAVAKACGNMFQLLQGHSAETSGGLLICLPREQAAAYC 359
              L AH K+     +  + ++ + A    ACG+  +LLQ  +AE   GL +      +AA  
Sbjct: 711415 LVAHSKRGSGVTV--MALVDQAACNKSACGDYVKLLQDLTAENLKGLNVQPTTNPSAAVG 711588

Query: 360    KDIEKQEG 367
                 E+Q G
Sbjct: 711589 AHFERQRG 711612


>scaffold_15
          Length = 987846

 Score = 29.6 bits (65), Expect = 7.2
 Identities = 18/55 (32%), Positives = 28/55 (50%)
 Frame = +1

Query: 344    GGLLICLPREQAAAYCKDIEKQEGYQAWIIGIVEKGNKTARIIDKPRVIEVPAKD 398
              G  ++  P    A   +D+E   G    ++ +V+ G + ARI   PRV+ VPA D
Sbjct: 479773 GAGVVVDPGHVEARVPRDVEVGAGMMVVMMIVVDPGRREARI---PRVLRVPAAD 479928


>scaffold_70
          Length = 76094

 Score = 29.3 bits (64), Expect = 9.5
 Identities = 14/41 (34%), Positives = 23/41 (56%)
 Frame = +3

Query: 213   AVVGDVLVLTKPLGTQVAVNAHQWIDQPERWNRIKLVVSEE 253
             AV+ D    + PL    A+ + QW+  P  +NR++L +S E
Sbjct: 27978 AVL*DAFFSSTPLRFSTALRSCQWLPYPIIYNRLELDISSE 28100


>scaffold_56
          Length = 168398

 Score = 29.3 bits (64), Expect = 9.5
 Identities = 14/41 (34%), Positives = 23/41 (56%)
 Frame = -1

Query: 213   AVVGDVLVLTKPLGTQVAVNAHQWIDQPERWNRIKLVVSEE 253
             AV+ D    + PL    A+ + QW+  P  +NR++L +S E
Sbjct: 82745 AVL*DAFFSSTPLRFSTALRSCQWLPYPIIYNRLELDISSE 82623


>scaffold_5
          Length = 3000003

 Score = 29.3 bits (64), Expect = 9.5
 Identities = 19/77 (24%), Positives = 35/77 (45%)
 Frame = -2

Query: 230    AVNAHQWIDQPERWNRIKLVVSEEDVRKAYHRAMNSMARLNRVAARLMHKYNAHGATDIT 289
              AV   + +D+ E   R    ++  D+ +   R + + A    V  RL H  +AHG   + 
Sbjct: 235955 AVEVQRAVDRREL--RAVAPLARRDLERVALRVLLAHAVGRAVRGRLEHLVHAHGHDQVR 235782

Query: 290    GFGLLGHAQTLAAHQKK 306
                G LGH + +   +++
Sbjct: 235781 VRGALGHGRAVGRRRRR 235731


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.321    0.136    0.409 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 26,085,409
Number of Sequences: 1185
Number of extensions: 395271
Number of successful extensions: 1987
Number of sequences better than 10.0: 11
Number of HSP's better than 10.0 without gapping: 76
Number of HSP's successfully gapped in prelim test: 32
Number of HSP's that attempted gapping in prelim test: 988
Number of HSP's gapped (non-prelim): 1349
length of query: 398
length of database: 18,886,859
effective HSP length: 106
effective length of query: 292
effective length of database: 18,761,249
effective search space: 5478284708
effective search space used: 5478284708
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.9 bits)
S2: 64 (29.3 bits)
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: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             69   1e-11
scaffold_1000                                                          32   1.0  
scaffold_2                                                             29   8.6  

>scaffold_6
          Length = 2755913

 Score = 68.9 bits (167), Expect = 1e-11
 Identities = 84/316 (26%), Positives = 126/316 (39%), Gaps = 45/316 (14%)
 Frame = -1

Query: 57      DCAVIPLKRHKDYL--LIQTVDFFYPMVNDPELLGRIALANVLSDVYAVGVTQFDTVEMI 114
               D AV+     KD    ++ TVDFF   V+DP + G+IA  + LSD  A+G      + ++
Sbjct: 1819235 DGAVVAYGPEKDQRTQVVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALV 1819056

Query: 115     VSTSTSFSEKERDVVIGLVMKGFQNSLKANGYRNTPLIIRQLKINPWCIIGGIATSVCRS 174
               V    S +    D +  L M G + +L A G                 ++GG     C  
Sbjct: 1819055 VVPYGS-NRCVEDTLFQL-MSGARVALDAAGCT---------------LVGG---HTCEG 1818936

Query: 175     EEIIL--PSNAQP--------------GDVLVLTKPLGGQMAMDAHLWQLNQTEKYKKLL 218
                E+ L    N  P              GD L+LTKPLG        +  L++  + ++ L
Sbjct: 1818935 AELGLGFAVNGVPDGGAAFLSKEGMVAGDALILTKPLGTGAVFAGDMRALSRAAEVRECL 1818756

Query: 219     SECSDADIKETFEIAVKSMTYLNKNAALLMHKYQAHCATDITGFGLLGHANNLAQFQ-KE 277
               +                SMT  N+ AA  + ++     TD+TGFGLLGH   + +     
Sbjct: 1818755 A----------------SMTQSNRAAATELKEHGCTSCTDVTGFGLLGHLAEMCRASGGA 1818624

Query: 278     KVLFQINKLPIIKNVLK------FSTLVGQSTKFRSGRSVE------------------- 312
               +V   +   P+     +      FS+L  Q    R+ R+VE                   
Sbjct: 1818623 RVELSLAATPVFAGAKRLARDGVFSSL--QDENLRAKRAVENHGEAVELDAAAYALLFDP 1818450

Query: 313     -TSGGLLICLPADAAD 327
                T+GGLL  +PAD AD
Sbjct: 1818449 QTAGGLLATVPADRAD 1818402


>scaffold_1000
          Length = 1887

 Score = 32.3 bits (72), Expect = 1.0
 Identities = 13/41 (31%), Positives = 22/41 (53%)
 Frame = +3

Query: 166  GIATSVCRSEEIILPSNAQPGDVLVLTKPLGGQMAMDAHLW 206
            G  T   +   +++P NA+PGD L +  P  G+ A++   W
Sbjct: 1629 GSRTEKQKIVRVLVPKNAKPGDSLAVATPTTGKFAVEVPTW 1751


>scaffold_2
          Length = 3465789

 Score = 29.3 bits (64), Expect = 8.6
 Identities = 38/181 (20%), Positives = 68/181 (37%), Gaps = 24/181 (13%)
 Frame = +2

Query: 208     LNQTEKYKKLLSECSDADIKETFEIAVKSMTYL-----NKNAALLMHKYQAHCATDITGF 262
               L   E  + +  E  +A  +   ++  +  TYL     +K    +M K  A  A  +   
Sbjct: 1214933 LGLVELVENVEQEADEATGRMPTDVLAQGKTYLLVLSQDKYTTEMMAKRDADVACGL--- 1215103

Query: 263     GLLGHANNLAQFQKEKVLFQINKLPIIK-----------NVLKFSTLVGQSTKFRSGRSV 311
                     +NL    K+KVL  +NK+ +++           +V  F +L G     + G   
Sbjct: 1215104 ----FCDNLVNLTKDKVLAWLNKVSVLQSLGEPALRKLADVAHFRSLPGGKVVLKQGDRA 1215271

Query: 312     E-----TSGGLLICLPADAADKFCRDFEEATNGEQKSFQIGHVTAAN---ESDAVLCEDV 363
               E      SG L +   +D         +  T+  +K   +G++   +   ES  +  ED 
Sbjct: 1215272 ERLHVILSGSLDVACKSDMRKTIVGQRDSTTHDGKKQVVVGNLHEGDFFGESALLSTEDS 1215451

Query: 364     E 364
               E
Sbjct: 1215452 E 1215454


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
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: 17,538,640
Number of Sequences: 1185
Number of extensions: 205827
Number of successful extensions: 868
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 36
Number of HSP's successfully gapped in prelim test: 18
Number of HSP's that attempted gapping in prelim test: 669
Number of HSP's gapped (non-prelim): 444
length of query: 370
length of database: 18,886,859
effective HSP length: 106
effective length of query: 264
effective length of database: 18,761,249
effective search space: 4952969736
effective search space used: 4952969736
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: 64 (29.3 bits)
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= SPP00000055_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Anopheles gambiae # Complete
         (396 letters)

Database: genome.fa 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             86   1e-16
scaffold_9                                                             35   0.13 
scaffold_25                                                            33   0.85 
scaffold_15                                                            30   4.2  
scaffold_400                                                           30   5.5  
scaffold_12                                                            30   5.5  

>scaffold_6
          Length = 2755913

 Score = 85.5 bits (210), Expect = 1e-16
 Identities = 85/303 (28%), Positives = 132/303 (43%), Gaps = 30/303 (9%)
 Frame = -1

Query: 102     MVQTTDFFYPIVDDPYMMGKIACANVLSDLYAMGVTEVDNMLMLLAV--STKMIEKERDV 159
               +V T DFF   V DPY+ G+IA  + LSD  AMG      + +++    S + +E   D 
Sbjct: 1819187 VVHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALVVVPYGSNRCVE---DT 1819017

Query: 160     VIPLIMRGFKDSALEAGTSVTGGHSVVNPWCTIGGVATTICQQN-EFIVPDNAVVGDVLV 218
               +  L M G + +   AG ++ GGH+       +G     +      F+  +  V GD L+
Sbjct: 1819016 LFQL-MSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALI 1818840

Query: 219     LTKALGTQVAVNAHQWLDQSERWNRIKLVVSEEDVRKAYHRAMDSMSRLNRVAARLMHKY 278
               LTK LGT  AV A            ++ +    +VR+     + SM++ NR AA  + ++
Sbjct: 1818839 LTKPLGTG-AVFA----------GDMRALSRAAEVRE----CLASMTQSNRAAATELKEH 1818705

Query: 279     NAHGATDVTGFGLLGH-AQTLASHQKNEVSFVIHNLPVIAKMAAVAKACGNMFQLLQGHS 337
                    TDVTGFGLLGH A+   +     V   +   PV A    +A+    +F  LQ  +
Sbjct: 1818704 GCTSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFAGAKRLAR--DGVFSSLQDEN 1818531

Query: 338     --------------------------AETSGGLLICLPREQAAAYCKDIEKQEGCQAWII 371
                                          +T+GGLL  +P ++A A    I +    +A +I
Sbjct: 1818530 LRAKRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIRRVGCARAAVI 1818351

Query: 372     GIV 374
               G V
Sbjct: 1818350 GTV 1818342


>scaffold_9
          Length = 1602328

 Score = 35.4 bits (80), Expect = 0.13
 Identities = 25/64 (39%), Positives = 32/64 (50%), Gaps = 1/64 (1%)
 Frame = +1

Query: 258    HRAMDSMSRLNRVAARLMHKYNAHGATDVTGFGLLGHAQTLASHQKNEVSFVIHN-LPVI 316
              HR    + RL R  AR      A GA  VTG G LGHA+    H +  V  VI + +P  
Sbjct: 652615 HRRQQLLPRLGRQGARDFCG-EAQGAAGVTGHGALGHARDGRGHGRGRVGHVIFSRIPRG 652791

Query: 317    AKMA 320
              A++A
Sbjct: 652792 AELA 652803


>scaffold_25
          Length = 552203

 Score = 32.7 bits (73), Expect = 0.85
 Identities = 23/74 (31%), Positives = 34/74 (45%), Gaps = 3/74 (4%)
 Frame = -2

Query: 224    GTQVAVNAHQWLDQSE-RWNRIKLVVSEEDVRKAYHRAMDSMSRLNRVAARLMHKYNAH- 281
              G + A   H + D++  R  R+      E  R    R +D++ R  R   RL  + +A  
Sbjct: 309094 GRRRAAQGHDFPDEAAPRGLRLLGAPQREQRRNVLRRPVDALPRFGRRPLRLRPRASAPR 308915

Query: 282    -GATDVTGFGLLGH 294
               G  DV GFG +GH
Sbjct: 308914 AGQDDVRGFGAVGH 308873


>scaffold_15
          Length = 987846

 Score = 30.4 bits (67), Expect = 4.2
 Identities = 19/55 (34%), Positives = 28/55 (50%)
 Frame = +1

Query: 342    GGLLICLPREQAAAYCKDIEKQEGCQAWIIGIVEKGSRTARIIDKPRVIEVPAKD 396
              G  ++  P    A   +D+E   G    ++ +V+ G R ARI   PRV+ VPA D
Sbjct: 479773 GAGVVVDPGHVEARVPRDVEVGAGMMVVMMIVVDPGRREARI---PRVLRVPAAD 479928


>scaffold_400
          Length = 7388

 Score = 30.0 bits (66), Expect = 5.5
 Identities = 16/42 (38%), Positives = 21/42 (50%)
 Frame = -1

Query: 338 AETSGGLLICLPREQAAAYCKDIEKQEGCQAWIIGIVEKGSR 379
           AE  G  +   P  +  A CKD EKQE    W+   +E+ SR
Sbjct: 791 AERDGVQISIYPCPKITASCKDEEKQEEQDQWLHDALERESR 666


>scaffold_12
          Length = 1412611

 Score = 30.0 bits (66), Expect = 5.5
 Identities = 14/34 (41%), Positives = 19/34 (55%)
 Frame = +1

Query: 88      LDCSVIPLRHGGLCMVQTTDFFYPIVDDPYMMGK 121
               LD  V P+RHGG+ +V   D   P+VD     G+
Sbjct: 1329898 LDDEVAPVRHGGVVVVLADDEVAPVVDGRRRRGR 1329999


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
Lambda     K      H
   0.322    0.136    0.410 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 25,981,269
Number of Sequences: 1185
Number of extensions: 380853
Number of successful extensions: 1505
Number of sequences better than 10.0: 6
Number of HSP's better than 10.0 without gapping: 45
Number of HSP's successfully gapped in prelim test: 30
Number of HSP's that attempted gapping in prelim test: 1007
Number of HSP's gapped (non-prelim): 865
length of query: 396
length of database: 18,886,859
effective HSP length: 106
effective length of query: 290
effective length of database: 18,761,249
effective search space: 5440762210
effective search space used: 5440762210
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.9 bits)
S2: 64 (29.3 bits)
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 
           1185 sequences; 56,660,578 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

scaffold_6                                                             92   1e-18
scaffold_3                                                             32   1.0  
scaffold_1                                                             31   2.2  
scaffold_7                                                             30   3.8  
scaffold_2                                                             30   3.8  
scaffold_60                                                            30   6.5  
scaffold_5                                                             29   8.5  

>scaffold_6
          Length = 2755913

 Score = 92.0 bits (227), Expect = 1e-18
 Identities = 99/360 (27%), Positives = 156/360 (43%), Gaps = 32/360 (8%)
 Frame = -1

Query: 23      GXGSKVPQDVLNRLLAGVYGEQLGDKDGKGPNKEEGVGIGL--DSSVIALKHD----LFL 76
               G G+KV    L R LA +  + L   +G  P  +  VG+    D +V+A   +      +
Sbjct: 1819361 GCGAKVGASALTRALANL-PKSLTFDNGTPPPCDVVVGVDAPDDGAVVAYGPEKDQRTQV 1819185

Query: 77      VQSVDFFYPLIDDPFMLGKIALANVVSDVFAVGATEIDQIKLIVTAPTEFTEQEREVVVP 136
               V +VDFF   + DP++ G+IA  + +SD  A+GA     + L+V      + +  E  + 
Sbjct: 1819184 VHTVDFFRSFVSDPYVFGQIAANHALSDCDAMGALPASALALVVVPYG--SNRCVEDTLF 1819011

Query: 137     MVMQGFLEAAKACNAPVQIGSIAENPWCVIGGAASAVCHRSELIMPYNAQ-PGDALVLTK 195
                +M G   A  A    +  G   E     +G A + V       +       GDAL+LTK
Sbjct: 1819010 QLMSGARVALDAAGCTLVGGHTCEGAELGLGFAVNGVPDGGAAFLSKEGMVAGDALILTK 1818831

Query: 196     PLGTQLATNAYIWMGEQQSDSWAQLRERFTVADIEQTYRIALESMSRLNKTGAELMKKYG 255
               PLGT       ++ G+ ++ S A        A++ +     L SM++ N+  A  +K++G
Sbjct: 1818830 PLGT-----GAVFAGDMRALSRA--------AEVRE----CLASMTQSNRAAATELKEHG 1818702

Query: 256     AHAATDVTGFGLYGH-AENLASHQTADVDFHLDTLPIIKNVREIAE----TLGRGAKLLA 310
                 + TDVTGFGL GH AE   +   A V+  L   P+    + +A     +  +   L A
Sbjct: 1818701 CTSCTDVTGFGLLGHLAEMCRASGGARVELSLAATPVFAGAKRLARDGVFSSLQDENLRA 1818522

Query: 311     GKAVE--------------------TSGGLLICLPREAAAGFCEEYRRCTKHEAWIVGHV 350
                +AVE                    T+GGLL  +P + A       RR     A ++G V
Sbjct: 1818521 KRAVENHGEAVELDAAAYALLFDPQTAGGLLATVPADRADAALAAIRRVGCARAAVIGTV 1818342


>scaffold_3
          Length = 3001400

 Score = 32.3 bits (72), Expect = 1.0
 Identities = 18/45 (40%), Positives = 27/45 (60%), Gaps = 3/45 (6%)
 Frame = -2

Query: 33      LNRLLAGVYGEQLGDKDGKG---PNKEEGVGIGLDSSVIALKHDL 74
               L+R+L      +L D + +G   P +  GVG+GLD  V+AL+ DL
Sbjct: 1942105 LHRVLLVALEHRLVDLEERGVDVPREVAGVGVGLDLEVLALRRDL 1941971


>scaffold_1
          Length = 3520325

 Score = 31.2 bits (69), Expect = 2.2
 Identities = 17/55 (30%), Positives = 29/55 (52%)
 Frame = -3

Query: 27      KVPQDVLNRLLAGVYGEQLGDKDGKGPNKEEGVGIGLDSSVIALKHDLFLVQSVD 81
               +V ++ L  LL  V G  + D+   G +  E VG+G  ++ +  +HD  L + VD
Sbjct: 1443480 RVDREDLVALLEAVRGLAVDDRREAGDDAREDVGLGHAAAAVRQRHDARLAELVD 1443316


>scaffold_7
          Length = 2177197

 Score = 30.4 bits (67), Expect = 3.8
 Identities = 23/84 (27%), Positives = 39/84 (46%)
 Frame = +3

Query: 224    FTVADIEQTYRIALESMSRLNKTGAELMKKYGAHAATDVTGFGLYGHAENLASHQTADVD 283
              F +A +EQ   I++E         AE+  K+      D+  +      +   +HQ+AD D
Sbjct: 225354 FLLAPVEQPGYISVEEFVIFGLVVAEIEPKFH-QVKEDLAAYRAASGTQRW-THQSADRD 225527

Query: 284    FHLDTLPIIKNVREIAETLGRGAK 307
                ++TLP   + R IA+  G G +
Sbjct: 225528 PSIETLPKDGSQRAIAQLAGTGTR 225599


>scaffold_2
          Length = 3465789

 Score = 30.4 bits (67), Expect = 3.8
 Identities = 21/73 (28%), Positives = 31/73 (42%), Gaps = 15/73 (20%)
 Frame = +2

Query: 163    WCVIGGAA-SAVCHRSELIMPYNAQ-----------PGDALVLTK---PLGTQLATNAYI 207
              WC + GA+  A   R  ++    A            P DAL + +   P+G  LAT   +
Sbjct: 919091 WCALPGASLEAAAGRDGVVAKVEASLSGAVDLSTASPRDALAVLQCAVPVGDDLATLGGV 919270

Query: 208    WMGEQQSDSWAQL 220
              W  +   D+WA L
Sbjct: 919271 WSPKAARDAWAPL 919309


>scaffold_60
          Length = 140414

 Score = 29.6 bits (65), Expect = 6.5
 Identities = 22/84 (26%), Positives = 37/84 (44%)
 Frame = +3

Query: 224   FTVADIEQTYRIALESMSRLNKTGAELMKKYGAHAATDVTGFGLYGHAENLASHQTADVD 283
             F +A +EQ   I++E         AE+  K+  H   +              +HQ+AD D
Sbjct: 22458 FLLAPVEQPGYISVEEFVIFGLVVAEIEPKF--HQVKEDLAAHRAASGTQRWTHQSADRD 22631

Query: 284   FHLDTLPIIKNVREIAETLGRGAK 307
               ++TLP   + R +A+  G G +
Sbjct: 22632 PSIETLPKDGSQRAMAQLAGTGTR 22703


>scaffold_5
          Length = 3000003

 Score = 29.3 bits (64), Expect = 8.5
 Identities = 23/67 (34%), Positives = 37/67 (55%), Gaps = 9/67 (13%)
 Frame = -3

Query: 251    MKKYGAHAATDVTGF-GLYGH--AENLASHQTA----DVDFH--LDTLPIIKNVREIAET 301
              ++ +GAH    V  F G YG    E+ A+ + A    D  F+  LD++P++    ++AE 
Sbjct: 391093 IRDHGAHGVVAVLTFAGGYGPFPLEDAATVECAGGAPDDAFNRVLDSIPVVAVHEQVAEA 390914

Query: 302    LGRGAKL 308
              LGRG K+
Sbjct: 390913 LGRGDKV 390893


  Database: genome.fa
    Posted date:  Feb 16, 2009  1:41 PM
  Number of letters in database: 56,660,578
  Number of sequences in database:  1185
  
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: 22,698,897
Number of Sequences: 1185
Number of extensions: 346464
Number of successful extensions: 2311
Number of sequences better than 10.0: 8
Number of HSP's better than 10.0 without gapping: 112
Number of HSP's successfully gapped in prelim test: 40
Number of HSP's that attempted gapping in prelim test: 949
Number of HSP's gapped (non-prelim): 1855
length of query: 369
length of database: 18,886,859
effective HSP length: 106
effective length of query: 263
effective length of database: 18,761,249
effective search space: 4934208487
effective search space used: 4934208487
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)
