TBLASTN 2.2.17 [Aug-26-2007]


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.

Reference for composition-based statistics:
Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
Eugene V. Koonin, and Stephen F. Altschul (2001), 
"Improving the accuracy of PSI-BLAST protein database searches with 
composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

Query= SPP00000033_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Homo sapiens # Complete
         (448 letters)

Database: genome.fa 
           23,491 sequences; 86,045,674 total letters

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



                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gb|AAXJ01000359.1| Perkinsus marinus ATCC 50983 strain ATCC 5098...   238   2e-62
gb|AAXJ01002698.1| Perkinsus marinus ATCC 50983 strain ATCC 5098...    31   3.8  
gb|AAXJ01000431.1| Perkinsus marinus ATCC 50983 strain ATCC 5098...    30   7.4  

>gb|AAXJ01000359.1| Perkinsus marinus ATCC 50983 strain ATCC 50983 gcontig_1104296309226,
             whole genome shotgun sequence
          Length = 32957

 Score =  238 bits (607), Expect = 2e-62,   Method: Composition-based stats.
 Identities = 162/416 (38%), Positives = 231/416 (55%), Gaps = 12/416 (2%)
 Frame = -3

Query: 36    RPFEPQALGLSPSWRLTGFSGMKGUGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEE 95
             R ++P +LGLS  + LT FS  KG GCKV  + L   L         P +G G +  ++ 
Sbjct: 29805 RVWDPTSLGLSKDFVLTSFSERKG*GCKVQPDQLSSFL---------PQIGDGNMLTED- 29656

Query: 96    ASQEAGLPAGAGPSPTFPALGIGMDSCVIPLRHGGLSLVQTTDFFYPLVEDPYMMGRIAC 155
                          SP         D  ++P+  G   LV +TDF  P+ +DP+  G IA 
Sbjct: 29655 -------------SP---------DCSMVPVGDGERVLVSSTDFLSPISQDPFHQGWIAA 29542

Query: 156   ANVLSDLYAMGITECDNMLMLLSVSQSMSEEEREKVTPLMVKGFRDAAEEGGTAVTGGQT 215
             AN LSDLYA+G  + +N+LM+L+VS  MS  E+   T  M++GF    EE GT VTGGQT
Sbjct: 29541 ANTLSDLYAVGADDINNILMILTVSMKMSPLEQTVTTKKMIEGFNACCEEAGTKVTGGQT 29362

Query: 216   VVNPWIIIGGVATVVCQPNEFIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDNPERWNK- 274
             V NPW ++GG A      ++ +  D AV G+V+VLTKPLGTQV  N   W+ N  +W K 
Sbjct: 29361 VFNPWPLLGGCAMATMLQSQLVTTDRAVPGNVIVLTKPLGTQVTGNLALWMMNEAQWGKV 29182

Query: 275   VKMVVSREEVELAYQEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQR 334
             VK  +  E    A +    +M TLN+    LM  F A  ATD+TGFG++GH+ NLA+ Q+
Sbjct: 29181 VKAGLDVETCTRAVKMGDLSMMTLNKNGCRLMREFGALGATDVTGFGLVGHADNLARVQK 29002

Query: 335   NEVSFVIHNLPII--AKMAAVSKASGR-FGLLQGTSAETSGGLLICL-PREQAARFCSEI 390
             +  S  I  +P I   K  ++ + +GR +GLL+G +AETSGGLL+ L  +E A +F + +
Sbjct: 29001 SGASIEIDTVPTIQGCKTLSMVRPNGRDYGLLEGRAAETSGGLLVMLRDKETAEKFIAAL 28822

Query: 391   KSSKYGEGHQ---AWIVGIV---EKGNRT-ARIIDKPRVIEVLPRGATAAVLAPDS 439
             ++    +  +    WIVG V   +KG++  A+++  P VIEV   GA A  +A  S
Sbjct: 28821 RACYSDDDPRREWGWIVGRVVETKKGDKQGAKVLPDPNVIEV--*GARAVSIAARS 28660


>gb|AAXJ01002698.1| Perkinsus marinus ATCC 50983 strain ATCC 50983 gcontig_1104296688020,
            whole genome shotgun sequence
          Length = 5752

 Score = 31.2 bits (69), Expect = 3.8,   Method: Composition-based stats.
 Identities = 18/50 (36%), Positives = 28/50 (56%), Gaps = 8/50 (16%)
 Frame = +2

Query: 241  SAVVGDVLVLTKPLGTQV--------AVNAHQWLDNPERWNKVKMVVSRE 282
            SA++G VL+ T+PL T +         +N+H+W  N  R+ K +   SRE
Sbjct: 3287 SALIG-VLLTTQPLNTAMLYVYMYM*VLNSHEWGQNLSRYRKSESTSSRE 3433


>gb|AAXJ01000431.1| Perkinsus marinus ATCC 50983 strain ATCC 50983 gcontig_1104296665852,
            whole genome shotgun sequence
          Length = 29653

 Score = 30.4 bits (67), Expect = 7.4,   Method: Composition-based stats.
 Identities = 22/67 (32%), Positives = 30/67 (44%)
 Frame = -3

Query: 45   LSPSWRLTGFSGMKGUGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPA 104
            L P+   +G S  +  GC    EAL  L   ++R   +PPL  GL         EA +P 
Sbjct: 5294 LPPTLISSGTSQRRPAGCTTGTEALRSLRQNISRHCGKPPLPPGLTVSAHITVYEA-MPQ 5118

Query: 105  GAGPSPT 111
             + P PT
Sbjct: 5117 HSIPRPT 5097


  Database: genome.fa
    Posted date:  Jan 16, 2009 12:51 PM
  Number of letters in database: 86,045,674
  Number of sequences in database:  23,491
  
Lambda     K      H
   0.317    0.133    0.392 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Sequences: 23491
Number of Hits to DB: 42,262,485
Number of extensions: 625437
Number of successful extensions: 2738
Number of sequences better than 10.0: 5
Number of HSP's gapped: 2736
Number of HSP's successfully gapped: 5
Length of query: 448
Length of database: 28,681,891
Length adjustment: 110
Effective length of query: 338
Effective length of database: 26,097,881
Effective search space: 8821083778
Effective search space used: 8821083778
Neighboring words threshold: 13
Window for multiple hits: 40
X1: 16 ( 7.3 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.7 bits)
S2: 43 (21.2 bits)
