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|148277097|ref|NP_001007284.2| type I iodothyronine
deiodinase [Danio rerio]
         (254 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_1484 dna:genescaffold genescaffold:gadMor1:GeneScaf...   170   1e-41
contig121407 dna:contig contig::contig121407:1:1277:1 REF             127   1e-28
GeneScaffold_4497 dna:genescaffold genescaffold:gadMor1:GeneScaf...   125   4e-28
GeneScaffold_2685 dna:genescaffold genescaffold:gadMor1:GeneScaf...    49   7e-05

>GeneScaffold_1484 dna:genescaffold genescaffold:gadMor1:GeneScaffold_1484:1:590555:1 REF
          Length = 590555

 Score =  170 bits (431), Expect = 1e-41
 Identities = 88/184 (47%), Positives = 120/184 (65%), Gaps = 26/184 (14%)
 Frame = -3

Query: 9      LRKLFVYISAVLMVCAAILQMSMLKLLSFISPGRMRKIHMKMGERSTMTQNPKFRYEDWG 68
              L+KL++Y+S   +    IL + ++++L  + P   +K+ +KMGE+ TMTQNPKFR+EDWG
Sbjct: 486297 LQKLWIYLSIGCLFFGQILTVFVIRVLYLVWPALAKKLVLKMGEQCTMTQNPKFRFEDWG 486118

Query: 69     PAFFSLAFIKTLFFVNWCSLGLEAFEGHAAPDSALITLDRQKTSVHRFL----------- 117
              P+F SL  +K +F   W SLG EAF G  APDS ++TL+   +++ + L           
Sbjct: 486117 PSFTSLTSLKIIFNHLWLSLGQEAFVGWDAPDSQVVTLEGNTSTILKHLNG**TDALIH* 485938

Query: 118    ---------------KGNRPLVLSFGSCTXPPFLYKLDEFKQLVKDFSNVADFLIVYLAE 162
                              GNRPLVLSFGSCT PPFL+KLDEFK+LVKD+ +VADFL+VY+AE
Sbjct: 485937 WVSSPIAFLM*NFISSGNRPLVLSFGSCT*PPFLFKLDEFKELVKDYRDVADFLVVYIAE 485758

Query: 163    AHAT 166
              AHAT
Sbjct: 485757 AHAT 485746



 Score = 77.8 bits (190), Expect = 1e-13
 Identities = 35/75 (46%), Positives = 49/75 (65%)
 Frame = -3

Query: 159    YLAEAHATDAWAFKNNVDISVHKNXXXXXXXXXXXXKEDPPCPVVVDEMNNITASKYGAL 218
              Y     ++D WAF NNV+I  H+               +P CPVVVD+M+++TA KYG+L
Sbjct: 485391 YWIHIPSSDGWAFVNNVEIKRHQTLGDRLVAAQTLVDLEPGCPVVVDDMSDVTACKYGSL 485212

Query: 219    PERLYVIQSGKVIYQ 233
              PERLYV+Q+GKV+Y+
Sbjct: 485211 PERLYVLQAGKVVYK 485167


>contig121407 dna:contig contig::contig121407:1:1277:1 REF
          Length = 1277

 Score =  127 bits (320), Expect = 1e-28
 Identities = 73/229 (31%), Positives = 126/229 (55%), Gaps = 2/229 (0%)
 Frame = +1

Query: 26  ILQMSMLKLLSFISPGRMRKIHMKMGERSTMTQNPKFRYEDWGPAFFSLAFIKTLFFVNW 85
           +L   ML +L F+     +K+ ++MG+      +P     D     F+L  ++ +++   
Sbjct: 106 LLAAVMLWILDFLCI--RKKVLLRMGQGEGGPDDPPVCVSD-SNKMFTLESLRAVWYGQK 276

Query: 86  CSLGLEAFEGHAAPDSALITL-DRQKTSVHRFLKGNRPLVLSFGSCTXPPFLYKLDEFKQ 144
                 A  GH+AP++ ++ + DRQ   +   ++G RPL+L+FGSC+ PPF+ +L  F++
Sbjct: 277 LDFFKAAHVGHSAPNTEVVAVQDRQTLRILDCMRGKRPLILNFGSCS*PPFMTRLAAFRR 456

Query: 145 LVKDFSNVADFLIVYLAEAHATDAWAFKN-NVDISVHKNXXXXXXXXXXXXKEDPPCPVV 203
           +   ++++AD L+VY+ EAH +D W   +    I  H+              E P   VV
Sbjct: 457 VASQYADIADSLLVYIEEAHPSDGWVSSDAPYQIPKHRCIEDRLRAAQLMLSEVPGSQVV 636

Query: 204 VDEMNNITASKYGALPERLYVIQSGKVIYQGGIGPWGYKPEEVRKVLEK 252
           +D M+N + + YGA  ERLY+++  +V+YQGG GP GY+  E+R  LE+
Sbjct: 637 IDTMDNSSNAAYGAYFERLYILRDQRVVYQGGRGPEGYRISELRDWLEQ 783


>GeneScaffold_4497 dna:genescaffold genescaffold:gadMor1:GeneScaffold_4497:1:397058:1 REF
          Length = 397058

 Score =  125 bits (315), Expect = 4e-28
 Identities = 64/155 (41%), Positives = 95/155 (61%), Gaps = 2/155 (1%)
 Frame = +1

Query: 95    GHAAPDSALITL-DRQKTSVHRFLKGNRPLVLSFGSCTXPPFLYKLDEFKQLVKDFSNVA 153
             G  AP++ ++ L D++++ V  + +G+RPL+L+FGSCT PPF+ +L  F  LV+   +VA
Sbjct: 18853 GCGAPNTEVVQLRDQKRSRVLDYAQGDRPLILNFGSCT*PPFMARLKAFHGLVQQHEDVA 19032

Query: 154   DFLIVYLAEAHATDAWAFKN-NVDISVHKNXXXXXXXXXXXXKEDPPCPVVVDEMNNITA 212
             D+L+VY+ EAH +D W   +    I  H+              E P C VVVD M N T 
Sbjct: 19033 DYLVVYIEEAHPSDGWVSSDAPYQIPKHRCLEDRLSAAQLMRLEVPGCAVVVDNMENATN 19212

Query: 213   SKYGALPERLYVIQSGKVIYQGGIGPWGYKPEEVR 247
             + YGA  +RLY++  GK++YQGG GP GY+  E+R
Sbjct: 19213 AAYGAYFDRLYILHDGKIVYQGGRGPEGYRISELR 19317


>GeneScaffold_2685 dna:genescaffold genescaffold:gadMor1:GeneScaffold_2685:1:119182:1
            REF
          Length = 119182

 Score = 48.9 bits (115), Expect = 7e-05
 Identities = 22/56 (39%), Positives = 35/56 (62%)
 Frame = -3

Query: 198  PPCPVVVDEMNNITASKYGALPERLYVIQSGKVIYQGGIGPWGYKPEEVRKVLEKS 253
            P C +V D M+N   + YG   ER+ ++Q  K+ Y GG GP+ Y  ++VR+ LE++
Sbjct: 4418 PQCRLVADSMDNNANAAYGVANERVCIVQRRKIAYLGGKGPFFYNLKDVRQWLEQT 4251


  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.324    0.138    0.423 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 251,204,350
Number of Sequences: 398859
Number of extensions: 3305438
Number of successful extensions: 12077
Number of sequences better than 1.0e-04: 4
Number of HSP's better than  0.0 without gapping: 9404
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): 12071
length of query: 254
length of database: 277,371,529
effective HSP length: 119
effective length of query: 135
effective length of database: 229,907,308
effective search space: 31037486580
effective search space used: 31037486580
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: 41 (22.0 bits)
S2: 114 (48.5 bits)
