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= SPP00000108_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Mus musculus # Complete
         (452 letters)

Database: /cursos/BI/genomes/project_2013/Mustela_putorius_furo/genome
.fa 
           4245 sequences; 2,400,181,698 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gb|JH418158.1| Mustela putorius furo unplaced genomic scaffold s...   518   e-145
gb|JH418181.1| Mustela putorius furo unplaced genomic scaffold s...   111   7e-23
gb|JH417973.1| Mustela putorius furo unplaced genomic scaffold s...    37   2.1  

>gb|JH418158.1| Mustela putorius furo unplaced genomic scaffold scaffold00186, whole genome
               shotgun sequence
          Length = 3810862

 Score =  518 bits (1334), Expect = e-145
 Identities = 256/322 (79%), Positives = 273/322 (84%)
 Frame = -1

Query: 111     PGSAFPSLSIGMDSCVIPLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGI 170
               P +  P+L IGMDSCV+PLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGI
Sbjct: 2928982 PSAPSPALGIGMDSCVVPLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGI 2928803

Query: 171     TECDNXXXXXXXXXXXXEKEREKVTPLMIKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVA 230
               TECDN            ++EREKVTPL+I+GFRD                NPW+I+GGVA
Sbjct: 2928802 TECDNMLMLLSVSQSMTDEEREKVTPLLIRGFRDAAEEGGTAVTGGQTVVNPWVIVGGVA 2928623

Query: 231     TVVCQQNEFIMPDSAVVGDVLVLTKPLGTQVAANAHQWLDNPEKWNKIKMVVSREEVELA 290
               TVVCQ +E IMPDSAVVGDVLVLTKPLGTQVA NAHQWLD+PE+W+K+ MVVSREEVELA
Sbjct: 2928622 TVVCQPSELIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDSPERWSKVSMVVSREEVELA 2928443

Query: 291     YQEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNLAKQQKNEVSFVIHNLPXX 350
               YQEAMF+MATLNRTAA LMHTFNAHAATD+TGFGILGHSQNLAKQQ+NEVSFVIHNLP  
Sbjct: 2928442 YQEAMFSMATLNRTAAALMHTFNAHAATDVTGFGILGHSQNLAKQQRNEVSFVIHNLPII 2928263

Query: 351     XXXXXXXXXXGRFGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQAWIVGIV 410
                         GRFGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQAWIVGIV
Sbjct: 2928262 AKMAAISKASGRFGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQAWIVGIV 2928083

Query: 411     EKGNRTARIIDKPRVIEVLPRG 432
               EKGNRTARIIDKPRVIEVLPRG
Sbjct: 2928082 EKGNRTARIIDKPRVIEVLPRG 2928017


>gb|JH418181.1| Mustela putorius furo unplaced genomic scaffold scaffold00209, whole genome
               shotgun sequence
          Length = 3329659

 Score =  111 bits (278), Expect = 7e-23
 Identities = 51/57 (89%), Positives = 55/57 (96%)
 Frame = -2

Query: 375     GGLLICLPREQAARFCSEIKSSKYGEGHQAWIVGIVEKGNRTARIIDKPRVIEVLPR 431
               GGLLICLPREQAARFC+EIKS KYGEGHQAWI+GIVEKGNRTARIIDKPR+IEV P+
Sbjct: 1462692 GGLLICLPREQAARFCAEIKSPKYGEGHQAWIIGIVEKGNRTARIIDKPRIIEVAPQ 1462522



 Score =  101 bits (252), Expect = 7e-20
 Identities = 52/72 (72%), Positives = 55/72 (76%), Gaps = 1/72 (1%)
 Frame = -3

Query: 305     AAGLMHTFNAHAATDITGFGILGHSQNLAKQQKNEVSFVIHNLPXXXXXXXXXXXXGR-F 363
               AAGLMHTFNAHAATDITGFGILGH+QNLAKQQ+NEVSFVIHNLP            G  F
Sbjct: 1464512 AAGLMHTFNAHAATDITGFGILGHAQNLAKQQRNEVSFVIHNLPVLAKMAAVSKACGNMF 1464333

Query: 364     GLLQGTSAETSG 375
               GL+ GT  ETSG
Sbjct: 1464332 GLMHGTCPETSG 1464297



 Score = 69.3 bits (168), Expect = 4e-10
 Identities = 30/41 (73%), Positives = 35/41 (85%)
 Frame = -1

