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: /users/rg/didac/GENOMES/C.fasciculata/genome.fa 
           1089 sequences; 37,048,325 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

Cf_Contig930 |  | 1 to 140823                                          30   2.4  
Cf_Contig1045 |  | 1 to 270110                                         29   3.2  
Cf_Contig551 |  | 1 to 16524                                           29   3.2  
Cf_Contig676 |  | 1 to 122494                                          28   5.4  
Cf_Contig888 |  | 1 to 77865                                           28   9.2  

>Cf_Contig930 |  | 1 to 140823
          Length = 140823

 Score = 29.6 bits (65), Expect = 2.4
 Identities = 11/38 (28%), Positives = 21/38 (55%)
 Frame = -3

Query: 117   LVLNFGSCTUPSFMFKFDQFKRLIEDFSSIADFLVIYI 154
             L ++ G C  PS  F+F   +++   + S   FL++Y+
Sbjct: 85426 LCVSLGMCKAPSLPFRFSHVEKMALPYISFFSFLLLYV 85313


>Cf_Contig1045 |  | 1 to 270110
          Length = 270110

 Score = 29.3 bits (64), Expect = 3.2
 Identities = 16/54 (29%), Positives = 23/54 (42%)
 Frame = -3

Query: 76     KVRWQRLEDTTELGGLAPNCPVVRLSGQRCNIWEFMQGNRPLVLNFGSCTUPSF 129
              K +  +LE T    GL  NC        +C IW ++  +   +L    CT P F
Sbjct: 171855 KFKEAQLELTVTASGLCVNCSTATKIDDQCAIWRYVDTSASWLL-VAFCTGPHF 171697


>Cf_Contig551 |  | 1 to 16524
          Length = 16524

 Score = 29.3 bits (64), Expect = 3.2
 Identities = 18/48 (37%), Positives = 26/48 (54%), Gaps = 1/48 (2%)
 Frame = +3

Query: 187   LLARSPQCPVVVDTMQNQSSQLY-AALPERLYIIQEGRILYKGKSGPW 233
             LLA  P CP      +    Q Y AA   +L + +EG+ILY+G+  P+
Sbjct: 13797 LLAVRPGCPAHRGASERGGPQAYRAAGCTQLKLQEEGQILYQGQEWPF 13940


>Cf_Contig676 |  | 1 to 122494
          Length = 122494

 Score = 28.5 bits (62), Expect = 5.4
 Identities = 13/29 (44%), Positives = 18/29 (62%), Gaps = 7/29 (24%)
 Frame = +2

Query: 54     PHFS-------HDNWIPTFFSTQYFWFVL 75
              PHFS       H N+I +FF+T +F+ VL
Sbjct: 107870 PHFSLFLFF*PHQNFISSFFTTSFFFLVL 107956


>Cf_Contig888 |  | 1 to 77865
          Length = 77865

 Score = 27.7 bits (60), Expect = 9.2
 Identities = 15/44 (34%), Positives = 22/44 (50%)
 Frame = +1

Query: 191   SPQCPVVVDTMQNQSSQLYAALPERLYIIQEGRILYKGKSGPWN 234
             S   P V+ T  + +S  YA LP R Y  + G  +  G+ G W+
Sbjct: 50671 SESFPSVICTTGSGASNFYA-LPPRYYADRSGWFIKMGRGGDWH 50799


  Database: /users/rg/didac/GENOMES/C.fasciculata/genome.fa
    Posted date:  Sep 29, 2011  8:09 PM
  Number of letters in database: 37,048,325
  Number of sequences in database:  1089
  
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: 11,763,721
Number of Sequences: 1089
Number of extensions: 170556
Number of successful extensions: 944
Number of sequences better than 10.0: 6
Number of HSP's better than 10.0 without gapping: 61
Number of HSP's successfully gapped in prelim test: 37
Number of HSP's that attempted gapping in prelim test: 829
Number of HSP's gapped (non-prelim): 270
length of query: 249
length of database: 12,349,441
effective HSP length: 99
effective length of query: 150
effective length of database: 12,241,630
effective search space: 1836244500
effective search space used: 1836244500
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: 60 (27.7 bits)