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|379656510|gb|AFD09496.1| glutathione peroxidase 4a
[Oncorhynchus kisutch]
         (200 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

contig411933 dna:contig contig::contig411933:1:1296:1 REF              96   4e-19
GeneScaffold_3230 dna:genescaffold genescaffold:gadMor1:GeneScaf...    87   1e-16
GeneScaffold_1266 dna:genescaffold genescaffold:gadMor1:GeneScaf...    81   1e-14
contig311944 dna:contig contig::contig311944:1:640:1 REF               73   2e-12
GeneScaffold_2821 dna:genescaffold genescaffold:gadMor1:GeneScaf...    72   4e-12
contig300910 dna:contig contig::contig300910:1:2822:1 REF              70   1e-11
GeneScaffold_168 dna:genescaffold genescaffold:gadMor1:GeneScaff...    64   1e-09
GeneScaffold_1983 dna:genescaffold genescaffold:gadMor1:GeneScaf...    57   1e-07
contig70176 dna:contig contig::contig70176:1:567:1 REF                 51   8e-06

>contig411933 dna:contig contig::contig411933:1:1296:1 REF
          Length = 1296

 Score = 95.5 bits (236), Expect = 4e-19
 Identities = 42/53 (79%), Positives = 46/53 (86%)
 Frame = +1

Query: 109 GSQEPGTEAQIKDFAKSYNAEFPMFSKIDVNGDNAHPLWKWLKEQPNGKGFLG 161
           G +EPGTE QIK FA SYNA F MFSK++VNG NAHPLWKWLKEQPNG+GFLG
Sbjct: 691 GLKEPGTETQIKAFAASYNARFDMFSKVNVNGANAHPLWKWLKEQPNGRGFLG 849



 Score = 87.8 bits (216), Expect = 8e-17
 Identities = 42/49 (85%), Positives = 46/49 (93%)
 Frame = +1

Query: 63  RGDVVIIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQ 111
           RG+VVII NVASK GKTPVNYSQFAEMHAKY+E+GLRILAFPSNQF +Q
Sbjct: 250 RGEVVIITNVASK*GKTPVNYSQFAEMHAKYSERGLRILAFPSNQFANQ 396


>GeneScaffold_3230 dna:genescaffold genescaffold:gadMor1:GeneScaffold_3230:1:185225:1 REF
          Length = 185225

 Score = 87.0 bits (214), Expect = 1e-16
 Identities = 37/56 (66%), Positives = 47/56 (83%)
 Frame = +2

Query: 111   QEPGTEAQIKDFAKSYNAEFPMFSKIDVNGDNAHPLWKWLKEQPNGKGFLGNGIKW 166
             QEPGT+A+IK+FAK ++A+F +FSKI+VNG+ AHPLWKW+K QP GKG LG   KW
Sbjct: 45020 QEPGTDAEIKEFAKGFDAKFDLFSKIEVNGEGAHPLWKWMKAQPKGKGLLG---KW 45178



 Score = 72.8 bits (177), Expect = 3e-12
 Identities = 36/49 (73%), Positives = 38/49 (77%)
 Frame = +1

Query: 63    RGDVVIIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQ 111
             RG V II NVASK GKT VNY+QFAEMH KYA+K L IL FP NQFG Q
Sbjct: 44266 RGCVCIITNVASK*GKTRVNYTQFAEMHVKYADKDLHILGFPCNQFGGQ 44412


>GeneScaffold_1266 dna:genescaffold genescaffold:gadMor1:GeneScaffold_1266:1:488504:1 REF
          Length = 488504

 Score = 80.9 bits (198), Expect = 1e-14
 Identities = 62/225 (27%), Positives = 90/225 (40%), Gaps = 65/225 (28%)
 Frame = -1

Query: 41     ATSIYDFSAKDIDGNEVSLEKYRGDVVIIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRI 100
              A + YD  A  ++G  V    +RG VV+I NVAS  G T  +YS+ +++ +KY  + L +
Sbjct: 477200 AKTFYDLKATTLEGATVDFNVFRGRVVLIENVASL*GTTTRDYSELSQLQSKYPHR-LVV 477024

Query: 101    LAFPSNQFGSQEPGT-------EAQIKDF------------------------------- 122
              L FP NQFG Q   T       +  I+ F                               
Sbjct: 477023 LGFPCNQFGYQVSSTAQ*DSTSDMPIEAFDL*PHLPLNALSKITVCVCVCHLQENLSNGE 476844

Query: 123    ----------AKSYNAEFPMFSKIDVNGDNAHPLWKWLKEQ-------PNG--------- 156
                           Y   F +F K +VNG N HP++ +LK++       PN          
Sbjct: 476843 ILNSLKYVRPGAGYEPGFTVFEKCEVNGANPHPVFAYLKDKLPYPEDDPNSLMQDPKHLI 476664

Query: 157    -KGFLGNGIKWNFTKFLINREGQVMKRYGPMDDPSVVEKDLPKYL 200
                    N + WNF KFL+  EG+  KRY        +E D+ + L
Sbjct: 476663 WTPICRNDVSWNFEKFLVGPEGEPFKRYSKKFPTISIEPDIQRLL 476529


>contig311944 dna:contig contig::contig311944:1:640:1 REF
          Length = 640

 Score = 73.2 bits (178), Expect = 2e-12
 Identities = 36/49 (73%), Positives = 39/49 (79%)
 Frame = -3

Query: 63  RGDVVIIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQ 111
           RG V II NVASK GKT VNY+QFAEMH KYA+K L IL FP NQFG+Q
Sbjct: 479 RGCVCIITNVASK*GKTRVNYTQFAEMHVKYADKDLHILGFPCNQFGAQ 333


>GeneScaffold_2821 dna:genescaffold genescaffold:gadMor1:GeneScaffold_2821:1:245547:1
            REF
          Length = 245547

 Score = 72.4 bits (176), Expect = 4e-12
 Identities = 37/93 (39%), Positives = 54/93 (58%), Gaps = 1/93 (1%)
 Frame = -1

Query: 66   VVIIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQEPGTEAQIKDFAKS 125
            V ++VNVAS  G T  +Y    ++   Y      +LAFP NQFG QEPG++ +I  F + 
Sbjct: 7941 VSLVVNVASDCGFTEEHYRDLQQLQKDYGPYHFNVLAFPCNQFGQQEPGSDKEIDSFVRR 7762

Query: 126  -YNAEFPMFSKIDVNGDNAHPLWKWLKEQPNGK 157
             Y   FP+FSKI V G  A+ ++K+L  + + K
Sbjct: 7761 VYGVSFPLFSKIAVVGTGANNVYKYLAGE*SNK 7663


>contig300910 dna:contig contig::contig300910:1:2822:1 REF
          Length = 2822

 Score = 70.5 bits (171), Expect = 1e-11
 Identities = 37/87 (42%), Positives = 49/87 (56%), Gaps = 1/87 (1%)
 Frame = -1

Query: 68   IIVNVASKXGKTPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQEPGTEAQIKDFAKS-Y 126
            ++VNVAS+  +T  NY    E+H +       +LAFP  QFG  EP    +I+ FAKS Y
Sbjct: 1274 LVVNVASQSEQTEENYRSLQELHRELGTSHFNVLAFPCGQFGLTEPWGTREIEAFAKSTY 1095

Query: 127  NAEFPMFSKIDVNGDNAHPLWKWLKEQ 153
               FP FSKI + G  A P +K+L  Q
Sbjct: 1094 GVTFPFFSKIKIMGSEADPAFKFLTGQ 1014


>GeneScaffold_168 dna:genescaffold genescaffold:gadMor1:GeneScaffold_168:1:1723178:1 REF
          Length = 1723178

 Score = 63.9 bits (154), Expect = 1e-09
 Identities = 37/93 (39%), Positives = 53/93 (56%)
 Frame = -2

Query: 19     GLVLVIGLYLFTMSAPTEDWQTATSIYDFSAKDIDGNEVSLEKYRGDVVIIVNVASKXGK 78
              G+ LV+    FTM         A+  ++ +AK + G  ++    +G VV+I NVAS  G 
Sbjct: 730324 GVCLVVTPSTFTM---------ASKFFEITAKRLTGEVLNFSSLKGKVVLIQNVASL*GT 730172

Query: 79     TPVNYSQFAEMHAKYAEKGLRILAFPSNQFGSQ 111
              T  +Y Q  E+HA+Y+ +GL ILA P NQFG Q
Sbjct: 730171 TSRDYLQMNELHARYSGQGLVILAVPCNQFGHQ 730073


>GeneScaffold_1983 dna:genescaffold genescaffold:gadMor1:GeneScaffold_1983:1:630567:1 REF
          Length = 630567

 Score = 57.4 bits (137), Expect = 1e-07
 Identities = 30/44 (68%), Positives = 33/44 (75%)
 Frame = +2

Query: 20     LVLVIGLYLFTMSAPTEDWQTATSIYDFSAKDIDGNEVSLEKYR 63
              LVL+ G   F + A  +D QTATSIYDFSA DIDGN VSLEKYR
Sbjct: 628622 LVLLNG*LCFLLQAAAQDPQTATSIYDFSATDIDGNVVSLEKYR 628753


>contig70176 dna:contig contig::contig70176:1:567:1 REF
          Length = 567

 Score = 51.2 bits (121), Expect = 8e-06
 Identities = 26/50 (52%), Positives = 33/50 (66%)
 Frame = +2

Query: 14  TLASCGLVLVIGLYLFTMSAPTEDWQTATSIYDFSAKDIDGNEVSLEKYR 63
           +L  CG+ L+  L      A   DW++A SIY+F A DIDGNEVSLEKY+
Sbjct: 203 SLTLCGV*LISSL----QCAQAVDWESAKSIYEFRATDIDGNEVSLEKYK 340


  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.318    0.136    0.415 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 173,680,452
Number of Sequences: 398859
Number of extensions: 2096895
Number of successful extensions: 7807
Number of sequences better than 1.0e-04: 9
Number of HSP's better than  0.0 without gapping: 6404
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): 7802
length of query: 200
length of database: 277,371,529
effective HSP length: 115
effective length of query: 85
effective length of database: 231,502,744
effective search space: 19677733240
effective search space used: 19677733240
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.7 bits)
S2: 112 (47.8 bits)
