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:       EMBOSS_001
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            global
Target Sequence      subseq(1321389,4689)
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 267.67 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

EMBOSS_001         1 MKAYYVVIILKVVLLLFFSYDFIFLKNEDNIIKGKGSSFLKSSYELLFV
                      K YY+ ++LK++LLLF+SYDFI LKN    +     +FLK SYE LFV
                     AKMYYLKLFLKLLLLLFLSYDFIVLKNGTREVVHGDRNFLKESYEFLFV
subseq(1321389,  997 gaattcattcatttccttatgtagtaagaagggcggaattagatgtttg
                     cataatatttatttttttgaattttaagcgattagagattaagaatttt
                     cagtcaaacgaaggcctacctcctggcgtgatgtacattaaattatacc


EMBOSS_001        50 QKDVLLPHIKKNNQISIFLCRS-QSQYIINKIKNYFDILNKGRETEIYF
                     QKD    +I ++NQISI+ C+S QS Y++++IK YFD+LN G +TEIYF
                     QKDTFPSYINRDNQISIYFCKSXQSHYVLSRIKKYFDLLNHGEDTEIYF
subseq(1321389, 1144 cagatcttaaagacatatttattctctgtacaaattgttacgggagatt
                     aaactccatagaaatctatgacgacaattggtaaatattaagaacatat
                     agcataacatgctacctttcacaatctgataagatctaactagtgaatt


EMBOSS_001        98 ENNNYIVDN----DQKYILFLMSIQPPLFFTLICYFILFIFILIALNIY
                     E ++Y  DN     +KY++FL +I+  L++ LIC F+LFI++LIA ++
                     EEDSYQADNSSNGEEKYVIFLRNIKAALIYRLICRFVLFIYVLIAASLI
subseq(1321389, 1291 gggatcggaaaagggatgattcaaaggtatccatctgctatgcaggata
                     aaagaacaaggagaaaattttgataccttagttggtttttatttccgtt
                     aactcaacctcttaagttatggccggaatcatacttctttccgtaatga


EMBOSS_001       143 LPKHCEFFIPSFLLNNTKFIETFEKMRKKKIIVSSIMFLAYNIIYSILC
                      P + + +IP +LLNN +F E FE M++KK++V  IMF +YN+IY +LC
                     FPHYLKPYIPPILLNNRRFNEVFENMQQKKLMVLGIMFFSYNVIYGMLC
subseq(1321389, 1438 tccttactaccaccaacataggtgaaccaacagcgaattttagatgatt
                     tcaatacatcctttaaggtaattaataaaattttgttttcaattagttg
                     cgttgagcaaatgactggctagtatgaaaaagatccgtttcccatgggc


EMBOSS_001       192 NTNEIHVYQNKNLIYKDYFNEYFFIQTLRTNLKNIHK
                     NT+ IHVYQN+NLIY DY+++++F+Q LR N+ ++
                     NTSVIHVYQNRNLIYDDYYSDHLFVQRLRENMTHVPG
subseq(1321389, 1585 aaagacgtcaaatatggtttgcctgccccgaaacgcg
                     acgttataaagattaaaaacaatttagtgaatcatcg
                     catgatgcatgtgtccttccttatggcccgtgttgtg


//
Gene 1
Gene 997 1695
  Exon 997 1695 phase 0
//
>subseq(1321389,4689).[997:1695].sp.tr
AKMYYLKLFLKLLLLLFLSYDFIVLKNGTREVVHGDRNFLKESYEFLFVQKDTFPSYINR
DNQISIYFCKS*QSHYVLSRIKKYFDLLNHGEDTEIYFEEDSYQADNSSNGEEKYVIFLR
NIKAALIYRLICRFVLFIYVLIAASLIFPHYLKPYIPPILLNNRRFNEVFENMQQKKLMV
LGIMFFSYNVIYGMLCNTSVIHVYQNRNLIYDDYYSDHLFVQRLRENMTHVPG
//
>subseq(1321389,4689).[997:1695].sp
GCCAAAATGTATTACCTAAAATTATTCCTGAAATTATTGTTGCTCCTCTTTTTAAGCTAC
GATTTCATCGTTTTGAAGAACGGGACTAGGGAAGTTGTGCATGGAGACAGAAATTTTTTA
AAAGAAAGTTATGAATTTTTATTCGTCCAAAAGGACACATTTCCATCATACATAAATAGG
GACAATCAAATCTCCATTTATTTTTGCAAATCCTGACAATCTCACTATGTGTTAAGTCGA
ATAAAGAAATATTTCGATTTATTAAACCATGGAGAGGATACGGAAATATATTTTGAAGAA
GACAGTTACCAAGCAGACAACAGTAGCAATGGTGAAGAAAAGTATGTTATATTTTTGCGG
AACATCAAGGCGGCATTAATTTACCGACTTATATGCCGTTTTGTCCTTTTTATTTACGTC
CTGATTGCAGCAAGTTTGATATTCCCGCATTATTTGAAACCGTACATACCACCAATTCTG
CTAAACAATCGGAGGTTCAATGAAGTGTTTGAAAATATGCAACAAAAAAAACTAATGGTA
CTTGGCATCATGTTTTTTTCTTACAACGTCATATATGGGATGTTGTGCAACACAAGTGTG
ATACATGTGTACCAAAATAGGAATTTGATTTACGACGATTATTACTCCGATCATCTATTT
GTGCAGCGCCTCCGCGAGAATATGACTCATGTGCCTGGG
//