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|159155224|gb|AAI54785.1| Thioredoxin reductase 1 [Danio
rerio]
         (602 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_2145 dna:genescaffold genescaffold:gadMor1:GeneScaf...   142   1e-36
GeneScaffold_2535 dna:genescaffold genescaffold:gadMor1:GeneScaf...    91   5e-17
GeneScaffold_2281 dna:genescaffold genescaffold:gadMor1:GeneScaf...    73   1e-11
contig443014 dna:contig contig::contig443014:1:589:1 REF               64   9e-09
scaffold01964 dna:scaffold scaffold:gadMor1:scaffold01964:1:8461...    58   5e-07
contig689731 dna:contig contig::contig689731:1:960:1 REF               55   2e-06
contig350844 dna:contig contig::contig350844:1:946:1 REF               55   2e-06
contig407116 dna:contig contig::contig407116:1:860:1 REF               54   5e-06
GeneScaffold_3769 dna:genescaffold genescaffold:gadMor1:GeneScaf...    52   3e-05

>GeneScaffold_2145 dna:genescaffold genescaffold:gadMor1:GeneScaffold_2145:1:113983:1 REF
          Length = 113983

 Score =  142 bits (359), Expect = 1e-32
 Identities = 66/77 (85%), Positives = 74/77 (96%)
 Frame = -3

Query: 283   SDDLFSLPYCPGKTLVVGASYVALECGGFLAGLGLDVTIMVRSILLRGFDQDMADRAGEY 342
             SDDLFSLPYCPGKTL++GASYVALEC GFLAGLGLDVT+MVRSILLRGFDQ+MA+RAGE+
Sbjct: 39926 SDDLFSLPYCPGKTLLIGASYVALECAGFLAGLGLDVTVMVRSILLRGFDQEMANRAGEH 39747

Query: 343   METHGVKFLRKFVPTKI 359
             ME HGVKFLRK+VP K+
Sbjct: 39746 MEEHGVKFLRKYVPVKV 39696



 Score =  112 bits (280), Expect(2) = 1e-36
 Identities = 58/93 (62%), Positives = 62/93 (66%), Gaps = 29/93 (31%)
 Frame = -2

Query: 132   KEAATLGKKVMVLDYVVPTPQGTAW-----------------------------GLGGTC 162
             +EAATLGKKVMVLDYVVPTP+GT W                             GLGGTC
Sbjct: 41289 QEAATLGKKVMVLDYVVPTPKGTTWGENASPCATLCVPPCVACSSF*SLCIWGTGLGGTC 41110

Query: 163   VNVGCIPKKLMHQTALLGTAMEDARKFGWEFAE 195
             VNVGCIPKKLMHQTALLGTAM+DARK+GWEF E
Sbjct: 41109 VNVGCIPKKLMHQTALLGTAMQDARKYGWEFDE 41011



 Score =  103 bits (257), Expect = 8e-21
 Identities = 52/64 (81%), Positives = 58/64 (90%), Gaps = 2/64 (3%)
 Frame = -2

Query: 406   DKAGVKIN--EKNGKVPVNDEEQTNVPHIYAIGDILEGKWELTPVAIQAGKLLARRLYAG 463
             DK G K +   +NGKVPV+DEEQTNVP+IYAIGDILEGKWELTPVAIQAGKLLARRLYAG
Sbjct: 36918 DKQGNKRSPLSRNGKVPVDDEEQTNVPYIYAIGDILEGKWELTPVAIQAGKLLARRLYAG 36739

Query: 464   ATMK 467
             A++K
Sbjct: 36738 ASLK 36727



 Score = 95.1 bits (235), Expect = 3e-18
 Identities = 42/48 (87%), Positives = 44/48 (91%)
 Frame = -1

Query: 197   VTHNWETMKTAVNNYIGSLNWGYRVSLRDKNVNYVNAYAEFVEPHKIK 244
             V HNW  MK A+NNYIGSLNWGYRV+LRDKNVNYVNAYAEFVEPHKIK
Sbjct: 40774 VKHNWVIMKEAINNYIGSLNWGYRVALRDKNVNYVNAYAEFVEPHKIK 40631



 Score = 92.4 bits (228), Expect = 2e-17
 Identities = 39/45 (86%), Positives = 44/45 (97%)
 Frame = -3

Query: 537   RVIGFHYLGPNAGEVTQGFGAAMKCGITKDQLDNTIGIHPTCAEI 581
             RVIGFHYLGPNAGEVTQG+G A+KCG+TK+QLDNTIGIHPTCAE+
Sbjct: 33770 RVIGFHYLGPNAGEVTQGYGVAIKCGVTKEQLDNTIGIHPTCAEV 33636



 Score = 74.3 bits (181), Expect = 5e-12
 Identities = 34/38 (89%), Positives = 36/38 (94%)
 Frame = -1

Query: 245   ATNKRGKETFYTAAQFVLATGERPRYLGIPGDKEFCIT 282
             AT KRGKET +TAA+FVLATGERPRYLGIPGDKEFCIT
Sbjct: 40447 ATTKRGKETMHTAARFVLATGERPRYLGIPGDKEFCIT 40334



 Score = 73.2 bits (178), Expect = 1e-11
 Identities = 31/38 (81%), Positives = 35/38 (92%)
 Frame = -2

Query: 467   KCDYVNVPTTVFTPMEYGSCGHPEEKAIQMYGQENLEV 504
             +CDYVNVPTTVFTP+EYG CG  EEKAI+MYG+ENLEV
Sbjct: 36447 QCDYVNVPTTVFTPLEYGCCGRSEEKAIEMYGEENLEV 36334



 Score = 69.7 bits (169), Expect = 1e-10
 Identities = 29/36 (80%), Positives = 33/36 (91%)
 Frame = -2

Query: 501   NLEVYHSLFWPLEFTVPGRDNNKCYAKIICNKLDNL 536
             +L +YHSLFWPLEFTV  RDNNKCY+KIICNKLDN+
Sbjct: 35661 HL*IYHSLFWPLEFTVANRDNNKCYSKIICNKLDNV 35554



 Score = 65.1 bits (157), Expect(2) = 1e-36
 Identities = 30/47 (63%), Positives = 35/47 (74%)
 Frame = -1

Query: 91    CDNTMKAHKDGVLQKLLGEGSEVYDYDLIVIGGGSGGLACSKEAATL 137
             C    +AH DG+LQ+LL   +E YDYDLIVIGGGSGGLACSK  A +
Sbjct: 41506 CPPPRQAHTDGLLQQLLNPDAETYDYDLIVIGGGSGGLACSKVCANM 41366



 Score = 62.8 bits (151), Expect = 2e-08
 Identities = 26/40 (65%), Positives = 36/40 (90%)
 Frame = -1

Query: 1     MPPIENDAGREQIRSKIKELIDSSAVVVFSKSFCPFCVKV 40
             MPP+E DAG+  ++++I+ELI+S+ VVVFSKS+CPFCVKV
Sbjct: 42652 MPPVEKDAGKTGLKTRIQELIESNQVVVFSKSYCPFCVKV 42533



 Score = 62.0 bits (149), Expect = 3e-08
 Identities = 27/35 (77%), Positives = 31/35 (88%)
 Frame = -1

Query: 60    EDGTNYQDLLHEMTGQKTVPNVFINKKHIGGCDNT 94
             ++GTNYQ+LL EMTGQKTVPNVFINK H+GG D T
Sbjct: 41848 DEGTNYQELLVEMTGQKTVPNVFINKTHVGGWDKT 41744



 Score = 57.4 bits (137), Expect = 6e-07
 Identities = 25/32 (78%), Positives = 30/32 (93%)
 Frame = -2

Query: 359   IEQLEAGTPGRIKVTAKSTESEEFFEGEYNTV 390
             +E+LEAGTPGR++VTAKSTES+E  EGEYNTV
Sbjct: 39060 VEELEAGTPGRLRVTAKSTESDELLEGEYNTV 38965


>GeneScaffold_2535 dna:genescaffold genescaffold:gadMor1:GeneScaffold_2535:1:381514:1 REF
          Length = 381514

 Score = 90.9 bits (224), Expect = 5e-17
 Identities = 37/49 (75%), Positives = 44/49 (89%)
 Frame = +1

Query: 153   GTAWGLGGTCVNVGCIPKKLMHQTALLGTAMEDARKFGWEFAEQVTHNW 201
             GT WGLGGTCVNVGCIPKKLMHQ ALLGTA++DARK+GW+  + V+H+W
Sbjct: 84484 GTKWGLGGTCVNVGCIPKKLMHQAALLGTAVKDARKYGWKIKDPVSHDW 84630



 Score = 52.4 bits (124), Expect = 2e-05
 Identities = 23/43 (53%), Positives = 28/43 (65%)
 Frame = +3

Query: 560    KCGITKDQLDNTIGIHPTCAEIFTTMEVTKSSGGDITQSGCXG 602
              +CG T   L  T+GIHPTCAE  T + +TK SG D T +GC G
Sbjct: 103470 RCGATLSHLKETVGIHPTCAEELTKVSITKRSGLDATVTGC*G 103598



 Score = 50.8 bits (120), Expect = 6e-05
 Identities = 24/34 (70%), Positives = 26/34 (76%)
 Frame = +1

Query: 300   GASYVALECGGFLAGLGLDVTIMVRSILLRGFDQ 333
             G S VALEC GFL G+G D  +MVRSI LRGFDQ
Sbjct: 90298 GDSDVALECAGFLTGIGQDTKVMVRSIALRGFDQ 90399



 Score = 50.4 bits (119), Expect = 8e-05
 Identities = 33/98 (33%), Positives = 44/98 (44%), Gaps = 37/98 (37%)
 Frame = +2

Query: 440   EGKWELTPVAIQAGKLLARRLYAGATMKCDYVNVPT------------------------ 475
             +G+ ELTP A++AG+LLARRL   +T   DY NV T                        
Sbjct: 97448 QGRPELTPTAVKAGRLLARRLAGQSTQLMDYDNVSTMRTQVYHTIHFTHTKSLILFWVFL 97627

Query: 476   -------------TVFTPMEYGSCGHPEEKAIQMYGQE 500
                          TVFTP+EY   G  EE+A+  +GQ+
Sbjct: 97628 SCNFSAYVYQVATTVFTPLEYSCVGLSEEEAVGRHGQD 97741


>GeneScaffold_2281 dna:genescaffold genescaffold:gadMor1:GeneScaffold_2281:1:51721:1 REF
          Length = 51721

 Score = 73.2 bits (178), Expect = 1e-11
 Identities = 39/82 (47%), Positives = 55/82 (67%), Gaps = 1/82 (1%)
 Frame = -1

Query: 281   ITSDDLFSLPYCPGKTLVVGASYVALECGGFLAGLGLDVTIMV-RSILLRGFDQDMADRA 339
             ++S++LFS P  PG+TLVVG S   +EC GFLAGLG +VT+M+  S LL GFD  M+ + 
Sbjct: 23959 LSSNELFSRPCPPGRTLVVGGSPEGVECAGFLAGLGHEVTLMLPTSDLLPGFDLKMSQKI 23780

Query: 340   GEYMETHGVKFLRKFVPTKIEQ 361
               +M  HGV FL +   +K+ +
Sbjct: 23779 ENHMLVHGVHFLYEHTLSKVSK 23714


>contig443014 dna:contig contig::contig443014:1:589:1 REF
          Length = 589

 Score = 63.5 bits (153), Expect = 9e-09
 Identities = 28/54 (51%), Positives = 37/54 (68%)
 Frame = +2

Query: 537 RVIGFHYLGPNAGEVTQGFGAAMKCGITKDQLDNTIGIHPTCAEIFTTMEVTKS 590
           RV+G H +GP+AG++ QGF AAM CG+TK QLD T G  PT A++     +T S
Sbjct: 317 RVVGLHLMGPHAGDIVQGFSAAMTCGLTKGQLDATAGQQPTAAQVGGRTRLTLS 478


>scaffold01964 dna:scaffold scaffold:gadMor1:scaffold01964:1:84619:1 REF
          Length = 84619

 Score = 57.8 bits (138), Expect = 5e-07
 Identities = 24/61 (39%), Positives = 39/61 (63%)
 Frame = -2

Query: 18    KELIDSSAVVVFSKSFCPFCVKVKDLFKELNVKCNTIELDLMEDGTNYQDLLHEMTGQKT 77
             +E++  + +V+FSK+ CP+C   K++F E+      IELD  +DG   Q+ L +MTG +T
Sbjct: 60975 QEVVSHNCIVIFSKTTCPYCRMAKNVFNEIGANYKVIELDEHDDGRRVQEALAQMTGART 60796

Query: 78    V 78
             V
Sbjct: 60795 V 60793


>contig689731 dna:contig contig::contig689731:1:960:1 REF
          Length = 960

 Score = 55.5 bits (132), Expect = 2e-06
 Identities = 21/46 (45%), Positives = 32/46 (69%)
 Frame = -1

Query: 460 LYAGATMKCDYVNVPTTVFTPMEYGSCGHPEEKAIQMYGQENLEVY 505
           +Y    ++CDY N+P  +FTP+EY +CG  EE+A   +G+ N+EVY
Sbjct: 249 IYGVLCLQCDYSNIPLVLFTPLEYSACGLSEERANATFGEANIEVY 112


>contig350844 dna:contig contig::contig350844:1:946:1 REF
          Length = 946

 Score = 55.5 bits (132), Expect = 2e-06
 Identities = 21/46 (45%), Positives = 32/46 (69%)
 Frame = -2

Query: 460 LYAGATMKCDYVNVPTTVFTPMEYGSCGHPEEKAIQMYGQENLEVY 505
           +Y    ++CDY N+P  +FTP+EY +CG  EE+A   +G+ N+EVY
Sbjct: 246 IYGVLCLQCDYSNIPLVLFTPLEYSACGLSEERANATFGEANIEVY 109


>contig407116 dna:contig contig::contig407116:1:860:1 REF
          Length = 860

 Score = 54.3 bits (129), Expect = 5e-06
 Identities = 21/37 (56%), Positives = 28/37 (75%)
 Frame = +2

Query: 502 LEVYHSLFWPLEFTVPGRDNNKCYAKIICNKLDNLRV 538
           ++VYHS +WPLE+TVP RD N CY K+IC+  D + V
Sbjct: 119 VQVYHSYYWPLEWTVPPRDKNSCYVKVICHVPDQVPV 229


>GeneScaffold_3769 dna:genescaffold genescaffold:gadMor1:GeneScaffold_3769:1:445740:1 REF
          Length = 445740

 Score = 51.6 bits (122), Expect = 3e-05
 Identities = 24/49 (48%), Positives = 29/49 (59%)
 Frame = +2

Query: 537   RVIGFHYLGPNAGEVTQGFGAAMKCGITKDQLDNTIGIHPTCAEIFTTM 585
             +V+G H  G    E+ QGF  A+K G TK   D TI IHPT +E F TM
Sbjct: 75953 QVVGLHMQGLGCDEMLQGFAVAVKMGATKADFDRTIAIHPTSSEEFVTM 76099


  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.137    0.417 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 620,251,958
Number of Sequences: 398859
Number of extensions: 9233806
Number of successful extensions: 34126
Number of sequences better than 1.0e-04: 9
Number of HSP's better than  0.0 without gapping: 16426
Number of HSP's successfully gapped in prelim test: 514
Number of HSP's that attempted gapping in prelim test: 4352
Number of HSP's gapped (non-prelim): 31499
length of query: 602
length of database: 277,371,529
effective HSP length: 128
effective length of query: 474
effective length of database: 226,317,577
effective search space: 107274531498
effective search space used: 107274531498
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: 118 (50.1 bits)
