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= SPP00000107_1.0 # Protein # Selenophosphate synthetase 1 (SPS1)
# Mus musculus # Complete
         (392 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...   542   e-152
gb|JH418181.1| Mustela putorius furo unplaced genomic scaffold s...   150   1e-34

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

 Score =  542 bits (1397), Expect = e-152
 Identities = 259/323 (80%), Positives = 294/323 (91%)
 Frame = -1

Query: 59      AVMPRLGIGMDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMMGRIACANVLSDLYAMGVTE 118
               A  P LGIGMD+CV+PLRHGGLSLVQTTD+ YP+V+DPYMMGRIACANVLSDLYAMG+TE
Sbjct: 2928976 APSPALGIGMDSCVVPLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGITE 2928797

Query: 119     CDNMLMLLGVSNKMTDRERDKVIPLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATT 178
               CDNMLMLL VS  MTD ER+KV PL+I+GF+DAAEE GT+VTGGQTV+NPW+++GGVAT 
Sbjct: 2928796 CDNMLMLLSVSQSMTDEEREKVTPLLIRGFRDAAEEGGTAVTGGQTVVNPWVIVGGVATV 2928617

Query: 179     VCQPNEFIMPDNAVPGDVLVLTKPLGTQVAVAVHQWLDIPEKWNKIKLVVTQEDVELAYQ 238
               VCQP+E IMPD+AV GDVLVLTKPLGTQVAV  HQWLD PE+W+K+ +VV++E+VELAYQ
Sbjct: 2928616 VCQPSELIMPDSAVVGDVLVLTKPLGTQVAVNAHQWLDSPERWSKVSMVVSREEVELAYQ 2928437

Query: 239     EAMMNMARLNRTAAGLMHTFNAHAATDITGFGILGHAQNLAKQQRNEVSFVIHNLPVLAK 298
               EAM +MA LNRTAA LMHTFNAHAATD+TGFGILGH+QNLAKQQRNEVSFVIHNLP++AK
Sbjct: 2928436 EAMFSMATLNRTAAALMHTFNAHAATDVTGFGILGHSQNLAKQQRNEVSFVIHNLPIIAK 2928257

Query: 299     MAAVSKACGNMFGLMHGTCPETSGGLLICLPREQAARFCAEIKSPKYGEGHQAWIIGIVE 358
               MAA+SKA G  FGL+ GT  ETSGGLLICLPREQAARFC+EIKS KYGEGHQAWI+GIVE
Sbjct: 2928256 MAAISKASGR-FGLLQGTSAETSGGLLICLPREQAARFCSEIKSSKYGEGHQAWIVGIVE 2928080

Query: 359     KGNRTARIIDKPRIIEVAPQVAT 381
               KGNRTARIIDKPR+IEV P+ AT
Sbjct: 2928079 KGNRTARIIDKPRVIEVLPRGAT 2928011


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

 Score =  150 bits (378), Expect = 1e-34
 Identities = 72/72 (100%), Positives = 72/72 (100%)
 Frame = -3

Query: 251     AAGLMHTFNAHAATDITGFGILGHAQNLAKQQRNEVSFVIHNLPVLAKMAAVSKACGNMF 310
               AAGLMHTFNAHAATDITGFGILGHAQNLAKQQRNEVSFVIHNLPVLAKMAAVSKACGNMF
Sbjct: 1464512 AAGLMHTFNAHAATDITGFGILGHAQNLAKQQRNEVSFVIHNLPVLAKMAAVSKACGNMF 1464333

Query: 311     GLMHGTCPETSG 322
               GLMHGTCPETSG
Sbjct: 1464332 GLMHGTCPETSG 1464297



 Score =  147 bits (372), Expect = 7e-34
 Identities = 71/71 (100%), Positives = 71/71 (100%)
 Frame = -2

Query: 322     GGLLICLPREQAARFCAEIKSPKYGEGHQAWIIGIVEKGNRTARIIDKPRIIEVAPQVAT 381
               GGLLICLPREQAARFCAEIKSPKYGEGHQAWIIGIVEKGNRTARIIDKPRIIEVAPQVAT
Sbjct: 1462692 GGLLICLPREQAARFCAEIKSPKYGEGHQAWIIGIVEKGNRTARIIDKPRIIEVAPQVAT 1462513

Query: 382     QNVNPTPGATS 392
               QNVNPTPGATS
Sbjct: 1462512 QNVNPTPGATS 1462480



 Score =  135 bits (340), Expect = 4e-30
 Identities = 73/108 (67%), Positives = 84/108 (77%), Gaps = 5/108 (4%)
 Frame = -2

Query: 1       MSTRESFNPETYELDKSFRLTRFTELKGTGCKVPQDVLQKLLESLQENHFQEDEQFLGAV 60
               MSTRESFNPE+YELDKSFRLTRFTELKGTGCKVPQDVLQKLLESLQENHFQEDEQFLGAV
Sbjct: 1487514 MSTRESFNPESYELDKSFRLTRFTELKGTGCKVPQDVLQKLLESLQENHFQEDEQFLGAV 1487335

Query: 61      MPRLGI-GMDTCVIPLRHGGL----SLVQTTDYIYPIVDDPYMMGRIA 103
               MPRLG+  + +  +P R        S+ Q  ++ +P+V  PY    I+
Sbjct: 1487334 MPRLGM*TIISFTLPKRRASWLQL*SMSQKPNFCHPVV--PYSTSLIS 1487197



 Score =  107 bits (266), Expect = 1e-21
 Identities = 51/52 (98%), Positives = 52/52 (100%)
 Frame = -2

Query: 136     ERDKVIPLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNEFIM 187
               ERDKV+PLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNEFIM
Sbjct: 1478394 ERDKVMPLIIQGFKDAAEEAGTSVTGGQTVLNPWIVLGGVATTVCQPNEFIM 1478239



 Score = 78.2 bits (191), Expect = 7e-13
 Identities = 35/35 (100%), Positives = 35/35 (100%)
 Frame = -1

Query: 65      GIGMDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMM 99
               GIGMDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMM
Sbjct: 1482088 GIGMDTCVIPLRHGGLSLVQTTDYIYPIVDDPYMM 1481984



 Score = 75.5 bits (184), Expect = 4e-12
 Identities = 35/36 (97%), Positives = 36/36 (100%)
 Frame = -2

Query: 100     GRIACANVLSDLYAMGVTECDNMLMLLGVSNKMTDR 135
               GRIACANVLSDLYAMGVTECDNMLMLLG+SNKMTDR
Sbjct: 1479381 GRIACANVLSDLYAMGVTECDNMLMLLGISNKMTDR 1479274



 Score = 72.4 bits (176), Expect = 4e-11
 Identities = 34/37 (91%), Positives = 35/37 (94%)
 Frame = -1

Query: 218     PEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTAAGL 254
               PEKWNKIKLVVTQEDVELAYQEAMMNMARLNRT  G+
Sbjct: 1472758 PEKWNKIKLVVTQEDVELAYQEAMMNMARLNRTGTGV 1472648



 Score = 65.1 bits (157), Expect = 6e-09
 Identities = 30/30 (100%), Positives = 30/30 (100%)
 Frame = -2

Query: 188     PDNAVPGDVLVLTKPLGTQVAVAVHQWLDI 217
               PDNAVPGDVLVLTKPLGTQVAVAVHQWLDI
Sbjct: 1474794 PDNAVPGDVLVLTKPLGTQVAVAVHQWLDI 1474705


  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.320    0.135    0.405 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 1,055,734,986
Number of Sequences: 4245
Number of extensions: 15947541
Number of successful extensions: 53022
Number of sequences better than 10.0: 2
Number of HSP's better than 10.0 without gapping: 70
Number of HSP's successfully gapped in prelim test: 83
Number of HSP's that attempted gapping in prelim test: 48210
Number of HSP's gapped (non-prelim): 15333
length of query: 392
length of database: 800,060,566
effective HSP length: 132
effective length of query: 260
effective length of database: 799,500,226
effective search space: 207870058760
effective search space used: 207870058760
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.8 bits)
S2: 78 (34.7 bits)