Query: 113     SAFPSLSIGMDSCVIPLRHGGLSLVQTTDFFYPLVEDPYMM 153
               + F    IGMD+CVIPLRHGGLSLVQTTD+ YP+V+DPYMM
Sbjct: 1482106 TTFIFAGIGMDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMM 1481984



 Score = 62.8 bits (151), Expect = 4e-08
 Identities = 28/37 (75%), Positives = 33/37 (89%)
 Frame = -1

Query: 272     PEKWNKIKMVVSREEVELAYQEAMFNMATLNRTAAGL 308
               PEKWNKIK+VV++E+VELAYQEAM NMA LNRT  G+
Sbjct: 1472758 PEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTGTGV 1472648



 Score = 59.3 bits (142), Expect = 4e-07
 Identities = 27/52 (51%), Positives = 33/52 (63%)
 Frame = -2

Query: 190     EREKVTPLMIKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVATVVCQQNEFIM 241
               ER+KV PL+I+GF+D                NPWI++GGVAT VCQ NEFIM
Sbjct: 1478394 ERDKVMPLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNEFIM 1478239



 Score = 52.8 bits (125), Expect = 4e-05
 Identities = 24/29 (82%), Positives = 25/29 (86%)
 Frame = -2

Query: 242     PDSAVVGDVLVLTKPLGTQVAANAHQWLD 270
               PD+AV GDVLVLTKPLGTQVA   HQWLD
Sbjct: 1474794 PDNAVPGDVLVLTKPLGTQVAVAVHQWLD 1474708



 Score = 49.3 bits (116), Expect = 4e-04
 Identities = 23/47 (48%), Positives = 30/47 (63%)
 Frame = -2

Query: 33      RSFSNYRPFEPQTLGFSPSWRLTSFSGMKGUGCKVPQETLLKLLEGL 79
               R+ S    F P++     S+RLT F+ +KG GCKVPQ+ L KLLE L
Sbjct: 1487520 RTMSTRESFNPESYELDKSFRLTRFTELKGTGCKVPQDVLQKLLESL 1487380



 Score = 49.3 bits (116), Expect = 4e-04
 Identities = 21/22 (95%), Positives = 22/22 (100%)
 Frame = -2

Query: 154     GRIACANVLSDLYAMGITECDN 175
               GRIACANVLSDLYAMG+TECDN
Sbjct: 1479381 GRIACANVLSDLYAMGVTECDN 1479316


>gb|JH417973.1| Mustela putorius furo unplaced genomic scaffold scaffold00001, whole genome
                shotgun sequence
          Length = 38702089

 Score = 37.0 bits (84), Expect = 2.1
 Identities = 27/63 (42%), Positives = 33/63 (52%), Gaps = 4/63 (6%)
 Frame = +3

Query: 29       WSAGRSFSN--YRPFEPQTLGFSPSWRLTSFSGMKGUGCKVPQETLLKL--LEGLTRPAL 84
                W+AG   S   +R F   +LGFS SW L+  S  +  G KVP   LLK+   EGL  P  
Sbjct: 32370258 WAAGSRTSEKGWRHFRITSLGFSSSWGLSLASHPRPLGPKVP-HYLLKIGPREGLATPLS 32370434

Query: 85       QPP 87
                 PP
Sbjct: 32370435 LPP 32370443


  Database: /cursos/BI/genomes/project_2013/Mustela_putorius_furo/geno
  me.fa
    Posted date:  Oct 24, 2012  8:21 PM
  Number of letters in database: 2,400,181,698
  Number of sequences in database:  4245
  
Lambda     K      H
   0.319    0.136    0.413 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 1,046,682,289
Number of Sequences: 4245
Number of extensions: 14561052
Number of successful extensions: 60495
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 101
Number of HSP's successfully gapped in prelim test: 77
Number of HSP's that attempted gapping in prelim test: 54053
Number of HSP's gapped (non-prelim): 19985
length of query: 452
length of database: 800,060,566
effective HSP length: 134
effective length of query: 318
effective length of database: 799,491,736
effective search space: 254238372048
effective search space used: 254238372048
frameshift window, decay const: 40,  0.1
T: 13
A: 40
X1: 16 ( 7.4 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.7 bits)
S2: 78 (34.7 bits)
