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= SPP00000001_1.0 # Protein # Iodothyronine deiodinase 1 (DI1) #
Homo sapiens # Complete
         (249 letters)

Database: genome.fa 
           113 sequences; 29,453,142 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gb|AAFD02000012.1| Thalassiosira pseudonana CCMP1335 chromosome ...    30   2.0  
gb|AAFD02000041.1| Thalassiosira pseudonana CCMP1335 chromosome ...    28   5.7  
gb|AAFD02000005.1| Thalassiosira pseudonana CCMP1335 chromosome ...    27   9.8  

>gb|AAFD02000012.1| Thalassiosira pseudonana CCMP1335 chromosome 9 THAPSchr_9_Cont12, whole
               genome shotgun sequence
          Length = 1191060

 Score = 29.6 bits (65), Expect = 2.0
 Identities = 11/43 (25%), Positives = 21/43 (48%)
 Frame = +2

Query: 195     PVVVDTMQNQSSQLYAALPERLYIIQEGRILYKGKSGPWNYNP 237
               P  V+ +    + +Y       YI++  R+ Y G +  W+Y+P
Sbjct: 1178762 PPTVNVLHIAQTSIYDGSESVAYILELARLTYGGSASSWHYSP 1178890


>gb|AAFD02000041.1| Thalassiosira pseudonana CCMP1335 chromosome 22 THAPSchr_22_Cont41,
            whole genome shotgun sequence
          Length = 263482

 Score = 28.1 bits (61), Expect = 5.7
 Identities = 11/34 (32%), Positives = 17/34 (50%)
 Frame = -2

Query: 97   VVRLSGQRCNIWEFMQGNRPLVLNFGSCTUPSFM 130
            ++  +  RC I+   QGN   +   G+ T P FM
Sbjct: 6765 LISFASNRCTIYVETQGNSISIFKLGAATVPDFM 6664


>gb|AAFD02000005.1| Thalassiosira pseudonana CCMP1335 chromosome 3 THAPSchr_3_Cont5, whole
               genome shotgun sequence
          Length = 2440052

 Score = 27.3 bits (59), Expect = 9.8
 Identities = 17/48 (35%), Positives = 23/48 (47%), Gaps = 1/48 (2%)
 Frame = -3

Query: 72      WFVLKVRWQRLEDTTELGGLAPNCPVVRLSGQRCNIW-EFMQGNRPLV 118
               WF+ K     + D  EL G   NC  V +SG    +W   M+  RP+V
Sbjct: 2381583 WFLSK----DVRDAAELFGEQMNCSFVSMSGLNSTLWCVMMR*MRPIV 2381452



 Score = 27.3 bits (59), Expect = 9.8
 Identities = 17/48 (35%), Positives = 23/48 (47%), Gaps = 1/48 (2%)
 Frame = -2

Query: 72      WFVLKVRWQRLEDTTELGGLAPNCPVVRLSGQRCNIW-EFMQGNRPLV 118
               WF+ K     + D  EL G   NC  V +SG    +W   M+  RP+V
Sbjct: 2376625 WFLSK----DVRDAAELFGEQMNCSFVSMSGLNSTLWCVMMR*MRPIV 2376494


  Database: genome.fa
    Posted date:  Jan 19, 2010  5:14 PM
  Number of letters in database: 29,453,142
  Number of sequences in database:  113
  
Lambda     K      H
   0.324    0.138    0.448 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 9,220,744
Number of Sequences: 113
Number of extensions: 134799
Number of successful extensions: 788
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 5
Number of HSP's successfully gapped in prelim test: 6
Number of HSP's that attempted gapping in prelim test: 726
Number of HSP's gapped (non-prelim): 178
length of query: 249
length of database: 9,817,714
effective HSP length: 97
effective length of query: 152
effective length of database: 9,806,753
effective search space: 1490626456
effective search space used: 1490626456
frameshift window, decay const: 40,  0.1
T: 13
A: 40
X1: 15 ( 7.0 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 40 (21.5 bits)
S2: 59 (27.3 bits)