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:       SelTryp;
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            local
Target Sequence      subseq(119362,6434)
Strand:              forward
Start/End (protein)  local
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 111.51 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

SelTryp;           8 ISFSLALLIAAFLAANAGSATNADNEATDYSSHRYKYLEMLKEADLKQM 
                     I  SLALL+  F A  A SA        + +   Y YL ML E DL+Q+ 
                     IVVSLALLF-VFAACRANSAHEVQ-PTQEGNKDPYAYLGMLDEVDLQQI 
subseq(119362,6 2469 aggtcgcct gtggtcgatgcggc cacggaagctgtcgacggggccca 
                     tttctcttt ttccggcaccaata ccaagaaacacatgttaatataat 
                     ctgaggttc gccgctatggtggg ctggacgcgcccctgctggcagga 


SelTryp;          57 LHEK--GEA-FHHLRSKRELILAVLKLEKREEALR               
                     LHEK  G+      RSK EL+ AV K+E+RE+A                 
                     LHEKTHGQVQVSAFRSKEELVAAVRKIEEREDAEV               
subseq(119362,6 2610 tcgaacgcgcgagtcaaggcggggcaaggcggggg               
                     taaacagatatgctggaaattcctgataagaacat               
                     ataggtaggggcccccgaggcgggtgtggcgcggc               


SelTryp;          89                           ARVTDT-----VQHEVRVEYCSG 
                                               A+ T +     V H++ + YC+G 
                               K:S[agc]        AKKTLSGRPLSVVHKLEILYCTG 
subseq(119362,6 2715 AGTTTCGA  Intron 1   CAGGCgaaatagccctggcacgacttag 
                      <1-----[2716 : 2967]-1>  caactgggctcttaatattagcg 
                                               gggggtgtctgcgtggatatcca 


SelTryp;         108 UGYRRHYEEVAESLLRSLPPELREQQKGKKPFIKFVGVVYSVGAFREFI 
                      GY +++E++ + L R LP             +++VG  Y     R +  
                     XGYAKYFEDMKQQLQRVLP---------NAGDVQIVGGTYPTPPARALA 
subseq(119362,6 3039 tgtgattggaacccccgcc         aggggcagggatcaccgcgcg 
                     ggacaataataaatagttc         acgatattggcacccccgctc 
                     aacggttatggggaacggc         cgttcatgttctccacccgcc 


SelTryp;         157 GNILSTGFLASIAISFFAPFLRGALPPHIAEWIEQHRGMVVGAGFMMNM 
                       + ST FLAS+ +++    L   LPP    +I Q RGM+VG GF++NM 
                     AKVCSTAFLASLGMALAGGQL-AFLPPAALNFIVQQRGMLVGTGFLLNM 
subseq(119362,6 3159 gagttagtcgtcgagcgggcc gtcccggcatagcccgacggagttcaa 
                     catgcccttcctgtctcggat cttcccctatttaaggtttgcgtttat 
                     gggcgactgccctgctcaagc gtgtcctgccttggtcggcctccggtg 


SelTryp;         206 VASSLLQSGAFEVYLNGSLIYSKLETGAVPTAETLADHILRQII      
                     + SSL Q+GAFEV L+G LI+SKL+ GAVPT E +   IL++ +      
                     IGSSLTQTGAFEVTLDGELIFSKLQAGAVPTVEVMRRIILQKTL      
subseq(119362,6 3303 agatcacaggtggatgggcattaccggggcagggaccaaccaac      
                     tggctcacgctatctagatttcatacgctcctattggtttaact      
                     ccctccgctccgtcgtcgtccggagtatagcggagctttgggat      


//
Gene 1
Gene 2469 3434 
  Exon 2469 2715 phase 0
  Exon 2968 3434 phase 1
//
>subseq(119362,6434).[2469:3434].sp
ATCGTTGTGTCACTGGCGCTTCTTTTCGTGTTCGCCGCGTGCCGTGCAAATTCGGCGCAT
GAGGTGCAGCCCACTCAGGAGGGAAACAAGGACCCGTACGCCTACCTCGGTATGCTCGAT
GAGGTGGACCTACAGCAGATATTACATGAAAAGACGCATGGACAGGTGCAGGTGAGCGCC
TTCCGCAGCAAGGAAGAGCTGGTCGCGGCGGTGCGTAAGATTGAGGAGCGCGAGGACGCG
GAGGTCAGCGCGAAGAAGACGTTGAGTGGGCGTCCCCTTTCGGTCGTGCATAAGCTGGAA
ATTCTATATTGCACCGGATGAGGATACGCGAAGTATTTTGAAGATATGAAGCAGCAGCTA
CAACGCGTGCTGCCCAACGCGGGTGATGTCCAAATTGTGGGTGGTACCTATCCCACCCCA
CCCGCCCGCGCGCTCGCCGCGAAGGTGTGCTCGACAGCCTTTCTGGCCTCCCTCGGTATG
GCCCTTGCCGGAGGACAGCTCGCGTTTCTGCCTCCCGCCGCTCTGAACTTCATTGTTCAG
CAGCGTGGCATGCTGGTCGGCACTGGCTTCTTGCTGAATATGATCGGCAGCTCTCTCACC
CAGACCGGTGCCTTCGAGGTTACCTTGGATGGCGAGCTTATCTTCTCGAAGCTACAGGCT
GGAGCTGTACCGACCGTGGAGGTAATGCGCCGTATTATTCTGCAGAAGACACTT
//