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= gi|12408011|gb|AAG53689.1|AF322072_1 selenoprotein Pb [Danio
rerio]
         (265 letters)

Database: /cursos/BI/genomes/project_2013/Gadus_morhua/genome.fa 
           398,859 sequences; 832,114,588 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

GeneScaffold_400 dna:genescaffold genescaffold:gadMor1:GeneScaff...    82   9e-27
contig56004 dna:contig contig::contig56004:1:2727:1 REF                72   6e-12

>GeneScaffold_400 dna:genescaffold genescaffold:gadMor1:GeneScaffold_400:1:166366:1 REF
          Length = 166366

 Score = 81.6 bits (200), Expect(2) = 9e-27
 Identities = 55/145 (37%), Positives = 70/145 (48%), Gaps = 56/145 (38%)
 Frame = +2

Query: 23    VEKESNGSRICKPAPQWEIDGKTPMKELLGNVVVVALLKASXHFCLTQAARLGDLRDKLA 82
             VE +S+ SRICKPAP WE+ G  PM ELLGNVVVVALLKAS  FCLTQA++  ++ + + 
Sbjct: 16319 VEGDSDPSRICKPAPHWEVGGIAPMPELLGNVVVVALLKAS*QFCLTQASK*AEIMNDVM 16498

Query: 83    N----GGLTNISFMVV-------------------------------------------- 94
             +      L+N+S  +V                                            
Sbjct: 16499 HF*GFVFLSNVSKPIVKQSPLRYCFVKSQKQLIRQIPL*LLL*LFLRIGHLREKLSRSNL 16678

Query: 95    --------NEQDSQSRAMYWELKRR 111
                     NE++ QSRAMYWELKRR
Sbjct: 16679 TDVAFMIVNEREFQSRAMYWELKRR 16753



 Score = 62.8 bits (151), Expect = 5e-09
 Identities = 26/39 (66%), Positives = 34/39 (87%)
 Frame = +3

Query: 144   RCGYLTFHIVLPFSFLHYPYIEAAIRATYHKNMCNCSLN 182
             RCG LTFHIVLP+SFL  PY+EAA+RAT  K++CNC+++
Sbjct: 17178 RCGLLTFHIVLPYSFLTQPYVEAAVRATSLKDICNCTVS 17294



 Score = 60.8 bits (146), Expect(2) = 9e-27
 Identities = 25/30 (83%), Positives = 29/30 (96%)
 Frame = +3

Query: 115   DIPVYQQSPLQNDVWEILEGDKDDFLVYDR 144
             ++PVYQQ PLQ+DVWEIL+GDKDDFLVYDR
Sbjct: 16764 NVPVYQQGPLQDDVWEILDGDKDDFLVYDR 16853


>contig56004 dna:contig contig::contig56004:1:2727:1 REF
          Length = 2727

 Score = 72.4 bits (176), Expect = 6e-12
 Identities = 32/70 (45%), Positives = 47/70 (67%)
 Frame = +1

Query: 77  LRDKLANGGLTNISFMVVNEQDSQSRAMYWELKRRTAQDIPVYQQSPLQNDVWEILEGDK 136
           LR KL N GL N+++MV+N Q  Q++ ++  L  R +  I +Y+Q     DVW+ L+G+K
Sbjct: 46  LRQKLENQGLVNVNYMVINHQGIQAQRLHSMLTTRLSSSISLYKQGDQDADVWQALKGEK 225

Query: 137 DDFLVYDRCG 146
           DDFL+YDR G
Sbjct: 226 DDFLIYDRSG 255


  Database: /cursos/BI/genomes/project_2013/Gadus_morhua/genome.fa
    Posted date:  Oct 24, 2012  8:34 PM
  Number of letters in database: 832,114,588
  Number of sequences in database:  398,859
  
Lambda     K      H
   0.315    0.131    0.402 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 304,726,333
Number of Sequences: 398859
Number of extensions: 4799658
Number of successful extensions: 77595
Number of sequences better than 1.0e-04: 2
Number of HSP's better than  0.0 without gapping: 38453
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 0
Number of HSP's gapped (non-prelim): 74990
length of query: 265
length of database: 277,371,529
effective HSP length: 119
effective length of query: 146
effective length of database: 229,907,308
effective search space: 33566466968
effective search space used: 33566466968
frameshift window, decay const: 40,  0.1
T: 13
A: 40
X1: 16 ( 7.3 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.6 bits)
S2: 114 (48.5 bits)
