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= SPP00000033_1.0 # Protein # Selenophosphate synthetase 2 (SPS2)
# Homo sapiens # Complete
         (448 letters)

Database: N.gruberi/genome.fa 
           784 sequences; 40,964,085 total letters

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

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gb|GG738852.1| Naegleria gruberi genomic scaffold NAEGRscaffold_...    82   1e-15
gb|GG738883.1| Naegleria gruberi genomic scaffold NAEGRscaffold_...    73   6e-13
gb|GG738907.1| Naegleria gruberi genomic scaffold NAEGRscaffold_...    30   6.2  

>gb|GG738852.1| Naegleria gruberi genomic scaffold NAEGRscaffold_8, whole genome shotgun
              sequence
          Length = 938756

 Score = 81.6 bits (200), Expect = 1e-15
 Identities = 89/387 (22%), Positives = 149/387 (38%), Gaps = 27/387 (6%)
 Frame = +1

Query: 50     RLTGFSGMKGUGCKVPQEALLKLLAGLTRPDVRPPLGRGLVGGQEEASQEAGLPAGAGPS 109
              +LT  +   G GCK+  + L  +L+ + +    P +   LVG +                
Sbjct: 697966 KLTSTTHGLGCGCKIRAQTLSSILSKIPKQPFDPNV---LVGTETS-------------- 698094

Query: 110    PTFPALGIGMDSCVIPLRHGGLSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGITE 169
                        D C +      L+L  T DFF P+V+ P   G IACAN LSD+YAMG  E
Sbjct: 698095 ----------DDCGVYKLTSELALCATVDFFTPIVDSPKDFGAIACANALSDVYAMG-GE 698241

Query: 170    CDNXXXXXXXXXXXXXXXREKVTPLMVKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVATV 229
              C                    +  L+++G ++                +     G   + 
Sbjct: 698242 CKVALNIVAFPVKRLSLD---ILQLILQGAQEKASEAGVTIIGGHSIEDNEPKFGMCVSG 698412

Query: 230    VCQPNEFIMPDSAVVGDVLVLTKPLGTQVAVNA-HQWLDNPERWNKVKMVVSREEVELAY 288
              +  P           GDVL+LTKP+GT V   A  + L  P      K +     +  + 
Sbjct: 698413 LVHPERIWKNVGCRNGDVLLLTKPIGTGVLTTAMKKGLIGPNDDCSKKCIEWMSTLNKSA 698592

Query: 289    QEAMFNMATLNRTAAGLMHTFNAHAATDITGFGILGHSQNL----------AKQQRNEVS 338
               E +  +  +N ++      +  HA TDITGFG++GH + +          A+ +  +V 
Sbjct: 698593 SEILRGLHRVNSSST----PYEIHACTDITGFGLMGHLREMITGDNQDNLYAQIESEKVP 698760

Query: 339    FV----------------IHNLPIIAKMAAVSKASGRFGLLQGTSAETSGGLLICLPREQ 382
              F+                I+N+  + K+   +++  +        A+TSGGLL  +   +
Sbjct: 698761 FMDKAYEMILSDIIPGGTINNMNFVQKITKFNESVSQIVRTLLCDAQTSGGLLFSVSPHE 698940

Query: 383    AARFCSEIKSSKYGEGHQAWIVGIVEK 409
              +    S+ K      G   W VG +++
Sbjct: 698941 SEEIISDFKRL----GISIWTVGSIKQ 699009


>gb|GG738883.1| Naegleria gruberi genomic scaffold NAEGRscaffold_39, whole genome
             shotgun sequence
          Length = 371401

 Score = 72.8 bits (177), Expect = 6e-13
 Identities = 60/205 (29%), Positives = 90/205 (43%), Gaps = 1/205 (0%)
 Frame = -2

Query: 131   LSLVQTTDFFYPLVEDPYMMGRIACANVLSDLYAMGITECDNXXXXXXXXXXXXXXXREK 190
             +++V T DFF P+V++PY  G IACAN LSD+YAMG                       +
Sbjct: 75543 VAIVSTIDFFAPIVDNPYDFGAIACANALSDVYAMG----GKPQTALSVVAFPSRTLPIE 75376

Query: 191   VTPLMVKGFRDXXXXXXXXXXXXXXXXNPWIIIGGVATVVCQPNEFIMPDSAVVGDVLVL 250
             +   +++G  D                +     G     +  P++         GD+L L
Sbjct: 75375 ILQEILRGSIDKANEAGVAIVGGHSIDDEEPKFGMSVNGIIHPDKIFRNVGCEEGDLLFL 75196

Query: 251   TKPLGTQVAVNAHQWLDNPERWNKVKMVVSREEVELAYQEAMFNMATLNR-TAAGLMHTF 309
             TK +GT V  +A            +K+     E E A Q A+ +M TLNR T+  L    
Sbjct: 75195 TKRIGTGVITSA------------MKLGNVPYESETA-QSAINSMKTLNRDTSYVLREMK 75055

Query: 310   NAHAATDITGFGILGHSQNLAKQQR 334
             +  + TD+TGFG+LGH + + K  R
Sbjct: 75054 SVSSCTDVTGFGLLGHLKEIMKGSR 74980


>gb|GG738907.1| Naegleria gruberi genomic scaffold NAEGRscaffold_63, whole genome
             shotgun sequence
          Length = 275174

 Score = 29.6 bits (65), Expect = 6.2
 Identities = 28/82 (34%), Positives = 33/82 (40%), Gaps = 3/82 (3%)
 Frame = +2

Query: 37    PFEPQALGLSPSWRLTGFSGMKGUGCKVPQE-ALLKLLAG--LTRPDVRPPLGRGLVGGQ 93
             PF    LGL P        G  G GC      +LL  + G  L  P  +PP G G +GG 
Sbjct: 65810 PFGGGTLGLLPK------DGSDGGGCGDSSC*SLLDGIGGAVLPPPPPKPPRGCGTLGG- 65968

Query: 94    EEASQEAGLPAGAGPSPTFPAL 115
                 +E G   G G  P  P L
Sbjct: 65969 --CGREEGTTGGCGFEPKPPPL 66028


  Database: N.gruberi/genome.fa
    Posted date:  Feb 9, 2011  9:02 AM
  Number of letters in database: 40,964,085
  Number of sequences in database:  784
  
Lambda     K      H
   0.320    0.136    0.412 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 15,333,426
Number of Sequences: 784
Number of extensions: 190779
Number of successful extensions: 579
Number of sequences better than 10.0: 3
Number of HSP's better than 10.0 without gapping: 13
Number of HSP's successfully gapped in prelim test: 8
Number of HSP's that attempted gapping in prelim test: 547
Number of HSP's gapped (non-prelim): 59
length of query: 448
length of database: 13,654,695
effective HSP length: 105
effective length of query: 343
effective length of database: 13,572,375
effective search space: 4655324625
effective search space used: 4655324625
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: 63 (28.9 bits)