FASTA-SWAP or FASTA-PAT Search Results

FASTA-SWAP or FASTA-PAT (FASTA-based Pattern database search tools) are modified versions of the FASTA sequence database search tool. See the FASTA-SWAP and FASTA-PAT Help Page for a detailed program description.

Istvan Ladunga, Brent A. Wiese, and Randall F. Smith (1995)
Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030.
Email


 FASTA-SWAP searches a protein pattern database
 FASTA-SWAP version 1.0, Dec. 1995, 
FASTA version 2.0u August, 1995
Please cite: I. Ladunga, B. Wiese & R.F. Smith (1996), submitted
 W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448

 /tmp/fastapat.seq.26533 :  345 aa

 >GI||SP|Q|DBP_HUMAND-BINDINGPROTEIN(DBP)(ALBUMINDB
 X-BINDINGPRTEINTAXREBMARPVSDRTPAxxxxxxxxxxxxxxxxxxGLRSLLQGTS
 KPKEPASCLLKEKERKxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
 xxxxxxERTLPFGDVEYVDLDAFxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
 xxxxxxxxxxxxxxxxGTASGHRAGLTSRDTPSPVDPDTVEVLMTFEPDPADLALSSIPG
 HETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRYKNNEAAKRSRDARRLK
 ENQISVRAAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL

 vs EC Database library
 searching /local/dot5/sl_home/beauty/seqdb/comb-db/ec 11 library


 one = represents 1 library sequences

   z-opt E()
< 20     0     0 :
  22     1     0 :=
  24     2     0 :==
  26     2     0 :==
  28     2     0 :==
  30     2     1 :*=
  32    10     4 :===*======
  34    11    11 :==========*
  36    21    22 :=====================*
  38    30    36 :==============================     *
  40    48    51 :================================================  *
  42    44    62 :============================================               *
  44    35    68 :===================================                        *
  46    28    70 :============================                               *
  48    31    67 :===============================                            *
  50    30    61 :==============================                             *
  52    45    53 :=============================================       *
  54    46    46 :=============================================*
  56    52    38 :=====================================*==============
  58    51    31 :==============================*====================
  60    44    25 :========================*===================
  62    44    20 :===================*========================
  64    37    16 :===============*=====================
  66    31    13 :============*==================
  68    21    10 :=========*===========
  70    21     8 :=======*=============
  72    11     6 :=====*=====
  74    14     5 :====*=========
  76     9     4 :===*=====
  78     5     3 :==*==
  80     5     2 :=*===
  82     3     2 :=*=
  84     1     1 :*
  86     1     1 :*
  88     0     1 :*
  90     0     1 :*
  92     0     0 :
  94     1     0 :=
  96     0     0 :
  98     0     0 :
 100     0     0 :
 102     0     0 :
 104     0     0 :
 106     1     0 :=
 108     0     0 :
 110     0     0 :
 112     0     0 :
 114     0     0 :
 116     0     0 :
 118     0     0 :
>120     0     0 :
 850743 positions in  1625 patterns
 statistics extrapolated from 740 to 740 patterns
 Kolmogorov-Smirnov statistic: 0.2202 (N= 28) at 52
 results sorted and z-values calculated from opt score
  740 scores better than 1 saved, ktup: 1, variable pamfact
 gap penalties: -100,-10
 joining threshold: 83, optimization threshold: 50, width: 32
  scan time:  0:02:01
The best scores are:                                     opt  z-sc  E(740)
942 glycerol-3-phosphate sn-glycerol-3-phosphate acylt   761 106.5  0.30
9_3 protein-tyrosine 1c 2c corkscrew phosphatase 3.1.3   634  93.8  1.50
573 oleoyl-acyl oleoyl-acyl-carrier carrier protein hy   567  85.6  4.35
1090 type ii deoxyribonuclease enzyme ec 3.1.21.4 endo   559  83.9  5.35
2_16 myosin myosin-light-chain light chain ec kinase 2   464  83.1  5.93
1063 glycosaminoglycan n-heparan n-acetylglucosaminyl    486  83.0  6.02
16_10 dipeptidyl-peptidase i cathepsin ec precursor c    524  82.8  6.22
295 rna polymerase alpha subunit nucleocapsid phosphop   514  81.4  7.39
12_3 glyceraldehyde glyceraldehyde-3-phosphate 3-phosp   519  80.6  8.18
483 aspartate-semialdehyde dehydrogenase ec asa 1.2.1.   514  80.6  8.20
828 protein-pii uridylyltransferase pii uridylyl-trans   450  80.0  8.83
3_23 interferon-induced serine/threonine-protein doubl   482  79.7  9.15
846 1l-myo-inositol-1-phosphate d-myo-inositol-3-phosp   475  79.4  9.57
712 aspartate-semialdehyde dehydrogenase asa ec 1.2.1.   498  79.3  9.69
614 flavonol udp-glucose 3-o-glucosyltransferase bronz   482  79.2  9.73
942 glycerol-3-phosphate sn-glycerol-3-phosphate acylt (828 aa)
initn: opt:  761 z-score: 106.5 E():    0.3
Smith-Waterman score: 778;    14.8% identity in 358 aa overlap
 19.0% noncontradicting positions,   4.2% class identity

               10        20        30        40        50        60
GI||SP XBINDINGPRTEINTAXREBMARPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSK
                             R  +  T                     L  L +  S 
942                   MTFCYPCRAFALLTRGFTSFMSGWPRIYYKLLNLPLSILVKSKSI
                              10        20        30        40     

               70                  80        90       100       110
GI||SP PKEPASCLLKEKER----------KXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
       P +PA  L  +  R----------                                    
942    PADPAPELGLDTSRPgIMYVLPYNSKADLLTLRAQCLAHDLPDPLEPLEIDGTLLPRYVF
                      s                                            
          50        60        70        80        90       100     

              120       130              140       150       160   
GI||SP XXXXXXXXXXXXXXXERTLPFGDVEYVDL-------DAFXXXXXXXXXXXXXXXXXXXXX
                      E    F D--Y DL-------D+                      
942    IHGGPRVFTYYTPKEESIKLFHD--YLDLHRSNPNLDVQMVPVSVMFGRAPGREKGEVNP
         110       120         130       140       150       160   

           170       180       190       200       210       220   
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPVDPDTVEVLM
                                         A  H    T     + V------  M
942    PLRMLNGVQKFFAVLWLGRDSFVRFSPSVSLRRMADEHGTDKTIAQKLARV------ARM
           170       180       190       200       210             

           230       240       250       260       270       280   
GI||SP TFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRY
        F----A   L+++    ---P R+      L    I K -----V +E +  K    + 
942    HF----ARQRLAAVGPRL---PARQDLFNKLLASRAIAKA-----VEDEARSKKISHEKA
           220       230          240       250            260     

                  290       300                  310       320     
GI||SP KNNE-------AAKRSRDARRLKE-------NQ----ISVRAAFLEKENALLRQEVVAVR
       + N -------AA  S +  RL +-------N ----I+V  A   ++ A   +E V V 
942    QQNAIALMEEIAANFSYEMIRLTDRILGFTWNRLYQGINVHNAERVRQLAHDGHELVYVP
         270       280       290       300       310       320     

         330       340                                             
GI||SP QELSHYRAVLSRYQAQHGAL                                        
          SH    L  Y   H  L                                        
942    CHRSHMDYLLLSYVLYHQGLVPPHIAAGINLNFWPAGPIFRRLGAFFIRRTFKGNKLYST
         330       340       350       360       370       380     

9_3 protein-tyrosine 1c 2c corkscrew phosphatase 3.1.3 (768 aa)
initn: opt:  634 z-score: 93.8 E():    1.5
Smith-Waterman score: 634;    18.5% identity in 65 aa overlap
 35.4% noncontradicting positions,  16.9% class identity

              220       230       240       250       260       270
GI||SP PSPVDPDTVEVLMTFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVP
                        d    +l      Et D r+  F EE    Q    k -  q  
9_3    aflqLrQPlnaTRinAAdInaRVrqLnKlaEseDkaKaGFWEEFEsLQqQdcrn-lhqRn
       tvvh k  fyt  vt  e en  ee s kq tt tr q       t  k eskd tfs l
          y    y        g  s  l           v g             v l  y  k
      210       220       230       240       250       260        

              280       290       300       310       320       330
GI||SP EEQKDEKYWSRRYKNNEAAKRSRDARRLKENQISVRAAFLEKENALLRQEVVAVRQELSH
       E qk E     RYkN      sR      e  i v          l            ==
9_3    EGqrqENrnKNRYrNILPfDHsRViLqdrdhneagadYINANiirnqldgddnNak  ==
         ykp  kg    k    y  t  k hggepsipvse     yvklefetkce ms    
               l               v l v s v           mpptl ply s     
               s                                                   
       270       280       290       300       310       320       

              340                                                  
GI||SP YRAVLSRYQAQHGAL                                             
       ===============                                             
9_3    ============================================================

573 oleoyl-acyl oleoyl-acyl-carrier carrier protein hy (398 aa)
initn: opt:  567 z-score: 85.6 E():    4.3
Smith-Waterman score: 680;     9.0% identity in 277 aa overlap
 26.0% noncontradicting positions,  17.0% class identity

                                    10        20        30         
GI||SP                      XBINDINGPRTEINTAXREBMARPVSDRTPAXXXXXXXX
                                     rt i    r      vs rtp       = 
573    SrgaP==AAaAtaaknnanaknnrgalphfrdigqarirrrrngairrrqarkraraP=l
        kpl     p vmlmlscsgtdhlhtfsfskspilciflnsdsiflcnsplsspkld  t
                   ttty   vv stfvyt  stssffspvk stlsvss spv  tvs  v
                             t         v sv   t y v     ts   w     
               10        20        30        40        50          

      40        50        60          70        80        90       
GI||SP XXXXXXXXXXGLRSLLQGTSKPKE--PASCLLKEKERKXXXXXXXXXXXXXXXXXXXXXX
                 g    L    k  E--==Sc l ek R                       
573    raaasadqnaqsagnLrqadaEiEe===ScslanrlRaGrLmEDGlSYKErFiiRcYEVG
       lsvrtnevgsgvlvg keeqk v k    stpgdkf l s t   y    k vv s    
       m  nygg pti p s    e             e   m                      
       f  i          v                                             
      60        70        80           90       100       110      

       100       110       120       130       140       150       
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXERTLPFGDVEYVDLDAFXXXXXXXXXXXXXXX
                                   ------------- D F               
573    INKTATiETIANLLQEVacNHaQgcGFS-------------TDGFATTlTMRKLHLIWVT
             v          gg  v kv                       p           
                              s                        t           
        120       130       140                    150       160   

       160       170       180       190       200       210       
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPV--D
                                               Asg   G---R T   V--+
573    aRMHIEIYrYPAWSDViEIETWcQgEGriGTRRDWIirDhAnnEVIG---RATSKWVMMN
       s       k       v     v s  kv       lk f sd                 
                                           m  s tg                 
                                              y                    
           170       180       190       200       210          220

         220       230       240       250       260       270     
GI||SP PDTVEVLMTFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKD
        DT-  L     d  d  L------- F Pr  R s  E-k    mKK  K++ P e   
573    qDT-RRLQrVnDdVRdEYL-------VFCPrelRLaFPE-eNnnSlKKIaKLEDPAqhSr
       e       k s e  e              ktp  s    k tk m   p      ey m
                 t                                s     s         t
               230              240       250        260       270 

         280       290       300       310       320       330     
GI||SP EKYWSRRYKNNEAAKRSRDARRLKENQISVRAAFLEKENALLRQEVVAVRQELSHYRAVL
            RR   +     +  a      +          E   +  +     Q      +  
573    LeLkPRRaDLDMNqHVNNVaYIGWaLESIPqEiidTHELqaITLDYRRECQrDDiVDSLT
        g s   s     k     t    v     k vly    et          q  v     
          v                          p         v          h        
             280       290       300       310       320       330 

1090 type ii deoxyribonuclease enzyme ec 3.1.21.4 endo (333 aa)
initn: opt:  559 z-score: 83.9 E():    5.4
Smith-Waterman score: 567;    13.5% identity in 327 aa overlap
 21.7% noncontradicting positions,   8.3% class identity

                  10        20        30        40        50       
GI||SP  XBINDIN--GPRTEINTAXREBMARPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQG
           N In-- P+ +i+ A      R v--------                 L   L  
1090     mpNKInFPePKqQiDFAFaLKRFRai--------YLQnALLETVRDMDIaaLDrQLae
         tt   v  p  l v    s     gv           s           ve  t  sk
                 10        20                30        40        50

        60        70        80        90       100       110       
GI||SP TSKPKEPASCLLKEKERKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
          p + As L     R----                                       
1090   YVnkADLAgTLAqYGLR----AELLFaVPalLEaNPhLLGYYRLLMGYSQKEFYGrDKGF
         pp    s   t             p  vv  t  f                  k    
               60            70        80        90       100      

       120       130       140           150       160       170   
GI||SP XXXXXXXXERTLPFGDVEYVDLD----AFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
               E     G v   dl ----AF                               
1090   naGCFKSMEeKGnaGKaAkadihdLChAFCatASaLLaGVGPLRiSRELLDDLTLLTVGP
       gv       v  ki  v vpklse  v   sv  s  q      v               
        110       120       130       140       150       160      

           180       190       200       210        220       230  
GI||SP XXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPVDPDT-VEVLMTFEPDPADL
                                  H-A    R++   V+  T-  VL+ F PDP-D+
1090   QLRGGANNQRGADGIVLVFEIIrEIVaH-AaAEVREnAIEVNSATGRnVLIEFAPDP-DI
                             k   s   v     s          p            
        170       180       190        200       210       220     

            240       250       260       270       280        290 
GI||SP ALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRY-KNNEAAKR
        + -----------------EE+  q----  R +-V  E K     S  +-+  EA K 
1090   IIR-----------------EEMend----nYRNV-VAIEVKSGTDVSNIHNRIGEAEKS
                              slq    h                             
                           230            240       250       260  

               300       310       320       330       340         
GI||SP SRDARRL--KENQISVRAAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL    
          AR l-- E    V +  L+   A--R E  +     S     L              
1090   HQKARqhGfTECWTVVNVgRLDMdKA--RkESPSTnRFYSiTaLSLRaGdqYaDFRRRVL
            gl y         s    v     s     d    l d    e ee e       
            270       280         290       300       310       320

                  
GI||SP            
                  
1090   SLTAIpAaPan
            s r tp
              330 

2_16 myosin myosin-light-chain light chain ec kinase 2 (1718 aa)
initn: opt:  464 z-score: 83.1 E():    5.9
Smith-Waterman score: 756;     9.3% identity in 343 aa overlap
 25.9% noncontradicting positions,  16.6% class identity

                                           10            20        
GI||SP                             XBINDINGPRTEIN----TAXREBMARPVSDR
                                     ++D  g  tE+n----T   E  ---VSD 
2_16   VYGiSEPSQESEltkVGEKqEEeLKdEeaElSDDdeKEpEVnYrTVTiNTEQK---VSDf
          t        vvt    p  p  e ve v   eg  t  d q   v           v
             1210      1220      1230      1240      1250          

       30        40        50        60        70                  
GI||SP TPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSKPKEPASCLLKEKE--------------R
                                 l g    k  a----KeKE--------------r
2_16   mnieEaLGgGKFGaVcrctErkTGlilAaKfiKaqsa----KdKEnirdEIeiMNcLnHr
       ydsk r  s    q ftlv ks  kkw g vf kytp     e  mvlq  gv  q h p
                                v                     ml  s        
      1260      1270      1280      1290          1300      1310   

           80        90       100       110       120       130    
GI||SP KXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXERTLPFGDV
       k                                                  e  l     
2_16   nLiQcyaAiEekanIilflEiieGGELFERIiDEDfqLTErdcikfmRQIcdGiefiHKq
       k v lvd f tphe vmvm mvs        v   ye   vetmvyv   se vlym  m
             s    s        y               h                       
          1320      1330      1340      1350      1360      1370   

          140       150       160       170       180       190    
GI||SP EYVDLDAFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
         v LD                                                      
2_16   rilHLDLKPENIlCVNkTGhriKiIDFGLARRlnnaekLKVnFGTPEFlaPEViNYdqIg
       gvv         m   t  tlv l        yepngs   l      vs   v  ep s
                           s             s                         
          1380      1390      1400      1410      1420      1430   

             200        210       220       230       240          
GI||SP XGT---ASGHRA-GLTSRDTPSPVDPDTVEVLMTFEPDPADLALSSIPGHETF----DPR
         T---+ G   - L S  +P   D Dt-E L        d ------  EtF----D  
2_16   daTDMWSiGVIcYiLlSGLSPFlGDnDn-ETLaNVlSanWdF------DdEaFdaiSDdA
       ye     l   t m v      m  d t    n  t gt y        e t eev  e 
        k                                                          
          1440      1450      1460       1470            1480      

        250        260       270          280       290       300  
GI||SP RHRFSEEELKPQ-PIMKKARKIQVPEEQ---KDEKYWSRRYKNNEAAKRSRDARRLKENQ
       +   S    K q-  m  a  +q P  q---kd K   rr k+    K+    Rr k n 
2_16   KDFiSNLikKdqgaRlnaaQCLaHPWLnNLAedaKrcnarllSqdrlKKYlarRrWqKng
          v   lt emkn mdce   q    q   kkt nmerkkk kilm   mmk k k tf
               v   ss  s t        m                                
       1490      1500      1510      1520      1530      1540      

            310       320       330       340                      
GI||SP ISVRAAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL                 
       i+VrAa   k  a             S    + +       ==                 
2_16   nAVrAanRlkkiamigalmargaSggSPTSPinAdrlEnEd====AFLEaVAEEKPHVKP
       h  s ig fssmssssgmsglks ts     lt ekv t e        e          
       i                     v                                     
       1550      1560      1570      1580          1590      1600  

1063 glycosaminoglycan n-heparan n-acetylglucosaminyl  (888 aa)
initn: opt:  486 z-score: 83.0 E():      6
Smith-Waterman score: 535;    14.2% identity in 134 aa overlap
 26.1% noncontradicting positions,  11.9% class identity

        170       180       190       200       210       220      
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPVDPDTVEVLMTFE
                              ====   + Sg  ag      p==v p  v    t  
1063   FcLfSmgflAYYlsgsnralEPl====aDaSeSdagadPgPaaPp==lPirPqrapapsR
        s v vfvs   vytwpkgk  s    g c s gcsgg p vr s  v pk vqpvett 
          30        40            50        60          70         

        230           240          250       260             270   
GI||SP PDPADL----ALSSIPGHE---TFDPRRHRFSEEELKPQPIM------KKARKIQVPEEQ
        dP  L----+l S  G+E---  +  R r+s E       M------   R + v  E 
1063   TdPlVLVFVESaYSQLGQEiVAILESSRFrYrTEiAPGrGDMPTLTDhgrGRfaLiIYEN
        e v       l       v         k s  l   k        kth  yv v    
      80        90       100       110       120       130         

           280       290       300       310       320       330   
GI||SP KDEKYWSRRYKNNEAAKRSRDARRLKENQISVRAAFLEKENALLRQEVVAVRQELSHYRA
       ------  +Y N +A  R    r   e  + +   F  +En+LL  +       L     
1063   ------iLKYVNLDAWnRELLDrYCVaYGVGIIGFFrArEnSLLSAQLKGFPLFLHSNLG
             l         s     k   e         k n h                   
           140       150       160       170       180       190   

           340                                                     
GI||SP VLSRYQAQHGAL                                                
                  L                                                
1063   LrDcqiNPkaPLLhlTRPSrlEkGpLPGdDWTiFQSNHSTYEPVLiAkhRpaE=SiP===
        k ysv  ss   yv    ev p v   e   v            l st ss   m    
           200       210       220       230       240             

16_10 dipeptidyl-peptidase i cathepsin ec precursor c  (468 aa)
initn: opt:  524 z-score: 82.8 E():    6.2
Smith-Waterman score: 595;    14.9% identity in 355 aa overlap
 21.7% noncontradicting positions,   6.8% class identity

               10        20        30        40        50        60
GI||SP XBINDINGPRTEINTAXREBMARPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSK
                            +  VS  TPA                    RS +  +--
16_10     MGPWTHSLRAALLLVLLGVCTVSSDTPANCTYPDLLGTWVFQVGPRHPRSHINCS--
                  10        20        30        40        50       

               70            80        90       100       110      
GI||SP PKEPASCLLKEKE----RKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
         EP ---l EK ----+K                                         
16_10  VMEPT---dEEKVVIHLKKLDTAYDEVGNSGYFTLIYNQGFEIVLNDYKWFAFFKYEVKG
               h                                                   
               l                                                   
          60           70        80        90       100       110  

        120       130       140       150       160       170      
GI||SP XXXXXXXXXERTLPFGDVEYVDLDAFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
                  T    Dv   +   F                                  
16_10  SRAISYCHETMTGWVHDyLGRNWACFVGKKMANHSEKVYVNVAHLGGLQEKYSERLYSHn
                        v                                         h
            120       130       140       150       160       170  

        180       190                  200       210       220     
GI||SP XXXXXXXXXXXXXXXXXXXGT-----------ASGHRAGLTSRDTPSPVDPDTVEVLMTF
                            ----------- SGH-+G   R  P+P+---T E+    
16_10  HNFVKAINSVQKSWTATTYrrYEKLSIRDLIRRSGH-SGRILRPKPAPI---TDEIQQQI
                          ee                                       
            180       190       200        210       220           

         230       240       250       260               270       
GI||SP EPDPADLALSSIPGHETFDPRRHRFSEEELKPQPI--MKKAR------KIQVPEEQKDEK
          P      ++ G     P R+  S         --M  AR------  Q P     E 
16_10  LSLPESWDWRNVRGINFVSPVRNQESCGSCYSFASiGMLEARIRILTNNSQTPILSPQEV
                                          l                        
      230       240       250       260       270       280        

       280               290       300          310       320      
GI||SP YWSRRYKNN--------EAAKRSRDARRLKEN---QISVRAAFLEKENALLRQEVVAVRQ
            Y   -------- A K + D     EN---  +  A    KEN L--------R 
16_10  VSCSPYAQGCDGGFPYLIAGKYAQDFGVVEENCFPYTATDAPCKPKENCL--------RY
      290       300       310       320       330               340

        330       340                                              
GI||SP ELSHYRAVLSRYQAQHGAL                                         
         S Y  V   Y   + AL                                         
16_10  YSSEYYYVGGFYGGCNEALMKLELVKHGPMAVAFEVHDDFLHYHSGIYHHTGLSDPFNPF
              350       360       370       380       390       400

295 rna polymerase alpha subunit nucleocapsid phosphop (422 aa)
initn: opt:  514 z-score: 81.4 E():    7.4
Smith-Waterman score: 514;     4.3% identity in 70 aa overlap
 42.9% noncontradicting positions,  38.6% class identity

      170       180       190       200       210       220        
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPVDPDTVEVLMTFEPD
                                     sg  aG t     Sp  pd   vl t  p 
295    aerqqar araapgnrdidangalaGqeerprlqaGaihaahrSrqqqanlnanaddaqa
       qtiienq sniqspdqgtqedirfh lkiesqsls vtlgvlq nlllddklvhvehvpq
        wvp e  vkt  sl s egstws  smmk gyp    ph pl gpspgssp l gi  l
          t g        f t gs  vy   t   k s     p    l   l t  p tp  k
                     x            x           x    x   x    x  x  x
          190        200       210       220       230       240   

      230       240       250       260       270       280        
GI||SP PADLALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRYKNNEA
       pad al  i   e  d R  r ---e k   imk    iq           s        
295    lannanaiinairaidaRlnri---dqqlDrilaqqaliqnirndiqqiKaglAtiEani
       fvdfvqemldllenleq mdql   egkv lvmlsgnmviqlkselsgl gsm vl gll
       p es l  memmkgmsl vqhv    t   k  ktsss pmm   mtsv tts  m  mm
       s iy k   stv    s  sk     y      s t t ttt   tvt   vv  v    
         x  x          x   x                                       
           250       260          270       280       290       300

      290       300       310       320       330       340        
GI||SP AKRSRDARRLKENQISVRAAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL   
       a           n             a l  e   v    s    +     aq  a    
295    amiKIlDPgnannsplrdiraearnhaiiiaengndgkladqdneialneLarPiahgre
       gtm  m  scpsgvsvdelnrgledvptlvpgpfdgpffsqeedllsmdk sq tqkpsk
       t v      t tk t esv ksfhq vvv s stspspstegggmmt  p    vlsstp
                v  v   s    tpss     t t  v  tvgssit y  s     pt v 
                            xx               x x xxx          x    
              310       320       330       340       350       360

12_3 glyceraldehyde glyceraldehyde-3-phosphate 3-phosp (344 aa)
initn: opt:  519 z-score: 80.6 E():    8.2
Smith-Waterman score: 519;     4.7% identity in 211 aa overlap
 26.5% noncontradicting positions,  21.8% class identity

               10        20        30        40        50        60
GI||SP XBINDINGPRTEINTAXREBMARPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSK
            +   +  iN   R--+ R V =R               =      l  llq  S 
12_3        laklKVGiNGFGR--IGRLVl=RaaenngnidiVaIND=lidlnnlaYllqhDSn
            mgms   v            t  lgivkpkvef g    ftppdimv mfky  t
               v                    lf ss   v       vs ly          
                    10           20        30         40        50 

               70        80        90       100       110       120
GI||SP PKEPASCLLKEKERKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
            s   k ke k                                             
12_3   hGklrgdqVdaendgiiinGnfIacfaernaanikWgaaGAdiiVEcTGlfTdmEgaqaH
       y pfhshs efkdgklvvd hp piyqvkdpkelp kdl  eyv  s  fy sp ksgk 
           pwlt k te s     ks tvvs  e sk    kk          v  ts   s  
              60        70        80        90       100       110 

                 130       140       150       160       170       
GI||SP XXXXXERTL---PFGDVEYVDLDAFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
             r l---P gd   v                                        
12_3   LdaGAKrViISAPaadadRmrmilmnVNddllDnnldtIiSnASCTTNCLAPiAKainDn
        qg   k l    sgepp vktlvyg  hekf gskkv v v          l  vlh s
                    tk       f v     ty p                  m       
             120       130       140       150       160       170 

       180       190       200         210       220       230     
GI||SP XXXXXXXXXXXXXXXXXXGTASGHRAGL--TSRDTPSPVDPDTVEVLMTFEPDPADLALS
                          t+ G r  l-- sR +   + P t     +      +L   
12_3   FGieeGlMTTiHAiTaTQklVDGPrgKdlRagRaAaqNIIPasTGAAKAialViPELdGK
         ltv t   v  m g  pt    sk lw ds g le    ht      vgk l   k  
          v         y    s           g          s               s  
             180       190       200       210       220       230 

         240       250       260       270       280       290     
GI||SP SIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRYKNNEAAKRSRDA
                       e   k     k      v  +              e    s   
12_3   LkGhAqRVPiangSVtdLTcrlaeaakaddinaalKQAsenneaGIlGYTdeqiVSsDin
        t m f   tkdv  ve  fkfekktsyeevcnvm   ttgplk  f   ehev  t fq
                vpt       y ts pvt k  kk v      s s               i
             240       250       260       270       280       290 

483 aspartate-semialdehyde dehydrogenase ec asa 1.2.1. (403 aa)
initn: opt:  514 z-score: 80.6 E():    8.2
Smith-Waterman score: 514;     2.2% identity in 274 aa overlap
 37.6% noncontradicting positions,  35.4% class identity

               10        20        30        40        50        60
GI||SP XBINDINGPRTEINTAXREBMARPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSK
             n        a rebm     drt----         =      G  ==== t k
483     MarqlnaailGATGagrqrmrnqLadrn----headelrll=aSeRaaGrk====irfr
         gqgkhigvf    qvgeqfiql qqse    lphfkvkff s k se kt    lnwn
         sgkfiv  v    s  tt let ee t    f mts fv    p    tv    yqyk
          mmtt            v  lv           v t t     s          ve  
            x                x                                  x  
                10        20            30         40            50

               70        80        90       100       110       120
GI||SP PKEPASCLLKEKERKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
        ke   --  e e  =====                                        
483    dqdird--pqnaeda=====deddleacnialfgadadaaaaqakraaaaGaaiidNaka
       qkeleh  veditqi     ekegfkdidvvfsslggelsrqiepefmer iivvs ssn
       gttvl     e  tt     kpffwsqvh        sm qkl  i vk  vt    t h
           t     l   v     ssss  e          ty gvf  l      v       
                 x         x     x           x x x  x              
                 60             70        80        90       100   

              130       140       150       160       170       180
GI||SP XXXXXERTLPFGDVEYVDLDAFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
            e  L    V   dldaf       ==========                     
483    fRrndniPLiVpeVNaaaraaranrLaka==========IIanpNCsTaaalaaLePird
       w qdqdt  v sv  eedidefrqn dtn            ccs  t iqlmpp k lhq
       y keeev        pfhlhivghk v g             i     mgmvvv    ie
         m p           s k v  ks                                 yk
         x x           x x                                         
           110       120       130                 140       150   

              190       200       210       220               230  
GI||SP XXXXXXXXXXXXXXXGTASGHRAGLTSRDTPSPVDPDTVE------VLMTFEPDP--ADL
          =            + sG  a  t      pv  d ve------ lm  epdp--ad+
483    aaG=inaintsslQaiSGaGnaainelarqaaaigndnedleadeaaadaaagnakgadI
       ef  ldrlhvtty sv  s qegvdttygeirelidggvqfkeigqdiqpgdpdkfhqq 
       kw   qkvi           is  e  vk tqkvl hv epvhvhpelmypetqppsip 
        y   e  t           l   k   s vkp v k        s  f sgvgvyvs  
        x   x              x                           x xx     x  
            160       170       180       190       200       210  

            240       250       260       270       280       290  
GI||SP ALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRYKNNEAAKRS
       a + ip    f       + eelkp    k===Ki  pee k      R       a   
483    adNaiaqadneqDndegETdedqalraEdr===KilndediqVaaqCnRlaiqdgHaeai
       lf clpligqle dgsf  kkelkmin sq    lmeepell ngh v vpvlfs slcv
        y iv f  gfm kmtw  f im plw tk     fgi s k spt      fs  t s 
          v  y  kyt   y   y     t          hl   m          st    t 
                x x             x          xx   x           x      
            220       230       240          250       260         

            300       310       320       330       340            
GI||SP RDARRLKENQISVRAAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL       
           = k      ra  l ke ==============s    v    q     l       
483    nari=nraaadadraqaaiaaa==============agaelqrdadrqiqP qaaaaagg
       hief dqddipiedikncleqf              dtihvedqikqplf  ldrnrdsi
       sl t skpppsv qtmdi get              q vk legpsh  y  mpld g k
       yv      v t  ev el k y              p  t f v t       v   v s
                       xx x                x  x x                  
     270        280       290                     300        310   

828 protein-pii uridylyltransferase pii uridylyl-trans (902 aa)
initn: opt:  450 z-score: 80.0 E():    8.8
Smith-Waterman score: 473;     6.6% identity in 290 aa overlap
 24.1% noncontradicting positions,  17.6% class identity

            30        40        50        60        70        80   
GI||SP PVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSKPKEPASCLLKEKERKXXXXXXXX
           R                         L  = s  kep+  ++k+  R         
828    LWriRfALHllagRadnRLLFDrQrqiAqrLnY=dgdaneaVERmMqdYfRVtrriaELn
         kv y   mvvs yed     l lsv gl g  ssegklp   f kk y  vlgvs  t
             280       290       300        310       320       330

            90       100       110       120       130       140   
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXERTLPFGDVEYVDLDAFX
                                                  rt  ------  L  F 
828    qlllQhFdEaILagdAaerprPindrFQlRdglidlrdaglFqrqPe------AiLriFy
       emvf l e v  plp dgkve ldee  v gtylevthetv met s       l em v
              340       350       360       370             380    

           150       160       170       180       190       200   
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRA
                                                           g    hR 
828    llarrpaIrGiradTiRqLRdaRrhiddalrrdPeargLFieilrhrqairRnLrrMnRh
       tmvhnse t vyst l l  hh yllqqpfcni hntk  lslfkspggvh g lp h y
          390       400       410       420       430       440    

           210       220       230       240       250       260   
GI||SP GLTSRDTPSPVDPDTVEVLMTFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIMKK
       g+  r  P---------------------  s I G+  fD  -H ++ +e    -+mkk
828    giLGaYlP---------------------qfgHIVGQMQhDLF-HaYTVDaHTir-liln
       sv  r m                      ews       f     i    e  ln vmkk
          450                            460        470        480 

           270           280                290        300         
GI||SP ARKIQVPE-EQKD---EKYWSRRYKN---------NEAAK-RSRDARRLKENQISVRAAF
        rk   pE- qk --- k w r  k --------- + aK-R  D   L-  q  v  a 
828    LrklaseElaqrhPLasdlierLPkpELIfiAaLfHDIaKGRGGDHSeL-GAaDalaFaq
        esffwp treky  cvkvwpk  st   yl g y   g        i    q vvh ce
             490       500       510       520       530        540

     310       320        330       340                            
GI||SP LEKENALLRQEVV-AVRQELSHYRAVLSRYQAQHGAL                       
       l   na   Q vv- VRq L        r      ++                       
828    rHdLnardTQLiaWLVRnHLLMStTAQRrDiqDPqVIkdFAqeVrdqnrLrYLycLTVAD
       l g psme   vv    q     v    k ls  e  fq  ek qtety d  vv     
              550       560       570       580       590       600

3_23 interferon-induced serine/threonine-protein doubl (554 aa)
initn: opt:  482 z-score: 79.7 E():    9.1
Smith-Waterman score: 580;     9.3% identity in 259 aa overlap
 27.8% noncontradicting positions,  18.5% class identity

          30        40        50        60           70            
GI||SP SDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSKPKEP---ASCLLKE---KERKXXXX
         +                           +G SK kE ---A+ L  e---kE K    
3_23   kYqELpnSGPPHDRRFTFQViIDerEFgEaeGRSKhkEArNAAAKLAVdILnnEnKadch
       t k  st             l  gk  p gk    lt  k        e  dk k vvsp
          30        40        50        60        70        80     

      80        90       100       110       120        130        
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXER-TLPFGDVEYVD
       ====                                            -    G  eY  
3_23   ====TsaceqGLfmGNYIGLiNriAQKKrLsVNYEQCapnsegPqrFhcKCKiGQkeYgi
            tnsse  sv      v sf    k t      esgvhl eg iy   m  tm st
              90       100       110       120       130       140 

      140       150       160       170       180       190        
GI||SP LDAFXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTA
                             ======         ========             ++
3_23   GsGsTKQEAKQLAAKeAYqqiL======KSdplsaGs========SnSlVTSTlageSSS
        t v           l  lkl         pykts t         s v    fssf   
             150       160             170               180       

      200       210       220       230        240                 
GI||SP SGHRAGLTSRDTPSPVDPDTVEVLMTFEPDPADLALSSI-PGHETFDPRRHR--------
       sg  +  tS == +p    +  vlm    +    a  s+-P==== D +r +--------
3_23   egdfSadtSq==SapdSlnSenllmNGLrnNqRKagrklaP====dDmqrnKYTlDaRFn
       ssmt ngv e   nsg fs ssvft   ge k  skvsvs     p vket   v k  g
       190         200       210       220           230       240 

       250       260       270       280        290       300      
GI||SP --FSEEELKPQPIMKKARKIQVPEEQKDEKYWSRRYK-NNEAAKRSRDARRLKENQISVR
       --F e El       ++ K +   + K--+Y  +R K-NnE A r   A         V 
3_23   mDFedIEeIGlGGFGQVFKAKHRIDGK--rYaIKRVKYNnEKAErEVqALAeLnHVNIVq
       s  ke  l  s                  t v       t    h  k   k d     h
             250       260         270       280       290         

        310       320       330       340                          
GI||SP AAFLEKENALLRQEVVAVRQELSHYRAVLSRYQAQHGAL                     
                   +  ============= s                              
3_23   YngCWdGfDYDPEhS=============mknsSRsKTrCLFIQMEFCDKGTLEQWirnRr
        hs  e v     t              ssdt  y  k                 mek n
     300       310                    320       330       340      

846 1l-myo-inositol-1-phosphate d-myo-inositol-3-phosp (542 aa)
initn: opt:  475 z-score: 79.4 E():    9.6
Smith-Waterman score: 494;     8.1% identity in 186 aa overlap
 23.1% noncontradicting positions,  15.1% class identity

             30        40        50        60        70        80  
GI||SP RPVSDRTPAXXXXXXXXXXXXXXXXXXGLRSLLQGTSKPKEPASCLLKEKERKXXXXXXX
            r                        Sl Q- S  k   s ---e e         
846    LTagpranRegimWsfqeKAKGVqQaNYFGSlTQ-cSsirlGIda---nGndiYAPFnSL
         lwviyw ists vtkt     k p     m   s tlkv  sf   e eev    k  
          90       100       110        120          130       140 

             90       100       110       120       130       140  
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXERTLPFGDVEYVDLDAF
                                                  e   P   + y D  A 
846    LPMVnPdeiVfgGWDInnanLadAMqRaqVLdiDLQqrLrakMelmkPLPgIYnPDFIAA
           s khf vs    sdmd ye  g sk  ey   kq kpy ssvv   s  y      
             150       160       170       180       190       200 

            150       160       170       180       190       200  
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHR
                                                         ----TA+  R
846    NQdeRANncIngdeKqnVTTRGKWTHLQRIrrDIrnFKErnaldKViVLW----TANTER
         gs   qv klpk gq             id  qe   eekve  v             
             210       220       230       240       250           

            210       220        230       240       250       260 
GI||SP AGLTSRDTPSPVDPDTVEVLM-TFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIM
         -----  sp   DT+E L+-+ e D a++a Ssi             +          
846    Ys-----dlspGlNDTMENLLaaierDeaEIaPSsifAaAcIlEGVPfiNGSPQNTFVPG
        v     evvv v        qsvkn he  s  tly l s m    yv           
            260       270       280       290       300       310  

             270       280       290       300       310       320 
GI||SP KKARKIQVPEEQKDEKYWSRRYKNNEAAKRSRDARRLKENQISVRAAFLEKENALLRQEV
            i+       +   S   K   +         +K   I         +   L    
846    LiqlAehrnslIaGDDlKSGQTKlKSVLadFLVdAGIKPtSIaSYNHLGNNDGmNLSAPq
        vem ikegtf g   f      m    vq   g     v  v          y     k
            320       330       340       350       360       370  

712 aspartate-semialdehyde dehydrogenase asa ec 1.2.1. (376 aa)
initn: opt:  498 z-score: 79.3 E():    9.7
Smith-Waterman score: 500;    12.3% identity in 155 aa overlap
 25.8% noncontradicting positions,  13.5% class identity

       100       110       120       130       140       150       
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXERTLPFGDVEYVDLDAFXXXXXXXXXXXXXXX
                                     tL   D   + LD                 
712    VTaQGGDYTneifqaLRdaGWnGiWIDAASsLRMaDDAIIiLDPVNrnVIdaaLnrGirn
         c      ksvypk  es  q y      t   k     v     qd  tng kn vkt
                s       g                             p   d        
             80        90       100       110       120       130  

       160       170       180       190       200                 
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTS----------
                                             G-A   R  LT ----------
712    fiGGNCTVScMldglaGLFandladWmepapfqaagdgG-ArHMRELLTQmGhLnqaVad
       yv       l pmsyg   ntqvve tssmtysiqsgp   q        f l ygh kp
                               s v vf                              
            140       150       160       170        180       190 

        210       220               230       240       250        
GI||SP -RDTPSPVDPDTVEVLMTFEPDP--------ADLALSSIPGHETFDPRRHRFSEEELKPQ
       - dtPs    D      t   dp--------+ La s IP  + --      S EE K q
712    eLadPaaAILDIdRgVlarqrdgeLPqdnFGVPLaGnLIPWIDK--dLdnGqSREEWKaq
       l dt ss     e k ttltqsps  veh     g s         q gd m      ge
             200       210       220       230         240         

      260       270       280       290       300       310        
GI||SP PIMKKARKIQVPEEQKDEKYWSRRYKNNEAAKRSRDARRLKENQISVRAAFLEKENALLR
           K        =====       +       S          +s+             
712    AETNKILnrgaai=====PiDGLCVRiGAlRCHSQAlTIKLKrDVpideiEdLiAAgnqW
              gtssvf      v      v  m      f     k  slptv e l  htp 
     250       260            270       280       290       300    

614 flavonol udp-glucose 3-o-glucosyltransferase bronz (475 aa)
initn: opt:  482 z-score: 79.2 E():    9.7
Smith-Waterman score: 518;    10.8% identity in 176 aa overlap
 21.0% noncontradicting positions,  10.2% class identity

           40        50        60        70        80        90    
GI||SP XXXXXXXXXXXXXXXGLRSLLQGTSKPKEPASCLLKEKERKXXXXXXXXXXXXXXXXXXX
                         s           ASC L    R                    
614    VsCVVGDAFVWPaADAAaaAGAPWVaVWTAASCALLAHiRTDaLRrDVGDQAAnRaDelL
        t          t    ss      p            l   s  e       s v gp 
         130       140       150       160       170       180     

          100       110       120          130       140       150 
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXERT---LPFGDVEYVDLDAFXXXXXXXXX
                                       R ---LP     -V L+ F         
614    iaHaGLagYRVRDLPDGVVSGDFNYVInLLVHRqaQrLPraAaA-VALNTFPGLDPPDli
       vs p  gs                   s     mg c  ks t               vt
         190       200       210       220        230       240    

             160       170       180       190       200       210 
GI||SP XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXGTASGHRAGLTSRDTP
                           =        =                A g--+   S  T 
614    AALAaeLPNClPlGPYHLLl=AEddADTa=APADPHGCLAWLdRrPARg--VAYVSFGTn
           ei    v f      p   pt   n             g q   s          v
          250       260        270        280       290         300

             220       230       240       250       260       270 
GI||SP SPVDPDTVEVLMTFEPDPADLALSSIPGHETFDPRRHRFSEEELKPQPIMKKARKIQVPE
       +   PD    L +   d a   L S+         r            ++          
614    AcaRPDELrELAAGLEaSaAPFLWSLREdswaaaPrGFLdRAaGTGSGLVVPWAPQVaVL
        tp     q       d g         gvvphl p   e  p              g  
                                      tl                           
              310       320       330       340       350       360

Library scan:  0:02:01  total CPU time:  0:02:22

.
BCM HGSC