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= SPP00000014_1.0 # Protein # Selenoprotein H (SelH) # Homo
sapiens # Complete
         (122 letters)

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

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



                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gb|AAXJ01000722.1| Perkinsus marinus ATCC 50983 strain ATCC 5098...    41   4e-04
gb|AAXJ01000644.1| Perkinsus marinus ATCC 50983 strain ATCC 5098...    41   5e-04

>gb|AAXJ01000722.1| Perkinsus marinus ATCC 50983 strain ATCC 50983 gcontig_1104296790916,
            whole genome shotgun sequence
          Length = 20479

 Score = 41.2 bits (95), Expect = 4e-04,   Method: Composition-based stats.
 Identities = 25/56 (44%), Positives = 33/56 (58%), Gaps = 3/56 (5%)
 Frame = +3

Query: 32   EEATVVIEHCTSURVYGRNAAALSQALRLEAPELPVKVNP---TKPRRGSFEVTLL 84
            EEA V+IE C S  ++ R+AAAL  AL    P   V VN     KPR+G+F V ++
Sbjct: 7269 EEADVLIEACKS*GIFKRSAAALGVALEARRPGTTVVVNENGLAKPRKGTFAVKVI 7436


>gb|AAXJ01000644.1| Perkinsus marinus ATCC 50983 strain ATCC 50983 gcontig_1104296642754,
             whole genome shotgun sequence
          Length = 22778

 Score = 40.8 bits (94), Expect = 5e-04,   Method: Composition-based stats.
 Identities = 25/56 (44%), Positives = 33/56 (58%), Gaps = 3/56 (5%)
 Frame = +3

Query: 32    EEATVVIEHCTSURVYGRNAAALSQALRLEAPELPVKVNP---TKPRRGSFEVTLL 84
             EEA V+IE C S  ++ R+AAAL  AL    P   V VN     KPR+G+F V ++
Sbjct: 18159 EEADVLIEACKS*GIFKRSAAALGVALESRRPGTTVVVNENGLAKPRKGTFAVKVI 18326


  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.314    0.130    0.369 

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: 10,788,303
Number of extensions: 139490
Number of successful extensions: 516
Number of sequences better than 10.0: 8
Number of HSP's gapped: 516
Number of HSP's successfully gapped: 8
Length of query: 122
Length of database: 28,681,891
Length adjustment: 92
Effective length of query: 30
Effective length of database: 26,520,719
Effective search space: 795621570
Effective search space used: 795621570
Neighboring words threshold: 13
Window for multiple hits: 40
X1: 16 ( 7.2 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 42 (21.9 bits)
S2: 38 (19.2 bits)
