genewise $Name: wise2-2-0 $ (unreleased release)
This program is freely distributed under a GPL. See source directory
Copyright (c) GRL limited: portions of the code are from separate copyright

Query protein:       pstk_Drosophila_melanogaster
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            global
Target Sequence      subseq(404300,3000)
Strand:              forward
Start/End (protein)  global
Gene Paras:          human.gf
Codon Table:         codon.table
Subs error:          1e-05
Indel error:         1e-05
Model splice?        model
Model codon bias?    flat
Model intron bias?   tied
Null model           syn
Algorithm            623

genewise output
Score -83.47 bits over entire alignment
Scores as bits over a synchronous coding model

Warning: The bits scores is not probablistically correct for single seqs
See WWW help for more info

pstk_Drosophila    1 MRRICLVALIGLPAAGKSSLCSWLLGQQAALRVRHIVHLCYDDFLDATP 
                     M+  C     G P +GK +  + LL ++  L +   ++L +++++     
                     MK--CAFLFYGKPCSGKDTFINLLLTKRRQLEL--FLYL-FNNLIHKNE 
subseq(404300,3 1867 aa  tgtcttgactagagataacccaaaactgt  tctc taatacaag 
                     ta  gctttagacgggaacttatttcaggatat  ttat taattaaaa 
                     ga  cccgtttgcctggcccatacaggggaggg  tccc tttgttatg 


pstk_Drosophila   50 SADLAYKEQRGRIFKVIEKLISAIQEDTDWPPQVRRISSSGDYNSGRHL 
                     S  +  +E+   +F  + K    I++        RR +S    + G    
                     S-NVRLREK--IFFIQLIKYYYRIEKR-------RRKTSIFCGHGGETD 
subseq(404300,3 1999 t agctaga  attaccaatttcagac       aaaatattgcgggag 
                     c atgtgaa  ttttattaaaagtaag       ggaccttggaggaca 
                     c cgcagga  attcagtatctgtagg       gggtaattttgggtc 


pstk_Drosophila   99 ILCDDNFYYRSMRYKLYQLCRDSGCIFGQIYMASSLDSCLQANSLRSDA 
                           F  +   + L+ + R SG   G    +  + +      LR    
                     QKVSFRFLVQLTSHLLH-IGRSSGWDRGA--GSKGVKNGKTKEQLRRGL 
subseq(404300,3 2116 cagttcttgccaactcc agaatgtgcgg  gtaggaagaaagccccgc 
                     aatctgtttatcgatta tgggcggaggc  gcagtaagacaaatgggt 
                     gggtcctgcaatttgcc tggctagcgtg  atagtatagagaacatat 


pstk_Drosophila  148 TRVPVDVVRQMNERLEVPDTSEAWERNSLTLNGLDMDTT------GSAL 
                        P  +  Q+   L      +   R +LT    D+              
                     P--PRGITHQLYPTLRKVVQRKPLRRTTLTNQSDDFHLICYLYRHKDEW 
subseq(404300,3 2254 c  cagaaccctcaccaggcaaccaaaacaactggtccattttacaggt 
                     c  cggtcaatacctgattagactggcctcaacaatattgatagaaaag 
                     c  aaatctattacgcaacagatcgatgaccaactctatccgcacatag 


pstk_Drosophila  191 LAFIASLL---DLPAMETTLD-LPPAVARGLRQDQSLVHRLDLLLRTRI 
                     +    S L    +  +  +LD +   + RG   + + ++   ++ R  + 
                     MYHYHSYLRRNKIKTFNVSLDFIEKQLYRGHLPEVTSLNPSKIVQR-NL 
subseq(404300,3 2395 atctcttccaaaaaatagttgtagacttcgcccggatcactaagca ac 
                     taaaacatggaatactatctattaaataggatcatcctaccattag at 
                     gctttccatgtaaattcatatcaaggattacacaatcatggattag tt 


pstk_Drosophila  236 GALLDGLRQ--GEDKRLAGQTLNSRRKDILTKFRTEIGGKQIDGDVEGD 
                       LL+      G +K   G      +K  +      +         +G  
                     HLLLNRTYTVYGREKHYRGVYYPIVKKTTIKGGHVMLIRALNHKKEKGP 
subseq(404300,3 2539 cctcaaatagtgagactcggttcagaaaaaaggcgacaagcacaagagc 
                     atttagcactaggaaaaggtaacttaacctaggattttgctaaaaaagc 
                     tgagcagccttgggattgtccctcaaacccagtctggtgtcccaaaaag 


pstk_Drosophila  283 ALDYYVNLLN                                        
                      +   VN+                                          
                     NMTEAVNIAM                                        
subseq(404300,3 2686 aaagggaaga                                        
                     atcactatct                                        
                     cgaaagcaag                                        


//
Gene 1
Gene 1867 2715 
  Exon 1867 2715 phase 0
//
>subseq(404300,3000).[1867:2715].sp.tr
MKCAFLFYGKPCSGKDTFINLLLTKRRQLELFLYLFNNLIHKNESNVRLREKIFFIQLIK
YYYRIEKRRRKTSIFCGHGGETDQKVSFRFLVQLTSHLLHIGRSSGWDRGAGSKGVKNGK
TKEQLRRGLPPRGITHQLYPTLRKVVQRKPLRRTTLTNQSDDFHLICYLYRHKDEWMYHY
HSYLRRNKIKTFNVSLDFIEKQLYRGHLPEVTSLNPSKIVQRNLHLLLNRTYTVYGREKH
YRGVYYPIVKKTTIKGGHVMLIRALNHKKEKGPNMTEAVNIAM
//
>subseq(404300,3000).[1867:2715].sp
ATGAAATGCGCCTTCCTGTTTTATGGTAAGCCCTGCAGTGGGAAGGACACCTTCATAAAT
CTACTCCTAACGAAGAGGAGGCAATTGGAGTTGTTTCTCTACCTCTTTAATAATTTGATT
CATAAAAATGAGTCCAACGTGCGCTTAAGGGAGAAAATATTTTTTATCCAACTGATTAAA
TATTACTATCGGATTGAAAAGCGGAGGAGGAAGACTTCAATATTTTGTGGTCATGGGGGG
GAGACTGACCAGAAGGTGTCTTTCCGCTTTTTGGTCCAACTAACTAGTCATTTGCTCCAC
ATTGGGAGGAGCTCTGGATGGGACCGGGGTGCGGGATCTAAAGGGGTTAAAAATGGAAAG
ACAAAGGAACAACTCCGACGTGGACTTCCCCCAAGAGGAATTACCCATCAACTTTATCCA
ACCCTGCGCAAAGTAGTCCAAAGGAAACCTCTCAGGAGAACTACGCTAACCAACCAATCA
GACGATTTCCATCTAATTTGCTACTTGTACAGACACAAAGATGAATGGATGTACCATTAT
CATTCCTACCTACGTAGGAATAAAATAAAAACTTTTAACGTATCTTTAGATTTCATAGAA
AAGCAGTTATATCGTGGACACCTACCCGAAGTAACTTCCCTAAATCCGTCGAAAATTGTT
CAAAGGAATCTTCATCTGTTACTGAACAGAACGTACACCGTTTATGGGAGGGAGAAACAT
TATCGGGGTGTCTACTACCCTATCGTAAAAAAAACCACCATCAAAGGGGGTCACGTTATG
CTGATTAGGGCTCTCAACCACAAAAAAGAAAAAGGACCGAACATGACAGAAGCAGTGAAC
ATAGCAATG
//