From 518ca89ce4ee69721810e1824ba711ce5d394281 Mon Sep 17 00:00:00 2001 From: Mattes Mrzik Date: Wed, 28 Oct 2026 10:59:05 +0100 Subject: [PATCH] working helper params with assert_eq(opti, clean) --- phyne/data/tkf/params/masa.fasta | 254 ++++++++++++++++++ phyne/data/tkf/params/tree.nwk | 1 + phyne/src/likelihood/mod.rs | 26 ++ phyne/src/optimisers/model_optimiser.rs | 71 ++++- phyne/src/optimisers/model_optimiser_tests.rs | 20 +- phyne/src/tkf_model/mod.rs | 16 ++ phyne/src/tkf_model/tkf92.rs | 90 ++++++- phyne/src/tkf_model/tkf_indel.rs | 40 +++ 8 files changed, 504 insertions(+), 14 deletions(-) create mode 100644 phyne/data/tkf/params/masa.fasta create mode 100644 phyne/data/tkf/params/tree.nwk diff --git a/phyne/data/tkf/params/masa.fasta b/phyne/data/tkf/params/masa.fasta new file mode 100644 index 00000000..8502c856 --- /dev/null +++ b/phyne/data/tkf/params/masa.fasta @@ -0,0 +1,254 @@ +>T19 +TAACCATA-CCTGCGCTGACAGTTTGACTGAT---CTCGGGAAGAGTTC-TAGGAAATATGCCACTTTGCGCCGAAGG------------------------------------------------------------------TGCATCTGCCTACACCGTGATGATTG----GT--------------------------GTCGGTAACTTTCACGCTCCGAATGGTCTTCACGCAATA-----------C-GGTCTGATGATATCCATGGCTCGTGG--TA-------GATG-CGATTGTCTACA-TAACACGTGTACGGCACT---GCGCACACAGGTTCTTCTGCGTGGGAGGTCCGACGCAAGAAGTTCAGAAGCCGGCATGGCAACTTGGTTTTCGCCTATGTTA---TGCCCGCACACTGAAACCCACTCTCTGCCTATCGAATCATT--------------TGTGTCGTAAGCTGCACAGTCTCTTTCACGAACGGCATACGGTCAAATGGTATTAAGCAACAATATCCGACTCCTTTCTGCA--------------G-----T-----CGACGA---G---------ACTAGCTGAA-GCAGAA-GTACAC-GTGTTGC-AGCCCTCAACTACG-----TATACCTATAGTCATCTTGGAC--------TTCCCGTAACGTCTCCGAGCGCCGTTGTG--------CCGGCATTTGAGTTATTTCGCC--------T----GAAGTAG------CC-TT-CAAGGCTGAGTCTTGTGTCACTGGCGTTG +>T38 +GTACCTTA-CCTCCGATGACAGTTTGACTGTG---CTCCGTACGAGCTC-TAGGGAAGTTAGCAC--------GTAGG------------------------------------------------------------------TGCATCTTCCAACACCCTGATGATAG----GA--------------------------TTCCCTTACTTTCTCTCTCCGAATGGGCTTCTCGCAATA-----------C-AGTATGATGATAGTCATGGCCTGTTG--TA-------GTTC-CGATTGACTACGCTTACACATTTACGTCACT---GGGCACGCAGTTACATCTGCGTTGCAGGTCCGACGCAAGAAACGCAGAGACCGGGATGATAACGTGATTCTCGTCTAGGTTA---TCCGCACAGTCTGAAACCAAATCTCCACCTATCAATTCAGT--------------TGTGTCGTACGCTGCGCAGTCTCTTGCACCCATGGCCTATGATCAGATGGTATTAAGCAAAAATACCCGACGCCTTTATGCA--------------G-----G-----GGACGA---C---------CCTAACTGAACGAAGAA-GTTCAG-GTAATGA-AGCACGCCCTGTCG-----TGTAGTAATAGACCTATTGGAC--------CTCCTCGAGCGTTTCCAACCGCCGTTGTA--------ACTGCCTTTGAGTTATTTTTCC--------T----GAAGTCG------CC-TT-CAACTCAAAGTCTTTTGTCACGGGCCTTG +>T9 +GCACTCTT-GATGCAACCACTGGTAGAGTTAT---CTTGGCAAGCCTAG-TTAGAAAAGTTTCAC--------GCGAA------------------------------------------------------------------TGTATCTTCCAAAACCATGAAAAATA----GGTAT-CTC-------------TGACGCATCACAGCCTCTCAAGATCCGAGTGGTCTCCTCGCAACA-----------C-AGTCTTCGGCTCGTCATGA-TTCTGG--CA-------TATC-GGACTTTCAATAAGTACGAGGATACTCCAGA---TGACACCCAGTTGCCGCTGCGTTTAAGTCA--AGACTCATATCGCGACGGCGCGCATGATAGTATGATACTCTTCGTGGCAA---ACCCAGCAGGGAGGTACCCTCTTTGTACTTAGCGTATGAGT--------------ATTCC-----------------------------------------------------------------------------------------------------------------------------------------CA-GTTTTG-GTA-CTA-GGACGGCCGGTACTCCTCTTGGCGCAAGAGACTCATTCGAT--------TTGGGCTGTCGTAGCCAAGCGCTAGTGGA--------CTGTCCTTTTCGCCATACCGCC--------T----GCCTTGC------CC-GG-AAGCTCCGCGTTTAACGTGCACGTCTTAG +>T64 +ACATCTTT-GATGCAATGACTGTATGATTGAT---CTTGGTAACCGTTG-GTGGAATACCGTCAT--------GCGCA------------------------------------------------------------------TCCATCGTCCAGCACCTTGACAAGAA----GATAC-ATG-------------AGAGGAATCACGGCCTCTCTAGAACCGACAGGTCTCCTCACAAAA-----------C-AGGATACGTATCATCGTGA-TTCTGA--GA-------TATC-CGGCCTCCTACCATTACACGGATATTCC------------CCAGTTGCGCCTGCGTGTACGGAA--AAACTCATTACGCAAAAACGGGCATGATGATTTTATCCCCGAGTAGTGAA---TCTGAGCCGACAGAGACCCACTCTGTACTTAGCGTAGCAGT--------------ATT-C--------------------------------------------------------------------------------------------------C----TGACGA---C---------GTGAATTGAACGACGCA-TTTTAG-GTAAGGA-GGACGGACAATCCGAATGTTAGTGAAATCGAACGCTTGACA--------CATGTCTGGCGTTGCCAAGAGCTAGTATA--------CCTT----TGCGCTGTTCCGCC--------T----GAATTGC------CC-GT--AAATCCCAGCTTTCTTTCACTGTCCTAC +>T26 +GCATGGTT-AAAGGAATGGCAGTTTAACTGTT---CTTGGTAATCGATC-TTGGAAATACGCCGC--------GTGAA------------------------------------------------------------------TTCATCTTCCAACATTTTGATAT--G----GATAC-CTGAAATGCTAGGTTGGGACGGATCAAAGTCACTCTAGATCCAGCGGGTCTCATCGCTACA-----------A-AGACTACAAGTATTTGTGA-TTCGGA--CC-------TATC-CGAATCGAAACCATGCCACAGATGTTCCAAC---GGACCGGCAGGCGATTGTGGGAGGTAGCCA--GGACTCATAACCTAAAACAGTGAAGGTTAAAATGGAACCCTTTTTAGGAA---AGCGCGCATCCCGATACACTCACTGTACTTAGAGGATTATA--------------ATTGT-------------------------------------------------------------------------------------------G-----T-----GGTCAC---C---------ATGAAGTGCATGCGTCC-GTTTAA-GTAACCA-AGACGGGCGTCACGAATTGTGCATGAATGGACCTCGTGCCT--------GAGCTCT-------CCGGTCGGTCCTGTA--------CCTATCTTTGATCTGGGTCCAC--------A----CAAGTCCAGCCAGCG-TT-GATCTCGGAGTTATTGATCACGGTCGTTT +>T45 +GCAACTTT-GATGCAATGACTGTTTGATTGAT---CTTGGTAAACGTTC-TTGGAATAATGCCAC--------GTGAA------------------------------------------------------------------TGCATCTTCCAACATCTTGATAAGTG----GATAG-CTG-------------AGACGAATCACGGACTCTCTAGAGCCGACAGGTCCCCTCGCAACA-----------C-AGACTTCTGGTAATCATGA-TTCGGA--CA-------TATA-CGAATTTAAACTATTCCACAGATGTTCCAAC---GGACACGCAGTTGCTGCTGCGTGGTGATAA--AGACTCATAACACAGAACCGTGCAGGATAAAATGATACCGTTTTAGGCAA---TGCGCGCAGACAGATACACTCTCTGTACATAGCGTATCAGT--------------ATTGC-------------------------------------------------------------------------------------------G-----T-----GGACAA---C---------ATGAATTGCACGAGGCA-GTTTAG-GTAACCA-AGACGGCCGGTGCGAATGTTGTAGGAATAGACCTCTTGGCT--------ATGCTCTGGCGTTTCCGAGCGCTAGTGTC--------CCTATCTTTGATCTGTTTCGAC--------T----GAAGTCC------CG-TT-GAAGTCCGAGTTTTATATCACTGTCTTAG +>T53 +GCCTCCTA-GACGCACTAACAGGTTTATTTAT---GATCGTAAGAGTTC-TTGGAAGGATACCAG--------GGGAA------------------------------------------------------------------TGCAAGTTCCTACACCAAGAATAGTA----GATAT-CTG-------------AGCCAACTCACTGTCCTTCCAGATTCGGCACGTTTCCTGGCAACA-----------C-A-GTATCGGATAAAGAAAA-CTCTGA--CA-------ATTT-CGAATTTAAGCCATAGCAGAGAGATGCCGAC---TGACGCACAATTGCTAATGCGTGGTAGTAA--TAACACAAGACGCAGAAACGTGCACGATCATCTGATGCCCATCTAGGCAG---TGCCAGCAGTCCGACACCCTTACTGTACTTAGCGTGTCTGT--------------CTTGC-------------------------------------------------------------------------------------------G-----T-----GCACTA---T---------ATTAACGGCAGGCGGCA-AATTAG-GTGACCA-AGACGGCCTGTACGAATGATTGAC-AAAAGACATCTTGGCG--------ATGATGTCGCGTTTCCATGCGCTTCTATA--------GCTGCTTATGAATCTTTCC---------------------------------------------------ACCAATCGGGAAG +>T58 +CG-CCTTC-ACTGGGGTGAAAGGAACGGGAGT---CTTAGTATGACGTC-TCGCTTAGACGGGGA--------GCACG------------------------------------------------------------------TGAATTCTCCACCTCCCTCAACATTA----TACCA-GAG-------------ACCTGA-TAATCAAGGGTCAAATGCAAACAGGGTTATTTAGCGGT-----------G-CGTACCTTT-TATCCATAGGACCTTTGACC-------ATTC-TGGCAGACTAAAATACTAGTCATACCAAAAT---GGACACTTATTCGGTACCTCTTGGTAGGGGGGATACCAACGAGCTAGAGACCGGTATGGTGGTTCAATCTCCGTCTAACCAT---CGTGTATACCCTGCACCCCGTGCGCTGCTCCTCGATTGATG--------------TATGT-------------------------------------------------------------------------------------------------------GAGA-T---T---------ACGTGCTGAATCGCCTA-TTTATG-GAATCCA-GATACGCCCTTAGAACTGGTACAACTCGAGTTCCAATAGAA--------CTAAACTTTCTTTTCCGGGGACCCTAATA--------CCCAACGCTGTGCTGATTTTAC--------T----TAGGTAC------CT-TG-CAACCCGGTGAAGATTCTTAATCTGTCAG +>T27 +GG-CCTTC-ACTGGGCTGAAAGGAACGCGAGT---CTTAGTGAGAGGTC-TCGCTTAGACGGGAA--------GCCCG------------------------------------------------------------------TGAATTCCCCATCTACGTCAACATTA----AACCA-GAG-------------CCATGA-TAATCAAGCATCAAATGCAGACAGGATTT-ATAGCGGT-----------G-CGTACGTTT-TATCCATAGGACCTTAGACC-------CTTC-TGGCAAACTAAAGGACTAGTCATACCAAAAT---GTAGGCTTATTCGGTACCTCGTGCTAGGGGGGATACAAACGACCTAGAGACCGGTATGGTGGTATAATCTCCGTCTACCCAT---CGGGAATACCCTGCTCCCCGTGCGGTGCTCCTCTATTGATG--------------TATGT-------------------------------------------------------------------------------------------------T-----GAGAGT---T---------ACGTGCTGAACCGCCTA-TTTAAG-GTAACCA-GAAACGCCGTTAGAACTGGTACAACCCTAGTTATCATAGAA--------CTGACATATCCTTTCTGGGCGCTCTCATA--------CCCAACACTGTGGTGATTTTACTCAGCACGT----TAGGTAC------CT-TG-GAATCCGGAGAAGAATCTTAATCTGTCAG +>T20 +AG-CCTTC-CCGGGGCTGAAAGGATCGCGAGT---CTTAGTAAGAGGTC-TTGCTTAGACGGGAA--------GTACA------------------------------------------------------------------TGTATTCTGCATCTACCTCAACATAC----AACCA-GAG-------------CCATGA-TAATCAAGTATTAAATGCAGACAGGATTA-ATAGCGGT-----------G-CGTACGTTT-TATCTATAGGACCTTAGATC-------ATTC-TGGCCCACAAAAGGTCTAGTCCAACCACAAT---GTACACTTATTCGGTACCTCGGGCTCTGGGGGATGCAAACACCCTAGAGACAGCTATGGTGGTATAATCTCAGTCTAACCAT---CGTGTATGCCCTGCTCCCCGTTTGGTGTTCCTCTATTGATG--------------AGTGT-------------------------------------------------------------------------------------------------T-----GCGAGT---T---------ACGTGCGGAACTGCCTA-TTTATG-GTAACCA-GAAACGCCGTTAGAACTCGTACAAACCGAGTTCTCATAGAA--------CTGATATATCCTTTCTGGGCACCCTCATA--------CTCAACCGTGTGCTGATCTTAC--------T----TTGGTAC------CT-AG-GAACTCGGAGAAGAATCTTAATCTGTCAG +>T28 +CGAAATTA-GCTGCTACAACAGGAATACTGAT---CTTGGAAAGAGTTG-TAAGGAGAATGGTGA--------GCTGG------------------------------------------------------------------TGGACCTTGTAAGAACTAAATCCGTA----CATCG-GTG-------------AGAAGAGGATCTAACTCTCCACTGTTGACATCTCGTCACGCCGCA-----------C-ATGATGCTT-TACGCATTGCTTCCCG--CGACAGTACCTTCAGGACGCTCTCAAAATTTAAGCCTACAATATC---GCGCGCTCAGTGGCTTCAGCGTGGTAGACTTGATACATAGATCGCAAGAAGCGGGATAGAAA------------------AC---GGAGCAGGGAGAGGAACGATGTCGGTCCGTATCAATTAAGT--------------GATGT-------------------------------------------------------------------------------------------T-----A-----GAACGA---G---------ACAAACTGCACGAGGTA-GATTCG-GTATCGC-AGGCCAACGGTAATAATTTAGGACCATTAGCCCTATTAGAT--------CTGCCATAGTGTTTCCACGT-CCACTGTA--------GACGGCTTTGTACGATTACGAC--------C----GAAGTAC------CC-GT-GAACTCAGTCATTAAAGTCATTCTGTAAT +>T55 +AGACCTTA-GCATTGTTGGCAGCACCGCTGGA---TATGGACGTAGA-C-TGAGAAAAGTGGTTA--------TCGAGTACTCCCAGGCGGTGTTGAGCCGCGAGAGCCTGAGAGTCCTTAAGTCTTGTAGTCCGAAAC-----TGGAGCATGCAAAAACCCCATGGGAA----AATTG-ATA-------------ACAAGAGTATTAGACTCCCTAGGGTCGGCAGGTCTTCATGCAGCA-----------C-AGGATGCAG-CAAGCTGGGCTTCTGG--AT-------TACC-ATACCGCATCACATTAAACCCTTAAGTCGAC---AAACACTTAGTGTAGACAGTGATTGATGAGGGATACCTAGCTCCCAGCGACCGGTATTGCAGATTGAGGCCCTTCCATGCAG---GGGACGTAGAGACAAACCCAAGCTTTACTTATGGATTGAGC--------------GATGG-------------------------------------------------------------------------------------------G-----T-GCGCGGGCCC---C-------T---TAACTG-ATCAGGCA-ATTGAG-CACACTA-CGCCCGAAGCTCCAACTG-------AGTAGACCTCTTAGAT--------CTAATCTATCCCTACTATGCGACGTCGGA--------CCAGCGTTCTAGTTGGTCC-GC--------T----TTAGTCC------CC-TA-TACCTAAGCCTTAAG-CTTAACGTCTTGG +>T60 +TGACCTTA-ACTTCGTTGGCAGTACCACTGGA---TATGGTCGTAGA-C-TGAGAAAAGTTGTCA--------TCGAG------------------------------------------------------------------TGGAACATGCAATAACCTCATGGTAA----AATTG-ATA-------------ACAAGAGTAATAGACTCCCTAGGGCCGGCAGGTCTTCTAGCAGCA-----------C-AGGATGCAG-TAAGCTGGGCTTCTGG--CT-------TAGC-GTAGCGTATCACATTAAACACTTACGACGAC---AAACGCATAGTCGAGACAGGGATGGAGGAGGGACACCTAACACGTAGCAACCGGTATTGCAAATTGAGGCCCTTCTATGCAG---AGCACGTAGAGACAAACCCAAGCGGTACTTATGGATCCAGC--------------GATGT-------------------------------------------------------------------------------------------G-----T-GCGCGGGCCG---C-------T---TAACTG-ATCAGCAA-ATTGAG-CACACGA-CGCCCGAAGCTCCAACGG-------AGTTGACCTCTTCGAT--------CTCATCTAGCCCTACTAAGCTCCGTCGGA--------CCCGCGTTCTAGCTGGTCC-GC--------T----TCAGTAC------CC-TA-GACCTAAGGCTTACA-CTCAATGTCTTAG +>T49 +TGGCCTTA-GATGCGTTCTCAGTACCACTGGG---CATGGTAGGCACTC-TGAAAAAAATCCCCA--------ACAAG------------------------------------------------------------------TGCAACATGCAACAACCACATGAGAA----AATTG-ATA-------------ACAAAAGTAACAGACTCCCTAGGGCCGTCAGGTCGTCTAGCAGCG-----------C-AGGAAGCAG-TAAGCATGCCTTCTGG--AA-------TAGC-GTACCGTCCAACATTAAACACTCACGACAAC---AGTCAAATAGTTGAGACGCGGTAGGACTAGAGACACCTGGAACGCAGCAACCGGTATTGCAAATTGTGGCCCTTCTATGAAT---AGCGCGGAGAGACAAACCCAAGCGGTACTTATCGATCCTGA--------------GATGA-------------------------------------------------------------------------------------------G-----T-GGGCGGTCCG---C-------T---AAACTGGATCACGAA-GTTTAG-GTAACGA-CGCCCGAATGCACAACGG-------AATGGACCTCTTAGAT--------CTGATCTAGCGCTTCCAAGCGCCGGGGAA--------CCCGCGTTCGAGCTGGTCT-AC--------A----TAAGTTC------AC-TT-GACCCAAGACTAATATCAAAATGCCTTGG +>T63 +TGGCCTTA-GGTACGTTCCCAGTACGACTGGG---CTGGGTAGCAAATC-TTAAAAAAATCACCA--------GCAAG------------------------------------------------------------------TGCAACATGCAACAACCTCATGAGAG----AATTG-ATA-------------AGAAAAGTAACAGACTCCCTAGGGCCGACAGGTCGTCTCGCCTCG-----------C-AGGAAGAAG-TAAGGGTGGCTTCTGG--AC-------AAGC-CTACCGTCCAACATTAAACACTTACTGCAAC---AGTTAAGTAGTTGAGACAGCGTAGGACTTGATTCACCTGGAACACAGCACCCGGTAATGTAAACTGTGTCCCTTCTATGAAT---AGCGCCGAGAAACCAAACCAAGCGGTACTTATCGATCCCTA--------------GTTGA-------------------------------------------------------------------------------------------C-----T-GGGCGCTCCG---C-------A---AAACTGGAACAGGAA-GTTTAG-GTCACAA-CGCCCAAATGTACAACGG-------AATTGAACTCGTAGAT--------CAGAGCTAGCGCTCCCAAGCGCCGGGGAA--------CCCGGGTTTCAACTGGTCT-CC--------A----TACGTAC------AC-TT-GGCCTTATTCTAATATCAACATGCCCTTT +>T11 +TGATACAC-AATGCGTTGGCAGGAT-TTATAT---CGT-------ATTCCTACGAACAATAGCCA--------GCTCG------------------------------------------------------------------CGGTTCTTTCAATCTCCTTATGCGTA----GATTG-CTC-------------GCAAGGGTAGCTGA---TCTAGCGCCTATAGGTCATAAGGCGCCA-----------G-AGAAAGCAG-TAAATATGGCTATTAG--CA-------ACTC-GGAGTAGGAAACCTTATTCACTTGCGACTCA---GGAGACTGCCTTGCTGCTGCGTGATCAGGGAGATATAAAGTGCCTACCAAC--GGCATGCAATTTTGGCACTTTTTATGCAA---GCGGCCCAGAGAGGAACGCCAGCCCT-CGTTACGGTCGCTT--------------CAGGT-------------------------------------------------------------------------------------------------T-----GGACAA---A-------CAACAAACTGGTTAAGGCA-GGCTGG-ATATGCA-AGACCGTAGGTACGAACT-------AATAGACCACTTGGTC--------GTGATCTGGCGCGTCCGGGCGCTATTGCA--------TCCAGCTTTGAGCTGATTCGAC--------T----TCAGTTC------GC-TG-GGTCTCACATTGTGAGCTCGATATGTCAG +>T16 +GGATACTC-GCTGCGCTGACAGGGT-TTATGG---CTT-------ATGTCTAGGAATCATTGCCA--------TCACG------------------------------------------------------------------GGTAACATAGTATATCATTATGCATA----GATTA-CTG-------------TAAAGCGTAGCTAT---TCTAGGCCCGCCAGGTACTATAGCACCT-----------C-ATAGTGGAG-TTAATATAACTTTGGG--GG--------CTC-TGACTAGCAAACACTTAACACTTGTGGCTGA---C-ACGCTTCATTGCAGCTGCGTGATCAGTGTGAGTTAAAGTATCTAGCATC--GTCAAGCAATTATGGGACATTCTATGCAA---GGCGCCCACAAAGTAAAGCTAGCCTC-CGTACCGATCGATT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGAGCA---C-------TTAAAAAGTGTTGAAGGAT-GCTTGC-ATATCGT-AGACCGTTGGGCCGTACC-------GATTGACCAATTGGTT--------ATGATCTATCGCGTAAGTGCCCTATTGCA--------CCCAGCTTTCAGTTGCTTCGTC--------C----TCAGTGC------GC-TA-GCACTCGCATTCTACGCCCAATGTCTCAC +>T17 +TAACCTGA-GGT-CATGTATGG-----------------------------------------------------------------------------------------------------------------------------------------GCTCACCCTGA-----AGTG-ATG-------------TCAAACGTAACTGCCTCTACAGGGAAGACAGCTATTGAGTTTCGC-----------A-GGTGCGGTC-ATAGGATGGATTGAGA--CA-------TCGA-ATGCTGTCTAACATTATTCTGTCGCGACCAC---AGGTAGTGAGTTGTATTGCGTTTAAGCTAAAGCCCCACAGTACCCAGCATCAGGTTTTCTTGCACAGGGCCA-TGTATCTTG---TGTTCGCTGTGAAAAAGACTACCTGTAGATATCTATCCACG--------------AATTT-------------------------------------------------------------------------------------------C-----T------GGACA---CATAA-GGTCCCAAAGTGC-CCAGCAA-GTTTGG-TAACCGC-AGGCCGGCGGACTATACG-------TGTAAAACTATAGGAT--------TTGTCTCACAGCGGCTTCCGACCGTTGTA--------CCCGAATTTGCCCTCGTTCTCT--------T----TCTTGAT------GC-TG-GAGCGAAGACGGAAAACTCGATTTCGTAC +>T7 +TAACCTTA-GGT-CTTTGATGGGATCACTGAT---CATGGTAAACGTTA-TGGGTCAAATGGCCA--------ACAAG------------------------------------------------------------------CGGAAGTAGCAGTCACCTCATCAGGA----GATTG-CTT-------------ACAACCGTAACCGACTCTACAGGGCCGACAGGTCTTCTATCTGGA-----------A-AGTGTGTAG-AGAGAATGGCTTTTGA--AT-------TAGC-GGCCTGTCTAACATTATACTCTGACAACAAC---AGCTACTTTGTTGTATCTGGGTGGGAATAGAGACTCGAAGTACGCAGCAGCAGGTATTGAGGCTCGGGTCCA-TCTGTGCAA---TGTTCGCTGAGAGGAAGACTAGCGGTACGTATCTATCCAGT--------------AATTT-------------------------------------------------------------------------------------------G-----T------GGCCA---T-------TCACATAACGA-CCAGGAA-GTTTGG-GAACCGC-AGCCCGGAGGTATAGACG-------CATAGACCTCTCGGAT--------CTGTCTTAGAGCTGCTAACCCCCAATGTA--------CCCGCATTTGCCCTCGTTCGCT--------T----TTTTCAT------CC-TT-GAGCTCAGACTTTTATCTCGATGTCTTAA +>T43 +TGACATTA-GGT-CCTTGATGGGATGACTGAT---CGTGGTAAACGTTA-TAGGACAAATGGCCC--------ACAAA------------------------------------------------------------------CGGAAGTAGCAGTCACTTCATCAGGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGGGCAGGCAGGTCTTCTATCTGCA-----------A-AGGGTGTTG-AGAGAATGGCTTCTGA--GT-------TAGC-GGCCTGTCTAACATTATACTCTGATGACAAC---AGCTACTTAGTTGTATCACGGTGGGAATAGCGACTCAAAGTACGGACCAGCAGGTATTGAGGTTCGGGTCCA-TCTTTGCAA---TGTTCGCTGAGAGAACGACTAGCGGTACGTATCTATCCAGT--------------AATTT-------------------------------------------------------------------------------------------G-----G------GGCCA---T-------TCAAAAAACGA-CCAGGAA-GTTTAG-GAACCGA-AAACCGGAGGTATAGACG-------CAGAGACCTCTCGGAT--------CTGTCTTAAAGCTCCCAACCCCCAATGTA--------CCCGCATTTGCGCTGGTTCGAT--------T----TTTGCAT------CC-TT-GAGCTCAGACTCTTATCTCGAAGTCTTAA +>T12 +TGACATTA-CGT-CCTTGATGGGATCACTGAT---CGTGGTAAACGTTA-TAGGACAAATGGCCA--------ACAAA------------------------------------------------------------------CGGAAGTAGCAGTCACTTAGTCAAGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGCGCCGGGAGATCTTCTATCTTCA-----------A-CTTGTATTG-TGACAATGGCTTCTGA--GT-------TAGC-CGCCTGTCTAACATTATACACTGATGACAAT---AGGAACTTATTTGTTTCACGGTGGGAATAGAGACTCAAAGTACGGAGCAGCAGGTATTGAGGTTCGGTTCCA-TCTGTTCAA---TATTGGGTGAGAGAACGGCTATCCGGACGCATGTATCCAGA--------------AATTT-------------------------------------------------------------------------------------------G-----G------GGACA---T-------TCAGAAAACTA-TCCGGAG-GTTGAG-GAAACGA-AATCCGGAGGTATAGACG-------CAGAGACCTCTCGGAT--------CTGTCTTAAAGCTCCAAAGCGCCAATGTA--------CCCGCATTAGCATTGGTTCGAT--------G----TTCGCAT------AC-TC-GAGCTCTGACTCTTATATCGAAGTCTTAA +>T61 +ATACAGTA-GGT-CATTTATGGATTAACTGAT---AGTGGTTAACTTTA-TAGGACAAATGACCA--------ACAAG------------------------------------------------------------------CAACCGTAGCAGTCACTGAGTCAAGA----------CTT-------------ACAAGCGTAACAGTCTCTACCGGGGCTGCAAGTTTTCTATCTTCA-----------A-CTAGTATTG-TGATGATGGCTTCTGA--GT-------TCGC-CGCCTGACTAGCATTATAGAGTGATGACTAA---AGTCACTTTGTTGTTTGTCGGTGACGATCTAAACTAAAAGGACCGGGCAGCAAGTATTGAGGTTCGGCTCCA-GGTT-AAAA---TATTTGGGGCGAGGACGTCTATTGGTGGACATGTCACCAGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGGG---T-------ACAGAAACCTA-TCCGCTT-GGTGCG-GATACTA-ATTCCGGAGGTATAGACG-------AAGAGAGCTTGCGGTT--------CTGTTTCTTAGCTCCAAAGTGGTTATGTA--------GCCGCATTAGCATTCGTTAGGT--------G----TACGGAT------AC-CT-GAGCTCGGACTCTAATAAGG-AGTATTAA +>T15 +ATACAGTA-GGT-CTTTTATGGAGAAATTGAT---AGAGGTCAACTTTA-TAGGACCGATCACCA--------ACAAG------------------------------------------------------------------CACCCGTATCAGCCACTGCGCCAGGA----------CTT-------------ACAAGCGATACAGTCTCGACAGGGGCTGCAACTCTTCTATCTTCC-----------A-CCAGTATTG-TGAAGAGGGCATCTGA--GT-------TCGC-TGCCTGTCAACCATGATAGAGGGATGCCTAA---AGTCACTGTGTTGTTTGTCGGTGACGATCTCTTCTAAAAGGACCGAGCAGCAAGGACTGAGGCTCGGTTCCA-GCTT-TTAA---TATTTGGATCGAGCACGCCTATCGCTTGACATCTCTGCAGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGCG---T--------------CCTA-TCCGTTG-GTTGAT-GTTGCTA-ATTCCAGAGGTATAGCCG-------ATGAGATCTAGCGGTT--------CTGGTTTTTAGCTCCAAAGAGGTTATGTA--------GCCGCATTCGCATTCGTTAGGT--------G----TACGCAT------AC-TT-GCACTCGGACTCTAATAACG-AGTGTTAG +>T2 +ATACGATT-TGC-CTTTTCTGAATTAACTGAT---AGAGTTCCAAGTTA-TAAGAGAGATGATCA--------ACAAG------------------------------------------------------------------CACCCGAAGCAGTCACTGTGTCACGA---------ACTT-------------ACAAGAGATACCGCCTCGACATCGGCTGCAACGCTTCTACCTTCC-----------A-CGAATATTG-TGCTGATGGTTTCTGA--AT-------TCCC-CGCCTGACTGTCATGATAGAGCGATGCCTAA---CGACACTTTGTTGATTGGCCATAACGAACTAAAGTAAAAGGACCGGGGAGCAAGGTTTGACGCCCGCTACCA-GTAA-TTAA---TATTTCGGTCGAACTCGCCCATCGGCGGAGATGTCTGCGGA--------------AA-----------------------------------------------------------------------------------------------------------------------------------CGG-TCGGATC-GTTGAGGGTTACTA-GTTCAGGAGGTATAGCAG-------AAGAGATCTTGCCGTT--------TTCGTTTTTAGCTGTAATGTTCTTATGTAGTCTGTATGCTGCATTAGGATTCTTTAGGT--------G----TAAGCGA------AC-CT-AACCTCGGAGTCTAATTATT-CGTGTTTG +>T62 +ATACAATT-TGA-CACTTATGAATTAACTTAT---AGAGGTCCACTTTA-TAGGACAGATGATCA--------ACAAG------------------------------------------------------------------CACCCGTAGCAGTCACTGAGTCAGGA----------CTT-------------ACAAGCGATACAGTCTCGACAGGGGCTGCAACGCTTCTATCATCA-----------A-CTAGTATTG-TGATGATGGCTTCTGA--GT-------TCGC-CGCCTGACTACCATGATAGAGTGATGCCTAA---AGACACTTTGTTGTTTGTCGGTGACGATCTAAACTAAAAGGACCGGGCAGCAAGGATTGAGGCTCGGTTCCA-GCAA-TTAA---TATTTGGGTCGAGCACGCCTATCGGCTGACATGTCTGCGGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGCG---C--------------CCTA-TCCGTTG-TTTGAG-GTTACTA-ATTCCGGAGGTATAGCCG-------AAGAGATCTAGCCGTT--------ATCGTTTTTAGCTCCAAAGTGCTTATGTA--------GCCGCATTGGGATTCGTTAGGT--------G----TACGCAT------AG-TT-GAACTCGGACTCTAATAACG-AGTGTTAG +>T59 +GGCCGTCA-GCGGCCATGACAGGATCACTGAT---CTTGGCAAGAGCTC-CGGGAAAAATGGCCG--------GGAAG------------------------------------------------------------------TGGATCTTGCAACAACCAGACCAGTA----GATCG-CTG-------------ACGAGAGTAACTGACTCTCTTGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TATGCATGGCTTCTAT--TA-------ACTC-GGCCTGCATAACTATACACACTTACGACAAC---GGACACTTAGTTGCTTCTGCGTGGTAGGAGAGGCACTAAGACCGCAGCCACCGGAATGGTAATAGGGATCCCCTCTATGCAA---GACGCGCAGAGAGTAGCCCTATCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGG---C---------ACAAACTGAACCAGGCA-GTTTAG-GTAACGA-AGAACGCCGGGACAACTGTTGGACCAATAGAGCTCTTGGAT--------CTCCACTAACGTATCCAAGCGCCATTGTA--------TCCGCTTCTGAGCTGTTTCGAC--------T----TAAGTAC------GA-TT-AAAGTCAGACTTTGATCTAAGGGTCCTAG +>T34 +ACTACTTT-TGGGACGCAACAGGATCACCGTA---GTTCTTAAGCGTTG-TGGCCGTACTATCCA--------GCCAT------------------------------------------------------------------GAAATCATGCAACTAACTGCCAATTA----GGAAG-CCG-------------ACGAGAGGAATGGGCCGGCTCGCCCCGGAAAGTACTCGTACAGCG-----------C-TAGCGTCAC--ATGGATGGCTTAAAT--CA-------AGTT-GTCATGCCAATCGAAGCAAAGTTACGAAAAT---CGCCATGTAGTATCTTCAGCCTGGTGCTAGGGGCGCTAAGACGGGAGCGTAGATGATAGTAACCCTGGCCCGGCCTATGCGT---GCCAGCGAGCGGGACATGCTATCGCGACGTACATCGTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCAAGA---C---------ACACAGTCAACCAGGCA-GTTTAA-GTTCCGA-CGCAAACGGATACATATGCTCGAAGAATAGAGCCGTCGTCG--------CTGAAGTAGAGTTTAGCATCGCTATCGTA--------CACGCTTATCTCTTGT-TGCAA--------T----AAAGTAC------CA-CT-GAAGTCACATAAGTACACCAGGAGCCCAT +>T35 +TCTACTAG-ATAGAACCAAGGGGGTCACTGTG---CTACTGAATCGATG-TGGCTCTCCTATCCT--------ACCTT-------------------------------------------------------------------CAATTAGCCAACAAACTGCATAGTC----CGATG-CCG-------------ACGAGAGGAGAGGCCTGGCTCTCCCGCGTTAGTGCTCGTACTGGG-----------A-TCGAATCTC---TGCATCGGTGATAT--CA-------GTGT-GTCATGCCAATTTCTGCAAACAGACTAAATA---GGTCATGTAGTGTCTTCAGCTTGATGGTAGGGGCTCAAAGTACGTAGGGTAGGAGATCGTAACTCCCGGACACCCTATGCAT---GGAAGCCAGTGGACCATACTAAAGGAACGTACTTAGTGAGT--------------TAATG-------------------------------------------------------------------------------------------------------TCCTGT---C---------GCCAAGCCAACCAGACA-CTCCAA-TTGACGA-CCCAAACGGATACATGGGGTCGAACTGTAGAACCTTTGTCG--------TTGTAGTATACGTTAGCACCGCTATCTAA--------CACGCTTACTACTGGA-TTCGA--------T----GAAGTAC------CT-CT-GAAGTCACATGAGTAGACCAGGAGCCCAT +>T18 +AGTACTGA-TTTGAAACAACAGGGTCACAGTG---CTTTTGAACCGTTA-TGGCTCTCCTATCCA--------ACCAG-------------------------------------------------------------------CAATTATCC-ACAAACTGCTCAGTG----CGGTG-CCG-------------ACGAGAGAGGAGGTCCGGGTCCCCCCATAGAGTACTCGTACAGCG-----------CCTCGACGAAC---TGCACCGGTTAAAT--CA-------AGTT-GTCATGCCATTCAGCGCGAAGATACGAACAT---GGCCATGTATTGTCTTCAGCCTGATGG-AGGGGCGCTATGAAGGGAGGGTAGAAGATCTTAACTCTCGGCCGGCCTATGCGT---GGAAGCCAGCCTAAGATTCTAAAAGAACGTACTTACTGAGT--------------GAATG-------------------------------------------------------------------------------------------------------ACTTGT---C---------ACCATGCCAGCCAGGCG-GTTAAA-TTTAATA-CGCAAAAAAATACCTATGGCCGAACTTTAGAGCCTTCGTCG--------TTATACAAGACCCTAGTACCGCTCTCGAA--------CACGCTTATCTCTGGT-TACCA--------T----AGAGTAC------CT-CT-GAAGTCGGATGAGGATACCAGCTGCCCAT +>T3 +AGCACTTT-GGGGACATTGCAGGACCACCTTA---CTACTTAATCGCGG-TCTGTGTACAACCCA--------GGAAT-------------------------------------------------------------------GGAGCGTGCAACTAACTGACCCTAA----GAAAA-CCG-------------ACGAGACGAACTGACCCTCTCAACCCGAAAAGAATGCGCACAGCG-----------C-GAGCTTCGC-CATGGACAGCTTCTAG--CC-------ACAC-ATCCCGCCAACCTACTGAAAATTACGAAAAT---CACAACATAGTATTTGTGGCCTGGTGGGAGTGGCACTAAGGCCGTAGCGGATGTAAGGGTAAGCCGTGGCGCGTCTATGTTA---AGTGTGCTGCGAGACATGATATCGCTTCGTATCCCCTCAAT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCCA---G---------ACAAAGTGAACTAGGGA-GATTGA-GACAGGT-AGGAAGCGGAAACATATGCTTGAACTATCCACCCGTGCTAG--------CTGATCTATCGCTCAGTAGGACCAACGTA--------CACGCT-AATAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCAAAAATGTATGTCAGGGGCCCAG +>T47 +ACGAATTT-TCCGAGGTTACCGGATCACGATA---CGTCTTCATCGCTG-TCTGTGTACTTCCCT--------GGAAT-------------------------------------------------------------------CGATGGTCCAAATAAGAGACCATAA----GAGCC-CCG-------------ACCAGGGGAACTGACCCTCTCGAACCGAGATGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACTACTTCTAT--CC-------ACTC-GTTCTGCCCATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGCGGCCTGGTGGGAGGGGCACCACGGTCGGAGCGAATGCAACTTTAAGCCGTGCCGCGTCCAGGCCA---GGTATTCAGCGAGACACGCTACGTCCTCGTAACTCCTCAGC--------------TATGA-------------------------------------------------------------------------------------------------T-----GATCGA---A---------ATAGAGGGAACTAGGGA-CATTGC-GACAGGT-AGAAAGCGGATACGAATGCTAGAAGGATCGACCCCTGGAAGCAGACAAGACGACCTAGCGCTTACTAGCGCGACCGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGAGGCCCAG +>T8 +AATAATTT-TCCGAGGTTACCGGATCACGATA---CGTCTTCATCGCTG-TCTGTGTACTTCCCT--------GGAAT-------------------------------------------------------------------CGATGGTCCAAATAAGAGACCATAA----GAGCC-CCG-------------ACCAGGCGAACTGACCCTCTCGAACCGAGATGAATGCGCATCGCG-----------C-TAGCTTCAC-CATGTACTACGTCTAT--CC-------TGAC-GTTCTGCCCATCGACGCAGAATCACGAAATT---CACAACTGAGTATCTGCGGCCTGGTGGGAGGGGCACCACGGTCGTAGCGATTGCAACTGTAGGCAGTGCCGCGTTCAGACCA---GGTATTCAGCGAGCCACGCTAGGTCCTCGTAACTCCTCAGC--------------TATGA-------------------------------------------------------------------------------------------------T-----GATCGA---A---------ATAGAGCGTACTAGGGA-CATTGC-GACAGGT-AGAAAGCGGATACGAATGCTAGAAGGATCTACCACTGGAAGCAGACAAGACGACCTAGCGCTTACTAGCGCTACCGTA--------CCCGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGGATGTCAGAGGCCCAG +>T33 +AACAATTT-TGCGAGATTACAGGATCACTATA---CTTCTTCATCCCTA-TCTGTGTACTACCCA--------GGAAT-------------------------------------------------------------------CGATGGTCCAAGTAACGGACCATAA----GAACG-CCG-------------ACGAGAGCAACTGACCCTTTCGACCCGAGACGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACGGATTCTAT--CC-------ACTC-GTTCTGCCAATCGACGGACAATGACGAAAAG---CACAACTGAGTATCTGCGGCCTGTTGGGAGGGGCATTTAGGTCGTAGCGAATGCAACGGTAAGCCGTGGCGCGTCTATGCTA---GGTATTCAGCGGGACATGCTATGTCCTCTTATCTCCTCAGC--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCGA---A---------ATAGAGGGAACCAGGGA-CTTTGA-GACAGGA-AGGAAGTGGATTCATATGCTTGAAGGATCGACCCCTGGTAG--------CCGATCTAGCGCTTAATAGCCGTACCGTA--------TACGCT-ATTAGACGC-TAGAT--------T----TAAATTC------CC-TC-AAAGTCATATATGTATGTGAGGGGCCCAG +>T4 +ACCAATTT-TGCGAGATTACAGGATCACCATA---CTTCTTCATCGCTG-TCTGTGTACTACCCA--------GGAAT-------------------------------------------------------------------CGATGGTCCAAGTAACTAACCATAA----GAACC-CCG-------------ACGAGAGGAACTGACCCTCTCGACCCGAGAAGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACAGCTTCTAT--CC-------ACTC-GTTCTGCCAATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGCGGCCTGGTGGGAGGGGCATTAAGGTCGTAGCGAATGCAACGGTAAGCCGTGGCGCGTCTATGCTA---GGTATTCAGCGAGACATGCTATGTCCTCGTATCTCCTCAGC--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCGA---A---------ACAGAGGGAACTAGGGA-CTTTGA-GACAGGA-AGAAAGCGGATACATATGCTTGAAGGATCGACCCCTGGTAG--------CCGATCTAGCGCTTACTAGCGCTACCGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGGGGCCCAG +>T6 +AGCACTTT-GGGGACGTAACAGGATCACCGAT---CGTCTGAAGCGCTG-TGGGACTACGAGCCA--------GCAAT------------------------------------------------------------------TGGATCGTGCAACTAACTGACCATTA----GAAAG-CCG-------------GCGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCTCAGCG-----------C-AAGCTTCAC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCAAGCGACGCACAATCACGAAAAT---CGCCACTTAGTTTCTTCAGCCTGGTGGGAGGGGCACTAAGACCGCAGCGAACATAATGGTAACCCGCGGCCCGTCTATGCTA---GGCGTGCAGCGAGACATCCTATCGCTACGTATCCACTCAGT--------------TTTGG-------------------------------------------------------------------------------------------------T-----GCACGA---G---------ACAGAGTGAACCAGGCA-GTTTAA-GACATGA-AGAAAGCGGATACATATGCTTGAACAATAGACCCGTCGTAG--------CTGATCTAGAGCTTAGCAGCGCCATCGTA--------CACGCTTATGAGTTGT-TAGAT--------T----TAAGTCC------CC-TT-AAAGTCACATATGTATCTCAGGGGCCCAG +>T57 +AACTCCCT-GGGGTAGTACTAGGATCACCGAT---CTTTGGAAAGGCTA-TGTGTGAACTAGCAA--------ACAAG-------------------------------------------------------------------AGACCCTGCATCTAACTGTCAAGTA----CTAAT-CCG-------------TCGACGGCAAATGAGCGTATCGCCATGTGTATTAACCGATCCGCT-----------T-AGGCTTCTC-CATGGATCGCCTCTAG--CT-------AATC-GTCTGGCCTAACGATCGTCAAGTACGAACAT---CGCCACATTGTGTTTTGA-GAACGTGGGAGGCCCACTAAGAGCGCAGCGCACGAAATGGTTATACGAGGTCTCTCGTTGCTC---TGGGTGCAGAGAGAGTTAGTATCGCTAGGTATCTACGCAGT--------------TATAC--------------------------------------------------------------------------------------------C----T-----GGATTA---T---------ACTGTGTGAACGACGCA-GTATAA-GCATCGA-AGCAAGCGGATACAGCTGTTTGTACAATCCAACATTCCTAG--------CCGAACTAAAGTTCAGTAATGTCATCGC---------CCCGCTTACGAGCTGC-TAGAC--------T----CAAGTCC------CCGAG-CCAGTCCAATATTTGTCCCAGGTGCCGAG +>T46 +AGCACCCT-GGGGTAGTACCAGGATCACCCAT---CTTTGGAAAGGCTG-TGTGAGAACTAGCAA--------ACAAG-------------------------------------------------------------------AGATCCTGCATCTAACTAACAAGTA----GTAAG-CCG-------------ACGACAGCAAATGAGCGTATCGCCCGGTGTATTAATCGCTCAGCT-----------C-AGGCTTCTC-CATGGATCGCCTCTAT--CT-------AATC-GTCTGGCCTAACGATCCTCAAGTACGAACAT---CGCCTCATAGTGTTTTTA-GAAGGTGGGAGGGGCACTAAGAGCGCAGCGAACGAAATGGTAACACGAGGTCTCTCGGTGCTC---TGGGTGCAGAGATACATAGTATCGGTAGGTATCTACGCAGT--------------TATAC-------------------------------------------------------------------------------------------------T-----GGATCA---T---------ACTGTGTGAACGACGCA-GTATAA-GCATCGA-AGCAAGCGGATACAGCTGTTTGTACAATCCTACAGTCCTAG--------ACGAACTAAAGTTCAGTAATGACATCGC---------CCCGCTTACGAGCTGC-TAGAC--------T----CAAGTCC------CCGCG-CCAGTCCGATATTTGTCTCAGGTGCCGAG +>T23 +AGCAACTT-GGGGTAGTAACAGCATCACCTAT---TTTTTGAAAGGCTG-TGAGAGAACTAGCAA--------GCAAG-------------------------------------------------------------------GGATACTGCAGCTAACTGTCAAGTA----GTAAG-CAG-------------ACGACAGGAAATGAGCGTATCGCCCGGTGTATTAAACGCTCAGCT-----------C-AGGCCCCTC-CATGGAGCGCCTCTAT--CC-------AATC-GTCTCGCCTAACGATCCACAAGTACGACCAT---CGCCACTTAGTTTTGTTAGGAAGGTGGGAGGGGCACTAAGACGGCAGCGAACGAAATGGTAACACGAGGTCTCTCGGTGCTC---TGGGTGCAGAGAGACAGCGTATCGCTAGGCATCTACACAGT--------------TATCC-------------------------------------------------------------------------------------------------T-----GGACCA---T---------ACTGTGTGAACGACGCA-GTATAA-CAAACGA-AGCAAGCGGATACAGCTGTTTGTACAATCCATCAGCCCTAG--------CCAAACTAAGGGTCAGTAATGACATCGC---------CCCGCTTATGCGCTGC-TAGAC--------T----TAAGTCC------CCGTG-AAAGTCCGATATTTGTCTCAGGTGCCCAG +>T30 +AGCACTTT-GGGGTCGTAGCAGGATCACCGAT---CTTGTGAAGGGCTG-TGGGAAAACTAGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCAACTAACTGACCAGTA----GAAAG-CCG-------------ACGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCCCAGCA-----------C-AGGCTTCTC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCTAACGATGTACAATTACGAAAAT---CGCCACTTAGTTTCTTCTGCCTGGTGGGATAGGCACTAAGACCGCAGCGAACGAAATGGTAACCCGGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGCTACGTATCTACTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GGACGA---T---------ACAGAGTGAACCAGCCA-GTTTAA-GAAACGA-AGAAAGCGGATACATCTGTTTGAACAATAGACCCGTCGTGG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CCCGCCTATGAGGTGT-TAGAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATTTATCTCAGGGGCCTAG +>T24 +AGCCGGTT-GGGGCCATAGCAGGATCACAGAC---CTTGCGATGATCTC-TGGGAAATCTGGCTA--------GCAAG------------------------------------------------------------------TGTCTCGTGAAATTACCTGACCAGTA----GATGA-CCG-------------CGGATAGTAACTGACTCTCTCGCGCCGGGAACTCTTCTCGCAGCA-----------C-AGGCGTCTC-CAGCCATAGCATCGCT--CA-------AGTC-GGCTGGCCTAACCCTGCACGCTCACGAAAAT---CGCTACTTATTTTCTTTTGCCCCGTAGGAGACGCAATAAGGTCGCAGCGACCCAAATGGTAACCCTGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGAGCCGATTCGACTCTGT--------------TATGA-------------------------------------------------------------------------------------------------T-----CGACGG---C---------AAAGAGTGCACCAGGCT-GTTTAA-GAAACGA-AGAATGCCGATACATCTGTTCGAACAAAAGTGTCCTTGGAG--------CTGAAATATAGATTAGCAGCCTCATCGAA--------ACCGCTTGTGAGGTGT-TCTAC--------G----TAAGTCC------CC-TT-AAAGTCAGACATCTATCTAAGGGGGGTTG +>T41 +AGCCGCTT-GGGGCCATACCAGGATCAAAGAT---CTTGTGATGATCTC-TGGGAAATCTGGCTA--------GCAAG-------------------------------------------------------------GCCCATGTATCGTGAAATTACCTGACCAGTA----AATGC-CCG-------------CGGATAGTAACAGACTCTCTCGTGCCGGGAAGTCTTCTCGCAGCA-----------C-ACTCGTGCC-CAGCCATAGCATCGCG--TA-------AGTC-GGCTGGCCTAACCCTGCACGCTCACGAAAAT---CGCTACTTATTATCTTTTGCCCCGTAGGAGACGCATTTAGGTCGCAGCGACGCAAATGGTAACCCTGGGCCCCTCTATTCTA---GGCGTGCAGAGAGACTTCCTATCGATACGATTCGTCTCTGT--------------TATGA-------------------------------------------------------------------------------------------------------CGACGC---C---------AAAAAGGACTCCAGGCT-GTTTAA-GAAACGA-AGAATGCCGCTAAATCTGTTCGAACAAAAGTGTCCTTGC-G--------CTGAAATAGAGCTTAGCAGCCTCATCGAA--------ACCGCTTGTGTGCTGT-TCTAC--------T----TAACTCC------CC-TT-AAAGTCAGATATCTATCCAAGGGGGCTTG +>T32 +AGCCCTTT-GGGCGCATAACAGGATCACGGTT---CCTGTGAAGCGATC-TGGCAAAGGTGGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCACCTGCCTGAACAGTG----GATAT-CCG-------------AGGAGAGTTGCTGACTCTTTCGCGCCGGGAAGTCCTTTCCAAGCA-----------C-AGGCTTCTC-TATGCATGGCATCGCT--TA-------TATT-TACCTGCCTAGCCATGCACCCTGACGAACAT---CGACACTTTTTTTATTCTGCCTTGTGGGTGAGGCTCTAAGATCGCAGCGACCGGACTTGTAACCCGGGGCCCCTCTATGATA---GACATGCAGAGAGACATACTATGGTTACGTATCCAAAATGT--------------TATGT-------------------------------------------------------------------------------------------------T-----GTACGG---C---------ACAAAGTCAACCAGGCAAGTTTTA-GAAGCGA-AGAAGGCCGATACATCTGTTCCCACAACAGACCCCTGGGAG--------GTGAACTAGAGCTTAGCAGCGCCATCGTA--------CTCGCTTATGAGGTGT-TCGAC--------T----CAAGCCC------CC-TT-AAAGTCAGATCTTTATCTCAGGGTTCTTG +>T52 +AGCCCTAT-TCGCGCTTCACCGTTTCACTGTA---CATTTGACGGGATC-CGCCAAAGCTGCCCA--------GCAAG------------------------------------------------------------------TAGATCGTGCTTCTGCCATAACGGTG----GACAT-CCT-------------AGGAGAGTCGCTGAATCTTTCGCCCAAGGAGGTTCTTCCATACCA-----------A-TTTCTTCTA-TATGC-TGGGATCGCT--TA-------TCTA-GCCCTGCACTTTCCTGCGGCCTGAGGAAAAT---CGACACGTAATTTGTTCTTCCTTGTGGCAGAAGGTTTAAGCTCGTAAGAACCGGCCTTTTACATAGGGGCCCCTCCATGTTG---GTCATGCTAAGAGTCAT-----GATTGCGTATACACAATGT--------------TGATT-------------------------------------------------------------------------------------------------T-----GTTCGG---A---------ACAATCTTCACATGGCAAGTTTCA-GAAGCGA-AGTAGGCAGATACATCGATCCCGACGATGGACCCCGGGAAG--------GTGAACTAGAGCATAGTAGCTCTTTGGTA--------CTCGCTTACGCCTTAT-TAGAC--------T----CAAGCCA------CC-TG-AACCTCCGAGAACTATCTAAGCGTTCTTG +>T14 +AGCCCTAT-TCGCGGCTCACCGTTTCACTGTA---CCTTTGACGGGATC-CGCCAAAGCTGCCCA--------GCAAG------------------------------------------------------------------TAG-TCGTGCTACTGCCATAACGGTG----GACAT-CCT-------------AGGAGAGTCGCTGAATCTTTCGCCCAAGGAGGTTCTTCCCGACCA-----------A-TTTCTTCTA-TATGC-TGGGATCGCT--TA-------TCTA-GGCCTGCAGTTCCCTGCGGCCTGAGGAAAAT---CGCCACGTAATTTATTCTTCCTTGTGGCAGAGGGTTTAAGCTCGTAAGAACCGGGCTTTTACATAGGGGCCCCTCCATGTTG---GTCATGATAAGAGTCAT-----GATTGCGTATACACAATGT--------------TGATT-------------------------------------------------------------------------------------------------T-----GTTCGG---A---------ACAATGTTCACATGGCAAGTTTCA-GAAGCGA-AGTAGGCAGATACATCGATCCCGACGATGGACCCCGGGAAG--------GTGAACTAGAGCATAGTAGCTCTTTGGTA--------CTCTCTTACGCCTTAT-TAGAC--------T----CAAGCCA------CC-TC-AATCTCCGAGAACTATCTAAGCGGTCTTG +>T10 +GGCTCTTC-CCGCGCCACACCCCTTCGCTGTA---CCTGCGACCGGATT-CGAGCAAGGTGGCCT--------ACAAG---------------------------------------------------------------------------------------------------------------------------------------------TCGAGGAGGTTCTGCCCGATCC-----------A-CCTCCCCTA-TATTC-CGGGAGCGCT--TA-------TGTA-GGCCTGCAGGTCCCCGCGGCAGGAGAAAAAT---CGCCATGGAATTTAATCTAGGTTGTGGCAGAGGGATTAGGCTCGTTAGATCCGAGCTACTACAGATGGGCCCCGATTTTGTG---GTAATGATAAGATTCAT-----TACTTCGTACATACAAGGT--------------TGATT---------------------------------------------------------------------------------------------GTACT-----GTGAGG---A---------ACAATTTGTACATGGAACGGTGAA-GAAAGGA-GGTAGTAATATACCGCGGTCCCGACGATGGAGCCCCGAAAG--------GTGCACTCGAGCATAGTAGCTGTTTGGTG--------CCCTCTTACGCCCTAT-TAGAC--------T----CAACCCT------AC-TC-TATCTCCGAGAACCCTCTAACCGCTTTTG +>T42 +GGCCCTAT-TCGCGCCACACCGTTTCACTGAA---CCTGTGACCGGATC-CGCCCAAGGTGGCCT--------ACAAG-------------------------------------------------------------------------------------------------------------------------------------------GCCCAAGGAGGTTCTTCCCGATCC-----------A-TCTCTCCTA-TATGC-TGGGATCGCT--TA-------TCTA-GGCCTGCAGTTCCCTGCGGCCTGAGAGAAAT---CGCCATGGAATTTGATCTACGTAGTGCCAGAGGGCTTAAGCTCGTAAGATCCGAGCTATTACTTATGGGCCCCTATTTGTTG---GTCATGATAATACTCAT-----TAATTCGTATATACAACGT--------------TGATT---------------------------------------------------------------------------------------------GTACT-----GTTCGG---A---------ACAATGTTCACATGGAAAGTTGCA-GAAGCGA-AGTAGGAAAATACCTCGGTCCCGATGATGGAGCCCGGAAAG--------GTGAACTCGAGCATAGTAGCTGTTTGGTA--------CTCACTTACGCACTAT-TAGAC--------T----GAACCCA------CC-TC-TATCTCCGAGAACCATCTAAGCGGTCTTG +>T5 +TGACCGTACACGTGAATGGCACGATCCCTGATATTTCTGAAAGAGGCGG-TAGCTGCCAAGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGTCTAACAGACT---------------------TGGCAGAA-----------C-GGTTTAGTG-TATGCATAGAGGCTGT--CC-------TGTA-CGACCGTCTAATATCATGCACCCACGTCCAT---GAAGAATTAGTTGCATCTCCGTGGTTAGTATGGCACAACAGATGCGGCAGCAGGGAGTGGAAGGTGGGTGCTTTCCCCGAAA---GGTGCGCAA-GAGCAAGCCTGCCGTGATGTATACGTTCTAT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATAAAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-TGAATGC-GAACCGCCCCTTCATAGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGTTGCATTCTA--------CCCGCCTATTGCCCCTATTGAT--------T----AAAATTA------CG-TC-GGAGCCAGACTCGTGTATCGATGTACTAT +>T21 +TGACCGTACACGTGAATAGCACGATCCCTGATATTTCTTAAAGAGGCGG-TAGCTGCCAAGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGTCTAACAGACT---------------------TGGCAGAA-----------C-GGTTTAGTG-TATGCATAGAGGCTGT--CC-------TGTA-CGACCGTCTAAAATCATGCACCCACGTCCAT---GAAGAATTAGTTGCATCTCCGTGGTTAGTATGGCACAACAGATGCGGCAGCAGGGAGTGGAAGGTGGGGGCTTTCCCCGAAA---GGTGCGCAA-GAGCAAGCCCGCGGTGATGTATACGTTCTAT--------------TGTAG-------------------------------------------------------------------------------------------------T-----ATAAAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-TGAATGC-GAACCGCCCCTTCATAGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGTTGCATTCTA--------CCCGCCTATTGCCCCTATTGAT--------T----AAAATTA------CG-TC-GGAGCCAGACTCGTGTATCGATGTACTAT +>T39 +TGTCCGTACACGTGGCTGGCACGATCCATGATATTTCCGGAACAAGCGG-TAGCTGCCAAGACCG--------GGCTG-------------------------------------------------------------------TGAACATACGACACGCTGGCCAGTTCACCGAAAG-CGA-------------ACTAGACTAACAGACT---------------------TTGAAAAA-----------C-GGTTAAGTG-TATGCATAGGGGTTGT--GC-------GGTA-CGAAAATCTAATGTCATGCACCCACGACCAT---GAAGAATTGGTTGCATCTCAGTGGTTAGTATGTCCCAACCAATTTGGCAGCAGGGAGTTGAAGGTGGGTGCTAACGCCGCAAGACGATGCGTAA-GGGCAAGCCCGCCGTGATGTATTCGTTCTCC--------------CGTAG-----------------------------------------------------------------------------A-------------G-----T-----ATATCATCAC----------------GAGCCCACAT-GTTTAG-GGAGTGC-GAACAGCCCCCTCTTTGCTTGCGGTACGAGTACACTTCGAT--------CTAATCGACGT-GTTCGAGTTACTTTCTG--------CCCTCCTATTGCCCCTATTGAT--------T----CAGATTA------CC-TG-GGAGCAAGACTCGTATATCGATGTACTAT +>T48 +GGTCCGTCATCGTGGAAGGCTCGCTCCCTCATAGTTCTTGGATAGGTGG-AACCTGCCAAGAGCA--------GTCTG------------------------------------------------------------------ATAAACATACAACAAACTGTTCAGATCAATGAGAG-CGA-------------AACAATCTGACAGACG----------------------CGCAGAA-----------C-GGCTACGTG-ATTGGATCGTGGTTGT--CG-------TGTA-CGACCGTATAATATCATGCACCCATGACCCC---GAAGAATTAGTTGAAGCTCGCTGGTAACAATAGAACATCGATTGACCCGGCAGGGAGTGGAACATCGGTGCTTTTACTGAGA---CGTACGGAA-GGGCCAGGCTTCCGTACCGTTTACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----AAACAATCAC---------ATGAATCGAGCGAACAT-GTATAG-GGAGTGC-GAAGCACCCGTTCCTGACTTACGGTAGGATAACACTTCGCT--------CTGGTCAGCTC-GATCGAGCAGGATTCCA--------CGCGCCGATTGCCGGAATTGAT--------T----AAAATTA------AG-AC-TGAGTCAGAGTGTGAAATGGATCCACTAT +>T36 +AGTCCGTAATCGTGGAAGGCTCCCTCCCTGATATTTCTTGGATAGGTGG-AACCTGCCAAGACCA--------GTCTG------------------------------------------------------------------ATAAACATACGACAAACTGTTCAGGTCAATGAGAG-CGA-------------AATAATCTAACAGACG---------------------TTGCACAA-----------C-GGCTACGTG-ATTGGATCGTGCTTGG--CG-------TGTA-CGACCGTATAATA-CATGCACCCATGACCCC---GACGAATTAGTTGAAGCTCGCTGGTTACTATAGCACATCGATTGACCCGGCAGGGAGTGGAACATCGGTGCTTTTACTGAGA---TGTACGCAA-GGGCCAGCCTTCCGTAAAGTTTACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----AAACAATCAC---------ATGCATCGAGCCAACAT-GTTCAG-GGAGTGC-GAACCACCCCTTCCTGACTCACGGTAGGATAACACTTCGCT--------CTGGTCAGCTT-GATCGAGCAGGATTCCA--------CCCGCCGATTGCCGGTATTAAT--------T----AAAATTG------AG-GC-TGAGTCAGAGTGTGAAATCGATCTAATAT +>T54 +GGTCCGTAATCGCGGAAGGCTCACTCCCTCATATTTCTTGGATAGGTGG-AACCTGCCAAGGCCA--------GTCTG------------------------------------------------------------------ATAAACATACGACAAAGTGTTCGGGTCAATGAGAG-CGA-------------AATAATCTAACAGACG---------------------TTGCGGAA-----------C-GGCGACGTG-TTTGTATCGAGGTTGT--CG-------TGTA-CGACCGTATAATATCATGCACCCACGACCCC---GACGAATTAGTTGAAGCTCGTTGGTAATTCTAGCAAATCGATTGACCCGGCAGGGAGTGCAAGCTCGGTACTTTTACTGAGA---CGTATGCAA-ATGCCAGCCTTCCGTCACCTTTACGTTCTCT--------------TTTAG-------------------------------------------------------------------------------------------G-----T-----AAACTATCAC---------ATGGATCGATCCTACGT-GTTTAG-GGAGTGC-GAACCACCCCTTCCTGACTTACCGTAGGATAACACTTCGCT--------CTGGTCAGCTT-GATCGAGCAGGATTCCA---------ACGCTGATTGTCGGTATTGAT--------T----AAAATTA------AG-AC-TGAGTCAGAGTCTGAAGTCTATCTACTAT +>T13 +TTGCCGAGTAAGTGGATCGATATCTCCAAAAGATACGTAGATAAGGTGT-TAGCTGCCGTAATCG--------GCCTT------------------------------------------------------------------AGGAACTTACTTGAGGGTTATCTGCCCACGGAACG-CGA-------------ACAAGTCTAACAGGTC---------------------TAGGAGAA-----------C-GATTTATGG-TATGGATGGGTCCTGA--CC-------AGTA-GTAAGGTATAATGTA-CGCGCATTCGGGCAA---GCATAACGCGTTGCATCTCAGTGGTGATTCTCGCACAAAAGATGCGGCATTAGGGGGTGGAAGGAGGCTGCTTTCGGC---A---TGTGGGC-A-TAGCATGGCTGCCGGGACGTATTCCTTCACA--------------TCTAG------------------------------------------------------------------------------AATAAAAGGCAGGG-----T-----ACCCAATCAT---------ACGAACATAGACAACAA-ATGTTG-TGAGCGA-TCCCAGCACCTTCATGGGGACCGGTAGGAGATCACTGCGAT--------CTATTGAACGT-GTTCGAGCTGCAATCCA--------GCCGCCTACTACTCGCATTCAT--------T----AGAAGAC------AT-TC-GGAGCAAAACTCTTTTATCGATCGACTAT +>T31 +TGGCCGTTCATGTGGATCGTAATCTCGACGAGATTTCTAGAAAAGATGT-TAGCTGCCCTAAGCG--------GCCTG------------------------------------------------------------------AAGAACTTACAAAAAGCCTATCAGTCCAATGCACG-CGA-------------ACCAGCCTAACAGGCC---------------------TCGCCGAA-----------C-GATTTCTCG-TATGCATAGGTGCTGA--CC-------TGAA-GGACTGTCTAATATC-TGCACAATCGGGCAC---GCTCAACGCGTTGCATCTCAGTGGTGATTCTCGCACAAAAGATGCGGCATTAGGGGGTGAAAGGAGGCTGCTTTCGGCGAAA---GGTGCGC-A-GAGCATGCCTGCAGTGCTGTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----AACCAATCAT---------ACGAACAAAGACAACAT-GTTTAG-GGAGCGC-GAACCGCCCCTTCATGGCGAACCGTAGGAGATCACCTCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCCA--------CCCGCATACTTCACGCATTGAT--------T----AGAAGAC------TT-TC-GGAGCATAACTCTTATATCGATCTACGAT +>T37 +TGGCCGTTCACGTGGATCGTATTAATCGCGAGATTTCTCAAATAGATGG-TAGCTGCAATGACCG--------GCCTG------------------------------------------------------------------ACGAACATACAACAAGCTGGTCACTCCACTGAACG-CGA-------------CAAAGGCTAACAGACC---------------------CCGGAGAA-----------C-GAATTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGAATGTCTCATATC-TGTACAATCGGGCAC---GCACAACTAGTTGCATCTCAGTGGTGATACTGACTCAAAAGATGCGGCATCAGGGGGTGGAAGGAGGCTGCATTAGGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTTACGTATGCCTTCACTGGATAGCCCTGCAATGTAC-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAT---------ACGAACATAGACAACAT-GTTTAG-GCAATAC-GAACCGTCCTATCATGGGGAGCGTTAGGAGAACACTGCGAT--------TTAATCAACGT-GTTCGAGCTGCATTCCA--------CCCGTCTACTGCCGGCATTGAT--------T----AGAAAAC------AT-TC-GGAGCAATACTCTTATATCGATCGACTAT +>T22 +GGGCCTGTCTCGCCCGTCGCAACGTCCACCTACTTTCTACAAAAATTAG-TAGCTGTCTCCACCG--------GTCAG-------------------------------------------------------------------TGCACA-AAGACGAGGTGATCAGTAGAGTCACTA-CGA-------------CCAGGCCTAGCGGACT---------------------TCTGCGAA-----------C-GCTTTATTG-TAGACGGCGAGGCAGA--CA-------TGTC-CGACTGTCAAATGTC-GTCAGCTTGGGCCAC---GCACAAGTCGTCGCTTGTCAGTTGTCGTTCCGACACATAAGCTATAGAACCAGGAGGTAGGCGGA--TGGGTGTCCGCGAAT---GGGGTGC-A-GAGCAAGCCTATCGTGACCCATTCCTTCGCT--------------AGATG-------------------------------------------------------------------------------------------G-----T-----ATCCCTTAAA---------ATGCATCATGACCAAGC-GTTTAG-TGAGCTC-GAACGGCTCCTTCCTGTCGAGCGCGCGGAGAACACAGTC-------------TTCAACGC-GTTCATGAG---------------------------------------------------AGATGGC------AT-AC-AGAGCAAACCCTTATGATCGGTCTCCTAG +>T40 +TTGGTTGACCTACCATTCTTCAAATGGTGCGTGATCAGGCTAGTTTTAG-TGCTTGACTCCC--G--------GGCAG--------------------------------------------------------------------CTACT-GGAAAGGTCAAATGAATACATTCAATA-TCA-------------CTAGGCCTTGGGGGCA---------------------TCTCCGATACA--------T-ACTTTGATA-TGAGCACCAAATATAG--TC-------TCAG-CCGCGGTTGGATAAC-TTTACCTTTGGCCTA---CCACCAATTGTGGATAGGATCCTTATCTGGTGGATCCAATTTCGGAGACAAGGGCGGTGGG---AGTCCAATGTCTCGG------GTACAGG-A-GGGGGGCATCA-CCTGTACCGTTCGTGCTCC--------------ATTCC-------------------------------------------------------------------------------------------C-----C-----TTCGCATAGA------------CACAGAGACGAACT-GCGTGT-TGCCCGT-CAACTGGCCCTGCATTGCGAGAGCGATGGGGCTTCAGTG-------------TTCGT--T-GTGCGCTCG---------------------------------------------------AGATCGC------CC-AGTACGTCTAACCCTTGAGTTCGGATTCCAAG +>T44 +TGGGTTGACCTACCCGTCTTCGTACGCTGCATCTTCCAACTAGTAATAG-TGGTTTACTCCC--G--------GGCAG--------------------------------------------------------------------CCACT-GTAAAGCTCAGATGAGTAGATGCACTA-TCA-------------CCAGGCCTAGCGGGCA---------------------TCTCCGACACC--------C-ACTTTTTTA-TGAACAACGAGGCTGG--GG-------GCTT-CGGCGGTAGAATAAC-GTCACCTTGGGCATA---TCACCAATCGTGGCTTGGATTTTTCTCTGCTGACTCAAAATACACAGAAACGGGCGGTAGG---AGACCGATGTCTCGGTAT---ATTGAGC-A-GGGGCTCATTAACCTGAACCGTTACTGCTGC--------------AATAG-------------------------------------------------------------------------------------------G-----C-----TTCACATAGA------------CACAAAGACGAACT-GATCAT-TGCCCGT-CAACTGCCCCTTCATGGCGAGAGCGAGGAGGATACATTG-------------ATAAT--C-GTGCGCTCG---------------------------------------------------AGGTCGC------AC-AGTACGGGCCACCGTTAAGTTAGGCTTCCTAG +>T25 +TGGCCCTTCAAGTTGATCGTTTTACCCTCTAGCTTTCTAGAAATACCAC-TTGCTGCTTCGGGGG--------TGCAG------------------------------------------------------------------ATGATCAATCAAAAATCTGATCGGTCTTCCCGCGC-CGA-------------GCAGGGCTAACAGACT---------------------TCGCAGCA-----------C-AATTTAGTG-TATACATGGGTGCTGA--CT-------CGCA-CGTTTCACTAAGACA-TGCACCCACGGCCTC---GCACAACTCGTTTCTTCTCAGTGGTGATCCTGACACCAAGGCTGCCCCCCCAGGCGGTGGAAAGAGGTTACTGACCGCTAAC---GCAGCGC-A-CAGGAAGCCTGCCTTGACCTATTACTTCACC--------------TATAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGATCTCAACCGCCAT-GTTTGG-TGCGCGC-GAACAGGCCCTTGCCATCGTGCGTTAGGATAACACTGCAGT--------CTAATCAACGT-TTTCGATCT-------------------------------------------------CTAGAAATC------AT-AC-GGAGCAATACCCTTAGATGAGTCGCCTAG +>T51 +TGGCCGTTCACGTTGATCGTCTTATCCTCGTGCTTTCTAGAAAAGCTAG-TAGCTGCCACGAGCG--------TCCTG------------------------------------------------------------------ATGAACACCCAACAATCTCATGAGTCCACTGACTC-CGA-------------GTAGGGCTATCAGACT---------------------TCGGAGCA-----------C-GATTTATTG-TATACATGAGTGCTGA--CC-------TGCA-CGCTTGCCTAAGATC-TGCACCCACAGCCAC---GCACAACTAGTTGCATCTCAGCGGTGATTCTGACACAAAGGATGCGGCATCAGGCGGTGGAAAGAGGTGACTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTACTTCACT--------------TCTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTCAAACGCCAT-GTTTAG-TGCGCGC-GATCCACCCCTTGCCGTCGAGCGTTAGGAGAACACTGCAAT--------CTAAGCGACGT-GTTCGATCT---------------------------------------------------AGAAGAC------TT-AC-GGAGCAATACCCTTAGATGGGTCGCCTAG +>T50 +TGACCGTTTACGTGGATAGTAGGATTCACGAGCTTTCTACAAAAGAGAG-TAGCTGCCACGACCG--------GCCTG-------------------------------------------------------------------TGAACATACCACAAGCTGAGCAGTCCCCTGAAAG-CGA-------------ACAATCCTAACAGACT---------------------TCGGAGAA---GCCGGTAGC-GATTTATGG-AATGCATGGGGGCTGA--CC-------TGAA-CGACCGTCTGATATC-TTCACCTTCGGCCAC---GCACAACGAGTTCCATCTCAGTGGTGAGTTTGACACAAAAGATGCGGCATCAGGGGGTG---GGAGGTTACTTTCCGCGATA---GGTTAGC-A-GAGCAAGCCTGTCGTGTTCCATCCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------------T-----ATCCAATCAC---------ACGAACTTAGACACCAT-GTCTAG-TGAGCGC-GAACCGTCCCCCCAGGTCGATAGGTAGCAGACCACTATAAT--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGGGAAAGACTCTTATATCGGTCGCCTCT +>T29 +TGACTGTTCTCGTGGATTGGATGAGTCGCGAGCTTTCTATAAAAGAGCG-AAGCCGCTACGACAG--------GCCTG-------------------------------------------------------------------TGCACATACCACAAACTGAGCAGTTCATTGATAG-AGA-------------ACAAACCAAACAGCTC---------------------TCGGAGAA---CCCGCTGGA-CGTTTTTGG-AATCCATGGTGGATGA--AG-------TTAA-ACACGGTATGTTGTT-TCAACCTTCCGCCTC---GCTGAACGGGATCCATCTCAGTCCTGAGTATGCCG--AAAGCTGCAGCTTCAGTGGGTG---GGCCGTGACTTTCCGCGAAA---GGTCAGC-A-GACCGAGCCTTACGTGTTCCAGAACTTTACT--------------TGTAG-------------------------------------------------------------------------------------------------T-----AACCAAACAC---------ACAACCTTCGTGACCAT-GTCTAG-TGAGCGC-GGATCGTCCGGCAAACTCGATTTGTCGCAGGCAACTAGAAT--------CTAATCCCCAT-GTTCGGCTT---------------------------------------------------AGAAGAC------AT-TC-GGGGAAAGGGAGTTATATCGGTCTCATCT +>T56 +TGGCCGTTCACGTGGCTCGTAGCAACCACGAGCTTTCGAAAAAATATAG-TAGCTGCCAGGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCCCTGAACG-CGA-------------ACAAGCCCAACAGACT--------------------------AAA-----------C-GATTCATGG-TATGCATGGGTGCTGT--CC-------TGAA-AGACTGTCTAATATC-TGCCCCCTCGGGCAC---GCACAACTAGTTGCATTTCAGTGGTGATTCTGACACAAAGGAAGCGGCATCAGGGGGTGGAAGGTGGTTGCTTTTCGCCATA---GTTGCCC-A-GAGCAATCCTGCCGTGACCTATCCCTTCACT--------------TGTAT-------------------------------------------------------------------------------------------G-----T-----ATCCACTAAC---------ACGAACTTAGACAACAT-GCTTAG-GGACCGC-GTACCGCCCCTTCATGTCGAGCGGTAGGAGAACACTGCAAC--------ATAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGAGCAAGACTCTTATATCGGTCGCCTAT +>T1 +TCAGCGTTAACTTGGAACGTACCATCCAGGAGCTTTCTAGAAAAGATGG-TGGCTGCCAGGACGG--------GCCTG------------------------------------------------------------------ATGAACATACAATAGGCTGGTCGGTACACTGATCG-CGA-------------ACAACCCTAACAGACT---------------------TCGGTAAG-----------C-GATTTATAG-TATGCATCGCTGCTGG--CG-------TGAC-CGACTGTATAATATC-TTCCACGTCGGGCAC---GCACAACTAGTTGCATCGGAGTGGAGATTCTGACACAAAGCTTACGGCTTTGGGGTGTGCGAGGTGTATGCTTCTCGAGTAT---CGTGCGC-A-GGGCAAGCCTGGCTTAACCGCTTCTTTTACT--------------TGTAG-------------------------------------------------------------------------------------------C-----T-----ATCCACTAAT---------AAGTACTTAGATAACAC-GTTTTG-TGACCGC-GTGCCGCCCTTTTCTAACCAGCGGTAGGAGAAGACCGCACC--------TTAATGAACGT-GCTCGAGCT---------------------------------------------------AGTAGAT------AC-TC-GGAGCAAGACTCTTATATCGGTGGCCTTT +>I62 +GGTACTTAAGCGTGGATGGCAGGATCCCTGATATTTCTGCAAATAGTTC-TAGCTACAATGACCA--------GCCTG------------------------------------------------------------------AGGAACATGCAACAATCTGTTCAGTTGACTGAAAG-CTA-------------ACTAGGCTAACTGACT---------------------TCGCAGAA-----------C-GGTTTACTG-TATGCATGGGGGCTGT--CA-------TGTA-CGACGGTCTAACATCATACACCTACGACCAC---GGACACCTAGTTGCATCTGAGTGGTAAGTGTGACACAAGAAATGCGGCACCGGGGGGTGTAATTTGGGTGCCTTCCACGAAA---GGCGCGCAGAGAGAAAGCCTGTCGTTACGTATCCATCCACT--------------GGTGT-------------------------------------------------------------------------------------------G-----T-----AGACGACCAC---------ACGAACTGAGCCAAGAT-GTTTAG-GGAGCGC-ACACCGCCCCTTCATGGGTTGCTGCAGGAGAACACTTCGAT--------CTAATCTAGGT-GTTCGAGCGGCATTCTA--------CCCGCCTATAAGCCGTTTCGAC--------T----TAAATAA------CG-TC-GGAGCCAGACTTTTATATCGATGTACTAT +>I44 +GGACCTTA-GCTGCGATGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGATCTTGCAACAACCTGATCAGTA----GATAG-CTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TAAGCATGGCTTCTGT--CA-------TATC-GGACTGTCTAACATTATACACTTACGACAAC---GGACACTTAGTTGCTTCTGCGTGGTAGGAGAGACACAAAGAACGCAGCAACCGGTATTGTAATTTGGGTCCCATCTAGGCAA---GGCGCGCAGAGAGAAACCCTGTCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ACAAACTGAACCAGGAA-GTTTAG-GTAACGA-AGACCGCCGGTACAACTGTTGCAGCAATAGACCTCTTGGAT--------CTGATCTAGCGTTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGTTTCGAC--------T----TAAGTAC------CC-TT-AAACTCAGACTTTTATCTCAATGTCTTAG +>I23 +GGACCTTA-GCTGCGATGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGATCTTGCAACAACCTGATCAGTA----GATAG-CTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TAAGCATGGCTTCTGG--CA-------TATC-GGACTGTCTAACATTATACACTTACGACAAC---GGACACTTAGTTGCTTCTGCGTGGTAGGAGAGACACAAAGAACGCAGCAACCGGTATTGTAATTTGGGTCCCTTCTAGGCAA---GGCGCGCAGAGAGAAACCCTGTCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ACAAACTGAACCAGGAA-GTTTAG-GTAACGA-AGACCGCCGGTACAACTGTTGCAGCAATAGACCTCTTGGAT--------CTGATCTAGCGTTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGTTTCGAC--------T----TAAGTAC------CC-TT-GAACTCAGACTTTTATCTCAATGTCTTAG +>I9 +GGACCTTA-GCTGCGATGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGATCTTGCAACAACCTGATCAGTA----GATAG-CTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TAAGCATGGCTTCTGG--CA-------TATC-GGACTGTCTAACATTATACACTTACGACAAC---GGACACTCAGTTGCTTCTGCGTGGTAGGAGAGACACAAAGAACGCAGAAACCGGTATTGTAATTTGGTTCCCTTCTAGGCAA---GGCGCGCAGAGAGAAACCCTGTCTGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ACAAACTGAACCAGGAA-GTTTAG-GTAACGA-AGACCGCCGGTACAAATGTTGCAGCAATAGACCTCTTGGAT--------CTGATCTAGCGTTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGTTTCGAC--------T----TAAGTAC------CC-TT-GAACTCAGAGTTTTATCTCAATGTCTTAG +>I5 +GGACCTTA-GCTGCGATGACAGTTTCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCAC--------GCAAG------------------------------------------------------------------TGCATCTTGCAACACCCTGATAAGTA----GATAG-TTG-------------ACAAGAGTCACTGACTCTCTAGCGCCGACAGGTCTTCTCGCAACA-----------C-AGAATGCTGATAAGCATGGCTTCTGG--CA-------TATC-CGATTGTCTACGATTACACAGATACGACAAC---GGACACACAGTTGCTTCTGCGTGGTAGGAGAGACACAAAGAACGCAGAAACCGGCATGGTAATTTGGTTCCCTTCTAGGCAA---TGCGCGCAGACAGAAACCCTCTCTGTACGTATCGAATCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ACAAACTGAACGAGGAA-GTTTAG-GTAACGA-AGACCGCCGGTACGAATGTTGGAGCAATAGACCTCTTGGAT--------CTGCTCTAGCGTTTCCAAGCGCCATTGTA--------CCTGCCTTTGAGATGTTTCGCC--------T----GAAGTCC------CC-TT-GAACTCAGAGTTTTATCTCAATGTCTTAG +>I0 +GAACCTTA-CCTGCGATGACAGTTTGACTGAT---CTCGGTAAGAGTTC-TAGGAAAGATGGCAC--------GTAGG------------------------------------------------------------------TGCATCTTCCAACACCCTGATGATTG----GA--------------------------TTCGCTTACTTTCTAGCTCCGAATGGTCTTCTCGCAATA-----------C-AGTATGATGATAGTCATGGCTTGTGG--TA-------GATC-CGATTGTCTACGCTTACACATTTACGTCACT---GGGCACACAGTTGCTTCTGCGTGGCAGGTCCGACGCAAGAAACGCAGAAGCCGGCATGGTAACTTGGTTCTCGTCTATGCTA---TGCGCGCAGACTGAAACCCAATCTCTACCTATCGAATCAGT--------------TGTGTCGTACGCTGCACAGTCTCTTGCACCCATGGCCTATGGTCAGATGGTATTAAGCAAAAATACCCGACACCTTTCTGCA--------------G-----G-----GGACGA---C---------ACAAACTGAACGAAGAA-GTTCAG-GTAATGA-AGCCCGCCACTACG-----TGTAGCAATAGACATCTTGGAC--------CTCCTCTAGCGTTTCCAAGCGCCGTTGTA--------CCTGCCTTTGAGTTATTTCGCC--------T----GAAGTCG------CC-TT-CAACGCAGAGTCTTATGTCAATGGCTTAG +>I4 +GCACCTTA-GATGCAATGACTGTTTGATTGAT---CTTGGTAAGCGTTC-TTGGAAAAATGTCAC--------GCGAA------------------------------------------------------------------TGCATCTTCCAACACCTTGAAAAGTA----GATAG-CTG-------------AGACGAATCACGGACTCTCTAGAGCCGACAGGTCTCCTCGCAACA-----------C-AGACTTCTGATAATCATGA-TTCTGA--CA-------TATC-CGAATTTAAACCATTACACAGATATTCCAAC---GGACACACAGTTGCTTCTGCGTGGTAGTAA--AAACTCATAACGCAGAAACGTGCAGGATAATATGGTACCCTTCTAGGCAA---TGCGAGCAGACAGATACCCTCTCTGTACTTAGCGTATCAGT--------------ATTGC-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ATGAATTGCACGAGGCA-GTTTAG-GTAACGA-AGACGGCCGGTACGAATGTTGGAGCAATAGACCTCTTGGCT--------TTGGTCTGGCGTTTCCAAGCGCTAGTGTA--------CCTGCCTTTGAGCTGTTCCGCC--------T----GAAGTCC------CC-GT-GAACTCCGAGTTTTATCTCAATGTCTTAG +>I3 +GCACCTTA-GATGCAATGACTGTTTGATTGAT---CTTGGTAAGCGTTC-TTGGAAAAATGTCAC--------GCGAA------------------------------------------------------------------TGCATCTTCCAACACCTTGAAAAGTA----GATAG-CTG-------------AGACGAATCACGGACTCTCTAGAGCCGACAGGTCTCCTCGCAACA-----------C-AGACTTCTGATAATCATGA-TTCTGA--CA-------TATC-CGAATTTAAACCATTACACAGATATTCCAAC---GGACACACAGTTGCTTCTGCGTGGTAGTAA--AAACTCATAACGCAGAAACGTGCAGGATAATATGGTACCCTTCTAGGCAA---TGCGAGCAGACAGATACCCTCTCTGTACTTAGCGTATCAGT--------------ATTGC-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ATGAATTGCACGAGGCA-GTTTAG-GTAACGA-AGACGGCCGGTACGAATGTTGGAGCAATAGACCTCTTGGCT--------TTGGTCTGGCGTTTCCAAGCGCTAGTGTA--------CCTGCCTTTGAGCTGTTCCGCC--------T----GAAGTCC------CC-GT-GAACTCCGAGTTTTATCTCAATGTCTTAG +>I1 +ACACCTTT-GATGCAATGACTGTTTGATTGAT---CTTGGTAAGCGTTG-TTGGAAAAATGTCAC--------GCGAA------------------------------------------------------------------TGCATCTTCCAAAACCTTGAAAAGTA----GATAG-CTG-------------AGACGAATCACGGCCTCTCAAGAGCCGACAGGTCTCCTCGCAACA-----------C-AGACTTCGGATAATCATGA-TTCTGA--CA-------TATC-CGACTTTCAACCATTACACGGATATTCCAAC---GGACACACAGTTGCTTCTGCGTGGAAGTAA--AAACTCATAACGCAAAAACGGGCATGATAATATGATACCCTTCTAGGCAA---TCCGAGCAGACAGATACCCTCTCTGTACTTAGCGTATCAGT--------------ATTGC-------------------------------------------------------------------------------------------------------TGACGA---C---------ATGAATTGCACGAGGCA-GTTTAG-GTAACGA-AGACGGCCGGTACGAATGTTGGAGCAATAGACCTCTTGGCT--------TTGGTCTGGCGTTTCCAAGCGCTAGTGTA--------CCTTCCTTTGAGCTGTTCCGCC--------T----GCATTCC------CC-GT-GAACTCCGAGTTTTATCTCAAGGTCTTAG +>I2 +GCAACTTT-GATGCAATGACTGTTTGATTGAT---CTTGGTAAGCGTTC-TTGGAATAATGCCAC--------GTGAA------------------------------------------------------------------TGCATCTTCCAACATCTTGATAAGTG----GATAG-CTG-------------AGACGAATCACGGACTCTCTAGAGCCGACAGGTCTCCTCGCAACA-----------C-AGACTTCTGGTAATCATGA-TTCGGA--CA-------TATA-CGAATTTAAACCATTCCACAGATGTTCCAAC---GGACACGCAGTTGCTTCTGCGTGGTGATAA--AGACTCATAACACAGAACCGTGCAGGATAAAATGGTACCCTTTTAGGCAA---TGCGCGCAGACAGATACACTCTCTGTACTTAGCGTATCAGT--------------ATTGC-------------------------------------------------------------------------------------------G-----T-----GGACAA---C---------ATGAATTGCACGAGGCA-GTTTAG-GTAACCA-AGACGGCCGGTACGAATGTTGTAGGAATAGACCTCTTGGCT--------ATGCTCTGGCGTTTCCGAGCGCTAGTGTC--------CCTATCTTTGATCTGTTTCGAC--------T----GAAGTCC------CG-TT-GAACTCCGAGTTTTATATCACTGTCTTAG +>I8 +GGACCTTA-GCTGCGATGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGATCTTGCAACAACCTGATCAGTA----GATCG-GTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TAAGCATGGCTTCTGG--CA-------TATC-GGACTGTCTAACAATATACACTTACGACAAC---GGACACTCAGTTGCTTCTGCGTGGTAGGAGTGACACAAAGAACGCAGAAACCGGTATTGTAATTTGGTTTCCTTCTAGGCAA---GGCGCGCAGAGAGAAACCCTGTCTGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C---------ACAAACTGAACCAGGAA-GTTTAG-GTAACGA-AGACCGCCGGTACAAATGTTGCAGCAATAGACCTCTTGGAT--------CTGATCTAGCGTTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGTTTCGAC--------T----TAAGTAC------CC-TT-GAACTCAGAGTTTTATCTCAATGTCTTAG +>I7 +AG-CCTTC-ACTGGGCTGAAAGGAACGCGAGT---CTTAGTAAGAGGTC-TCGCTTAGACGGGAA--------GCACG------------------------------------------------------------------TGAATTCTCCATCTACCTCAACATTA----AACCA-GAG-------------CCATGA-TAATCAAGGATCAAATGCAGACAGGGTTATATAGCGGT-----------G-CGTACGTTT-TATCCATGGGACCTTAGACC-------ATTC-TGGCAGACTAAAGGTCTAGTCATACCAAAAT---GGACACTTATTCGGTACCTCGTGGTAGGGGGGATACAAACGACCTAGAGACCGGTATGGTGGTTTAATCTCCGTCTAACCAT---CGTGTATACCCTGCACCCCGTGCGGTGCTCCTCGATTGATG--------------TATGT-------------------------------------------------------------------------------------------------T-----GAGAGT---T---------ACGTGCTGAACCGCCTA-TTTATG-GTAACCA-GAAACGCCGTTAGAACTGGTACAACCCGAGTTCTCATAGAA--------CTGATATATCCTTTCTGGGCGCCCTAATA--------CCCAACACTGTGCTGATTTTAC--------T----TAGGTAC------CT-TG-GAACCCGGAGAAGAATCTTAATCTGTCAG +>I6 +AG-CCTTC-ACTGGGCTGAAAGGAACGCGAGT---CTTAGTAAGAGGTC-TCGCTTAGACGGGAA--------GCACG------------------------------------------------------------------TGAATTCTCCATCTACCTCAACATTA----AACCA-GAG-------------CCATGA-TAATCAAGCATCAAATGCAGACAGGATTA-ATAGCGGT-----------G-CGTACGTTT-TATCCATGGGACCTTAGACC-------ATTC-TGGCAGACTAAAGGTCTAGTCATACCAAAAT---GTACACTTATTCGGTACCTCGTGCTAGGGGGGATACAAACGACCTAGAGACCGGTATGGTGGTATAATCTCCGTCTAACCAT---CGTGTATACCCTGCACCCCGTGCGGTGCTCCTCTATTGATG--------------TATGT-------------------------------------------------------------------------------------------------T-----GAGAGT---T---------ACGTGCTGAACCGCCTA-TTTATG-GTAACCA-GAAACGCCGTTAGAACTGGTACAACCCGAGTTCTCATAGAA--------CTGATATATCCTTTCTGGGCGCCCTCATA--------CCCAACACTGTGCTGATTTTAC--------T----TAGGTAC------CT-TG-GAACCCGGAGAAGAATCTTAATCTGTCAG +>I22 +TGACCTTA-GCTGCGTTGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGAACTTGCAACAACCTGATCAGTA----GATTG-CTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTAGCAGCA-----------C-AGGATGCAG-TAAGCATGGCTTCTGG--CA-------TAGC-GGACTGTCTAACATTATACACTTACGACAAC---GGACACTTAGTTGCATCTGGGTGGTAGGAGAGACACAAAGAACGCAGCAACCGGTATTGCAATTTGGGTCCCTTCTATGCAA---GGCGCGCAGAGAGAAACCCTAGCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C-------TCACAAACTGAACCAGGAA-GTTTAG-GTAACGA-AGCCCGAAGGTACAACCG-------AATAGACCTCTTGGAT--------CTGATCTAGCGCTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGGTTCGAC--------T----TAAGTAC------CC-TT-GAACTCAGACTTTTATCTCAATGTCTTAG +>I14 +TGACCTTA-GCTGCGTTGACAGGATCACTGAT---CTTGGTAAGAGTTC-TAGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGAACTTGCAACAACCTGATCAGTA----GATTG-CTG-------------ACAAGAGTAACTGACTCTCTAGCGCCGACAGGTCTTCTAGCAGCA-----------C-AGGATGCAG-TAAGCATGGCTTCTGG--CA-------TAGC-GGACTGTCTAACATTATACACTTACGACAAC---GGACACTTAGTTGCATCTGGGTGGTAGGAGAGACACAAAGAACGCAGCAACCGGTATTGCAATTTGGGTCCCTTCTATGCAA---GGCGCGCAGAGAGAAACCCTAGCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGA---C-------TCACAAACTGGACCAGGAA-GTTTAG-GTAACGA-CGCCCGAAGGTACAACCG-------AATAGACCTCTTGGAT--------CTGATCTAGCGCTTCCAAGCGCCATTGTA--------CCCGCCTTTGAGCTGGTTCGAC--------T----TAAGTAC------CC-TT-GAACTCAGACTTTTATCTCAATGTCTTAG +>I12 +TGACCTTA-GCTGCGTTGGCAGTACCACTGGT---CATGGTAGGAGTTC-TGAGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGAACATGCAACAACCACATGAGAA----GATTG-ATA-------------ACAAGAGTAACAGACTCCCTAGGGCCGGCAGGTCTTCTAGCAGCA-----------C-AGGATGCAG-TAAGCATGGCTTCTGG--CA-------TAGC-GTACCGTCTAACATTAAACACTTACGACAAC---AGACACATAGTTGAGACTGGGTGGGAGGAGAGACACCTAGAACGCAGCAACCGGTATTGCAAATTGAGGCCCTTCTATGCAA---GGCGCGGAGAGACAAACCCAAGCGGTACTTATCGATCCAGA--------------GATGT-------------------------------------------------------------------------------------------G-----T-GCGCGGCCCG---C-------T---AAACTGGATCAGGAA-GTTTAG-GTAACGA-CGCCCGAAGGTACAACGG-------AATGGACCTCTTAGAT--------CTGATCTAGCGCTTCCAAGCGCCGTGGTA--------CCCGCGTTCGAGCTGGTTC-AC--------T----TAAGTAC------CC-TT-GACCTAAGACTTATATCTAAATGTCTTGG +>I10 +TGACCTTA-GCTTCGTTGGCAGTACCACTGGA---TATGGTCGTAGA-C-TGAGAAAAGTGGTCA--------TCGAG------------------------------------------------------------------TGGAACATGCAACAACCTCATGGGAA----AATTG-ATA-------------ACAAGAGTAATAGACTCCCTAGGGCCGGCAGGTCTTCTAGCAGCA-----------C-AGGATGCAG-TAAGCTGGGCTTCTGG--CT-------TAGC-GTACCGTATCACATTAAACACTTACGACGAC---AAACACATAGTCGTGACAGGGATGGAGGAGGGACACCTAGCACGCAGCAACCGGTATTGCAAATTGAGGCCCTTCTATGCAG---GGCACGTAGAGACAAACCCAAGCGGTACTTATGGATCCAGC--------------GATGT-------------------------------------------------------------------------------------------G-----T-GCGCGGGCCG---C-------T---TAACTG-ATCAGGAA-ATTGAG-CACACGA-CGCCCGAAGCTCCAACGG-------AGTTGACCTCTTAGAT--------CTCATCTAGCCCTACTAAGCTCCGTCGGA--------CCCGCGTTCTAGCTGGTCC-GC--------T----TCAGTAC------CC-TA-GACCTAAGGCTTACA-CTTAATGTCTTGG +>I11 +TGGCCTTA-GCTGCGTTCTCAGTACCACTGGG---CATGGTAGGCACTC-TGAAAAAAATCCCCA--------GCAAG------------------------------------------------------------------TGCAACATGCAACAACCACATGAGAA----AATTG-ATA-------------ACAAAAGTAACAGACTCCCTAGGGCCGTCAGGTCGTCTAGCAGCG-----------C-AGGAAGCAG-TAAGCATGCCTTCTGG--AA-------TAGC-GTACCGTCCAACATTAAACACTTACGACAAC---AGTCAAATAGTTGAGACGCGGTAGGACTAGAGACACCTGGAACGCAGCAACCGGTATTGCAAATTGTGGCCCTTCTATGAAT---AGCGCGGAGAGACAAACCCAAGCGGTACTTATCGATCCTGA--------------GATGA-------------------------------------------------------------------------------------------G-----T-GGGCGGTCCG---C-------T---AAACTGGATCAGGAA-GTTTAG-GTAACGA-CGCCCGAATGTACAACGG-------AATGGACCTCTTAGAT--------CTGATCTAGCGCTTCCAAGCGCCGGGGAA--------CCCGCGTTCGAGCTGGTCT-AC--------A----TAAGTAC------AC-TT-GACCTAAGACTAATATCAAAATGCCTTGG +>I13 +TGATACGC-GATGCGTTGACAGGAT-AGATAT---CGT-------ATTCCTAGGAACAATTGCCA--------GCACG------------------------------------------------------------------TGGAACTTTCAATATCCTTATGCGTA----GATTG-CTG-------------TCAAGGGTAGCTGA---TCTAGCGCCGCTAGGTCCTATAGCACCA-----------C-AGAATGCAG-TAAGTATGGCTTTTAG--CA-------TATC-GGACTAGCAAACCTTATACACTTGCGACTGA---GGACACTTCGTTGCAGCTGCGTGATCAGTGAGATATAAAGTGCCCAGCAAC--GTCATGCAATTTTGGCACATTCTATGCAA---GGCGCGCAGAGAGGAACGCTAGCGTT-CGTAACGATCGAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACAA---A-------CCACAAACTGGTGAAGGAA-GGTTAG-ATATCGA-AGACCGATGGTCCGAACT-------AATAGACCACTTGGTT--------ATGATGTAGCGCGTCAGGGCGCTATTGTA--------CCCAGCTTTCAGCTGCCTCGAC--------T----TCAGTGC------GC-TG-GAACTCACATTCTGATCTCAATGTCTCAC +>I21 +TAACCTTA-GGT-CCTTGATGGGATCACTGAT---CGTGGTAAACGTTA-TAGGACAAATGGCCA--------ACAAG------------------------------------------------------------------CGGAAGTAGCAGTCACCTCATCAGGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGGGCCGACAGGTCTTCTATCTGCA-----------A-AGTGTGTAG-AGAGAATGGCTTCTGA--CT-------TAGC-GGCCTGTCTAACATTATACTCTGACGACAAC---AGCTACTTAGTTGTATCTGGGTGGGAATAGAGACTCAAAGTACGCACCAGCAGGTATTGAGGTTCGGGTCCA-TCTGTGCAA---TGTTCGCTGAGAGAAAGACTAGCGGTACGTATCTATCCAGT--------------AATTT-------------------------------------------------------------------------------------------G-----G------GGCCA---T-------TCACAAAACGA-CCAGGAA-GTTTAG-GAACCGA-AGCCCGGAGGTATAGACG-------CATAGACCTCTCGGAT--------CTGTCTTAGAGCTGCCAACCCCCAATGTA--------CCCGCATTTGCGCTGGTTCGAT--------T----TTTGCAT------CC-TT-GAGCTCAGACTTTTATCTCGATGTCTTAA +>I15 +TAACCTTA-GGT-CCTTGATGGGATCACTGAT---CATGGTAAACGTTA-TGGGTCAAATGGCCA--------ACAAG------------------------------------------------------------------CGGAAGTAGCAGTCACCTCATCAGGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGGGCCGACAGGTCTTCTATCTGGA-----------A-AGTGTGTAG-AGAGAATGGCTTCTGA--CT-------TAGC-GGCCTGTCTAACATTATACTCTGACGACAAC---AGCTACTTAGTTGTATCTGGGTGGGAATAGAGACTCAAAGTACGCAGCAGCAGGTATTGAGGCTCGGGTCCA-TCTGTGCAA---TGTTCGCTGAGAGAAAGACTAGCGGTACGTATCTATCCAGT--------------AATTT-------------------------------------------------------------------------------------------G-----T------GGCCA---T-------TCACAAAACGA-CCAGGAA-GTTTGG-GAACCGC-AGCCCGGAGGTATAGACG-------CATAGACCTCTCGGAT--------CTGTCTTAGAGCTGCTAACCCCCAATGTA--------CCCGCATTTGCCCTCGTTCGCT--------T----TTTTCAT------CC-TT-GAGCTCAGACTTTTATCTCGATGTCTTAA +>I20 +TGACATTA-GGT-CCTTGATGGGATCACTGAT---CGTGGTAAACGTTA-TAGGACAAATGGCCA--------ACAAA------------------------------------------------------------------CGGAAGTAGCAGTCACTTCATCAGGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGGGCCGGCAGGTCTTCTATCTGCA-----------A-AGTGTGTTG-AGAGAATGGCTTCTGA--GT-------TAGC-GGCCTGTCTAACATTATACTCTGATGACAAC---AGCTACTTAGTTGTATCACGGTGGGAATAGAGACTCAAAGTACGGACCAGCAGGTATTGAGGTTCGGGTCCA-TCTGTGCAA---TGTTCGCTGAGAGAACGACTAGCGGTACGTATCTATCCAGT--------------AATTT-------------------------------------------------------------------------------------------G-----G------GGCCA---T-------TCAAAAAACGA-CCAGGAA-GTTTAG-GAACCGA-AAACCGGAGGTATAGACG-------CAGAGACCTCTCGGAT--------CTGTCTTAAAGCTCCCAACCCCCAATGTA--------CCCGCATTTGCGCTGGTTCGAT--------T----TTTGCAT------CC-TT-GAGCTCAGACTCTTATCTCGAAGTCTTAA +>I19 +TGACATTA-GGT-CCTTGATGGGATCACTGAT---CGTGGTAAACGTTA-TAGGACAAATGGCCA--------ACAAA------------------------------------------------------------------CGGAAGTAGCAGTCACTTAGTCAAGA----GATTG-CTT-------------ACAAGCGTAACCGACTCTACAGGGCCGGGAGGTCTTCTATCTTCA-----------A-CTTGTATTG-TGACAATGGCTTCTGA--GT-------TAGC-GGCCTGTCTAACATTATACACTGATGACAAT---AGGAACTTAGTTGTTTCACGGTGGGAATAGAGACTCAAAGTACGGAGCAGCAGGTATTGAGGTTCGGTTCCA-TCTGTTCAA---TATTGGGTGAGAGAACGGCTATCGGTACGCATGTATCCAGA--------------AATTT-------------------------------------------------------------------------------------------G-----G------GGACA---T-------TCACAAAACTA-TCCGGAG-GTTGAG-GAAACGA-AATCCGGAGGTATAGACG-------CAGAGACCTCTCGGAT--------CTGTCTTAAAGCTCCAAAGCGCCAATGTA--------CCCGCATTAGCATTGGTTCGAT--------G----TTCGCAT------AC-TC-GAGCTCTGACTCTTATATCGAAGTCTTAA +>I18 +ATACAGTA-GGT-CATTTATGGATTAACTGAT---AGTGGTCAACTTTA-TAGGACAAATGACCA--------ACAAG------------------------------------------------------------------CAACCGTAGCAGTCACTGAGTCAAGA----------CTT-------------ACAAGCGTAACAGTCTCTACAGGGGCTGCAAGTCTTCTATCTTCA-----------A-CTAGTATTG-TGATGATGGCTTCTGA--GT-------TCGC-CGCCTGTCTAGCATTATAGAGTGATGACTAA---AGTCACTTTGTTGTTTGTCGGTGACGATCTAAACTAAAAGGACCGGGCAGCAAGGATTGAGGTTCGGTTCCA-GCTT-TAAA---TATTTGGGGCGAGCACGCCTATCGGTGGACATGTCTCCAGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGCG---T-------ACACAAACCTA-TCCGTTG-GTTGAG-GATACTA-ATTCCGGAGGTATAGACG-------AAGAGATCTAGCGGTT--------CTGTTTTTTAGCTCCAAAGTGGTTATGTA--------GCCGCATTGGCATTCGTTAGGT--------G----TACGCAT------AC-TT-GAGCTCGGACTCTAATAACGTAGTGTTAA +>I17 +ATACAGTA-GGT-CATTTATGGATTAACTGAT---AGAGGTCAACTTTA-TAGGACAGATGACCA--------ACAAG------------------------------------------------------------------CACCCGTAGCAGTCACTGAGTCAGGA----------CTT-------------ACAAGCGATACAGTCTCGACAGGGGCTGCAACTCTTCTATCTTCA-----------A-CTAGTATTG-TGATGATGGCTTCTGA--GT-------TCGC-CGCCTGTCTACCATGATAGAGTGATGCCTAA---AGTCACTTTGTTGTTTGTCGGTGACGATCTAAACTAAAAGGACCGGGCAGCAAGGATTGAGGCTCGGTTCCA-GCTT-TTAA---TATTTGGGTCGAGCACGCCTATCGGCTGACATGTCTGCAGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGCG---T--------------CCTA-TCCGTTG-GTTGAG-GTTACTA-ATTCCGGAGGTATAGCCG-------AAGAGATCTAGCGGTT--------CTGGTTTTTAGCTCCAAAGTGGTTATGTA--------GCCGCATTGGCATTCGTTAGGT--------G----TACGCAT------AC-TT-GAACTCGGACTCTAATAACG-AGTGTTAG +>I16 +ATACAATT-TGA-CACTTATGAATTAACTGAT---AGAGGTCCACTTTA-TAGGACAGATGATCA--------ACAAG------------------------------------------------------------------CACCCGTAGCAGTCACTGAGTCAGGA----------CTT-------------ACAAGCGATACAGTCTCGACAGGGGCTGCAACGCTTCTATCTTCA-----------A-CTAGTATTG-TGATGATGGCTTCTGA--GT-------TCGC-CGCCTGACTACCATGATAGAGTGATGCCTAA---AGACACTTTGTTGTTTGTCGGTGACGATCTAAACTAAAAGGACCGGGCAGCAAGGATTGAGGCTCGGTTCCA-GCAA-TTAA---TATTTGGGTCGAGCACGCCTATCGGCTGACATGTCTGCGGA--------------AATAT-------------------------------------------------------------------------------------------C-----G------GGGCG---T--------------CCTA-TCCGTTG-TTTGAG-GTTACTA-ATTCCGGAGGTATAGCCG-------AAGAGATCTAGCCGTT--------ATCGTTTTTAGCTCCAAAGTGCTTATGTA--------GCCGCATTGGGATTCGTTAGGT--------G----TACGCAT------AC-TT-GAACTCGGACTCTAATAACG-AGTGTTAG +>I43 +GGCCCTTA-GCGGCCATGACAGGATCACTGAT---CTTGGTAAGAGCTC-TGGGAAAAATGGCCA--------GCAAG------------------------------------------------------------------TGGATCTTGCAACAACCTGACCAGTA----GATAG-CTG-------------ACGAGAGTAACTGACTCTCTTGCGCCGACAGGTCTTCTCGCAGCA-----------C-AGGATGCTG-TATGCATGGCTTCTAT--TA-------ACTC-GGCCTGCCTAACTATACACACTTACGACAAC---GGACACTTAGTTGCTTCTGCGTGGTAGGAGAGGCACTAAGACCGCAGCCACCGGAATGGTAATTCGGGTCCCCTCTATGCAA---GGCGCGCAGAGAGAAGCCCTATCGGTACGTATCGATCCAGT--------------TATGT-------------------------------------------------------------------------------------------G-----T-----GGACGG---C---------ACAAACTGAACCAGGCA-GTTTAG-GTAACGA-AGAACGCCGGTACAACTGTTGGAGCAATAGACCTCTTGGAT--------CTGAACTAGCGTTTCCCAGCGCCATTGTA--------CCCGCTTCTGAGCTGTTTCGAC--------T----TAAGTAC------GC-TT-AAAGTCAGACTTTTATCTCAGGGTCCTAG +>I42 +AGCCCTTT-GGGGCCATAACAGGATCACCGAT---CTTGTGAAGAGCTC-TGGGAAAACTGGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCAACTACCTGACCAGTA----GATAG-CCG-------------ACGAGAGTAACTGACTCTCTCGCGCCGGGAAGTCTTCGCCCAGCA-----------C-AGGCTTCTC-CATGCATGGCTTCTAT--TA-------ACTC-GGCCTGCCTAACCATGCACACTTACGAAAAT---CGCCACTTAGTTTCTTCTGCCTGGTGGGAGAGGCACTAAGACCGCAGCGACCGAAATGGTAACCCGGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGCTACGTATCTACTCAGT--------------TATGT-------------------------------------------------------------------------------------------------T-----GGACGG---C---------ACAAAGTGAACCAGGCA-GTTTAA-GAAACGA-AGAAAGCCGATACATCTGTTCGAACAATAGACCCCTTGGAG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CCCGCTTATGAGGTGT-TCGAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATTTATCTCAGGGTCCTAG +>I35 +AGCACTTT-GGGGTCGTAACAGGATCACCGAT---CTTGTGAAGGGCTG-TGGGAAAACTAGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCAACTAACTGACCAGTA----GAAAG-CCG-------------ACGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCCCAGCA-----------C-AGGCTTCTC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCTAACGATGCACAATTACGAAAAT---CGCCACTTAGTTTCTTCTGCCTGGTGGGAGAGGCACTAAGACCGCAGCGAACGAAATGGTAACCCGGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGCTACGTATCTACTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GGACGA---T---------ACAGAGTGAACCAGGCA-GTTTAA-GAAACGA-AGAAAGCGGATACATCTGTTTGAACAATAGACCCGTCGTAG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CCCGCTTATGAGGTGT-TAGAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATTTATCTCAGGGGCCTAG +>I34 +AGCACTTT-GGGGTCGTAACAGGATCACCGAT---CTTGTGAAGGGCTG-TGGGAGAACTAGCCA--------GCAAT------------------------------------------------------------------TGGATCGTGCAACTAACTGACCAGTA----GAAAG-CCG-------------ACGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCCCAGCG-----------C-AGGCTTCTC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCTAGCGATGCACAATTACGAAAAT---CGCCACTTAGTTTCTTCAGCCTGGTGGGAGGGGCACTAAGACCGCAGCGAACGAAATGGTAACCCGGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGCTACGTATCTACTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GGACGA---T---------ACAGAGTGAACCAGGCA-GTTTAA-GAAACGA-AGAAAGCGGATACATCTGTTTGAACAATAGACCCGTCGTAG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CCCGCTTATGAGGTGT-TAGAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATTTATCTCAGGGGCCTAG +>I31 +AGCACTTT-GGGGACGTAACAGGATCACCGAT---CTTCTGAAGCGCTG-TGGGAGTACTAGCCA--------GCAAT------------------------------------------------------------------TGGATCGTGCAACTAACTGACCATTA----GAAAG-CCG-------------GCGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCACAGCG-----------C-AAGCTTCAC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCAAGCGACGCACAATTACGAAAAT---CGCCACTTAGTTTCTTCAGCCTGGTGGGAGGGGCACTAAGACCGCAGCGAACGTAATGGTAACCCGGGGCCCGTCTATGCTA---GGCGTGCAGCGAGACATCCTATCGCTACGTATCTACTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCACGA---G---------ACAGAGTGAACCAGGCA-GTTTAA-GACACGA-AGAAAGCGGATACATATGCTTGAACAATAGACCCGTCGTAG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CACGCTTATGAGTTGT-TAGAT--------T----TAAGTCC------CC-TT-AAAGTCACATATGTATCTCAGGGGCCCAG +>I30 +AGCACTTT-GGGGACGTAACAGGATCACCGTT---CTTCTGAAGCGCTG-TGGGAGTACTATCCA--------GCAAT------------------------------------------------------------------TGGATCGTGCAACTAACTGACCATTA----GAAAG-CCG-------------ACGAGAGGAACTGACCCTCTCGCCCCGGGAAGTATTCGCACAGCG-----------C-AAGCTTCAC-CATGGATGGCTTCTAT--CC-------ACTC-GGCCTGCCAATCGACGCACAATTACGAAAAT---CGCCACTTAGTATCTTCAGCCTGGTGGGAGGGGCACTAAGACCGCAGCGAATGTAATGGTAACCCGGGGCCCGTCTATGCTA---GGCGTGCAGCGAGACATCCTATCGCTACGTATCTCCTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCACGA---G---------ACAGAGTGAACCAGGCA-GTTTAA-GACACGA-AGAAAGCGGATACATATGCTTGAACAATAGACCCGTCGTAG--------CTGAACTAGAGCTTAGCAGCGCCAACGTA--------CACGCTTATGAGTTGT-TAGAT--------T----TAAGTCC------CC-TT-AAAGTCACATATGTATCTCAGGGGCCCAG +>I25 +ACTACTTT-TGGGACGCAACAGGATCACCGTA---GTTCTTAAGCGTTG-TGGCCGTACTATCCA--------GCCAT------------------------------------------------------------------GAAATCATGCAACTAACTGCCAATTA----GGAAG-CCG-------------ACGAGAGGAATGGTCCGGCTCGCCCCGGAAAGTACTCGTACAGCG-----------C-TAGCGTCAC--ATGGATGGCTTAAAT--CC-------AGTT-GTCATGCCAATCGAAGCAAAGTTACGAAAAT---CGCCATGTAGTATCTTCAGCCTGGTGGTAGGGGCGCTAAGACGGGAGCGTAGATGATAGTAACCCTGGCCCGGCCTATGCGT---GCCAGCGAGCGGGACATGCTATCGCGACGTACATCGTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCAAGA---C---------ACACAGTCAACCAGGCA-GTTTAA-GTTCCGA-CGCAAACGGATACATATGCTCGAAGAATAGAGCCGTCGTCG--------CTGAAGTAGAGTTTAGCATCGCTATCGTA--------CACGCTTATCTCTTGT-TGCAA--------T----AAAGTAC------CA-CT-GAAGTCACATAAGTACACCAGGAGCCCAT +>I24 +ACTACTTA-TTTGAACCAACGGGGTCACTGTG---CTTCTGAATCGTTG-TGGCTCTCCTATCCA--------ACCAT-------------------------------------------------------------------CAATTATCCAACAAACTGCCCAGTG----CGAAG-CCG-------------ACGAGAGAAGTGGTCCGGCTCTCCCCATATAGTACTCGTACAGCG-----------C-TCGAATCAC---TGCATCGGTGAAAT--CA-------AGTT-GTCATGCCAATCACCGCAAAGATACGAAAAT---GGCCATGTAGTGTCTTCAGCCTGATGGTAGGGGCGCTAAGACGGGAGGGTAGGAGATCGTAACTCTCGGCCGGCCTATGCGT---GGAAGCCAGCGGACCATACTAAATGAACGTACTTAGTGAGT--------------TAAGG-------------------------------------------------------------------------------------------------------TCCTGT---C---------ACCAAGCCAACCAGGCA-GTTAAA-TTTACGA-CGCAAACGGATACATATGGTCGAACTTTAGAGCCTTCGTCG--------TTGTAGTAGACCTTAGCACCGCTATCGAA--------CACGCTTATCTCTGGT-TTCCA--------T----AAAGTAC------CT-CT-GAAGTCACATGAGTACACCAGCAGCCCAT +>I29 +AGCACTTT-TGGGACGTTACAGGATCACCGTA---CTTCTTAATCGCTG-TCTGTGTACTACCCA--------GGAAT-------------------------------------------------------------------GGATCGTGCAACTAACTGACCATAA----GAAAG-CCG-------------ACGAGAGGAACTGACCCTCTCGACCCGAGAAGAATGCGCACAGCG-----------C-GAGCTTCAC-CATGGACAGCTTCTAT--CC-------ACTC-GTCCTGCCAATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGTGGCCTGGTGGGAGGGGCACTAAGGCCGCAGCGAATGTAATGGTAAGCCGTGGCGCGTCTATGCTA---GGTGTGCAGCGAGACATGCTATGGCTTCGTATCCCCTCAGT--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCGA---G---------ACAGAGTGAACTAGGGA-GTTTGA-GACAGGA-AGAAAGCGGATACATATGCTTGAACTATCGACCCGTGGTAG--------CCGATCTAGTGCTTAGTAGCGCCAACGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGGGGCCCAG +>I28 +ACCAATTT-TGCGAGATTACAGGATCACCATA---CTTCTTCATCGCTG-TCTGTGTACTACCCA--------GGAAT-------------------------------------------------------------------CGATGGTCCAAGTAACTGACCATAA----GAACC-CCG-------------ACGAGAGGAACTGACCCTCTCGACCCGAGAAGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACAGCTTCTAT--CC-------ACTC-GTTCTGCCAATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGCGGCCTGGTGGGAGGGGCACTAAGGTCGTAGCGAATGCAACGGTAAGCCGTGGCGCGTCTATGCTA---GGTATTCAGCGAGACATGCTATGTCCTCGTATCTCCTCAGC--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCGA---A---------ACAGAGGGAACTAGGGA-CTTTGA-GACAGGA-AGAAAGCGGATACATATGCTTGAAGGATCGACCCCTGGTAG--------CCGATCTAGCGCTTACTAGCGCTACCGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGGGGCCCAG +>I26 +ACGAATTT-TCCGAGGTTACCGGATCACGATA---CGTCTTCATCGCTG-TCTGTGTACTTCCCT--------GGAAT-------------------------------------------------------------------CGATGGTCCAAATAAGAGACCATAA----GAGCC-CCG-------------ACCAGGGGAACTGACCCTCTCGAACCGAGATGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACTACTTCTAT--CC-------ACTC-GTTCTGCCCATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGCGGCCTGGTGGGAGGGGCACCACGGTCGTAGCGAATGCAACTGTAAGCCGTGCCGCGTCCAGGCCA---GGTATTCAGCGAGACACGCTACGTCCTCGTAACTCCTCAGC--------------TATGA-------------------------------------------------------------------------------------------------T-----GATCGA---A---------ATAGAGGGAACTAGGGA-CATTGC-GACAGGT-AGAAAGCGGATACGAATGCTAGAAGGATCGACCCCTGGAAGCAGACAAGACGACCTAGCGCTTACTAGCGCTACCGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGAGGCCCAG +>I27 +ACCAATTT-TGCGAGATTACAGGATCACCATA---CTTCTTCATCGCTG-TCTGTGTACTACCCA--------GGAAT-------------------------------------------------------------------CGATGGTCCAAGTAACTGACCATAA----GAACC-CCG-------------ACGAGAGGAACTGACCCTCTCGACCCGAGAAGAATGCGCACCGCG-----------C-TAGCTTCAC-CATGTACAGCTTCTAT--CC-------ACTC-GTTCTGCCAATCGACGCACAATTACGAAAAT---CACAACTTAGTATCTGCGGCCTGGTGGGAGGGGCATTAAGGTCGTAGCGAATGCAACGGTAAGCCGTGGCGCGTCTATGCTA---GGTATTCAGCGAGACATGCTATGTCCTCGTATCTCCTCAGC--------------TATGG-------------------------------------------------------------------------------------------------T-----GCTCGA---A---------ACAGAGGGAACTAGGGA-CTTTGA-GACAGGA-AGAAAGCGGATACATATGCTTGAAGGATCGACCCCTGGTAG--------CCGATCTAGCGCTTACTAGCGCTACCGTA--------CACGCT-ATTAGATGC-TAGAT--------T----TAAATTC------CC-TT-AAAGTCATATATGTATGTCAGGGGCCCAG +>I33 +AGCACCTT-GGGGTAGTAACAGGATCACCGAT---CTTTTGAAAGGCTG-TGTGAGAACTAGCAA--------ACAAG-------------------------------------------------------------------GGATCCTGCAGCTAACTGTCAAGTA----GTAAG-CCG-------------ACGACAGGAAATGAGCGTATCGCCCGGTGTATTAATCGCTCAGCT-----------C-AGGCTTCTC-CATGGATCGCCTCTAT--CC-------AATC-GTCTGGCCTAACGATCCACAAGTACGAACAT---CGCCACTTAGTTTTTTTAGGAAGGTGGGAGGGGCACTAAGACCGCAGCGAACGAAATGGTAACACGAGGTCTCTCGGTGCTC---TGGGTGCAGAGAGACATCGTATCGCTAGGTATCTACACAGT--------------TATAC-------------------------------------------------------------------------------------------------T-----GGACCA---T---------ACTGTGTGAACGACGCA-GTATAA-GAAACGA-AGCAAGCGGATACAGCTGTTTGTACAATCCAACAGTCCTAG--------CCGAACTAAAGTTCAGTAATGACATCGC---------CCCGCTTATGAGCTGC-TAGAC--------T----TAAGTCC------CCGTG-AAAGTCCGATATTTGTCTCAGGTGCCGAG +>I32 +AGCACCCT-GGGGTAGTACCAGGATCACCGAT---CTTTGGAAAGGCTG-TGTGAGAACTAGCAA--------ACAAG-------------------------------------------------------------------AGATCCTGCATCTAACTGTCAAGTA----GTAAG-CCG-------------ACGACAGCAAATGAGCGTATCGCCCGGTGTATTAATCGCTCAGCT-----------C-AGGCTTCTC-CATGGATCGCCTCTAT--CT-------AATC-GTCTGGCCTAACGATCCTCAAGTACGAACAT---CGCCACATAGTGTTTTTA-GAAGGTGGGAGGGGCACTAAGAGCGCAGCGAACGAAATGGTAACACGAGGTCTCTCGGTGCTC---TGGGTGCAGAGAGACATAGTATCGCTAGGTATCTACGCAGT--------------TATAC-------------------------------------------------------------------------------------------------T-----GGATCA---T---------ACTGTGTGAACGACGCA-GTATAA-GCATCGA-AGCAAGCGGATACAGCTGTTTGTACAATCCAACAGTCCTAG--------CCGAACTAAAGTTCAGTAATGACATCGC---------CCCGCTTACGAGCTGC-TAGAC--------T----CAAGTCC------CCGCG-CCAGTCCGATATTTGTCTCAGGTGCCGAG +>I41 +AGCCCTTT-GGGGCCATAACAGGATCACCGTT---CTTGTGAAGAGCTC-TGGGAAAACTGGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCAACTACCTGACCAGTA----GATAG-CCG-------------AGGAGAGTAACTGACTCTCTCGCGCCGGGAAGTCTTCGCCCAGCA-----------C-AGGCTTCTC-CATGCATGGCTTCGCT--TA-------ACTC-GGCCTGCCTAACCATGCACACTTACGAAAAT---CGCCACTTAGTTTCTTCTGCCTGGTGGGAGAGGCACTAAGATCGCAGCGACCGAAATGGTAACCCTGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGCTACGTATCTACTCTGT--------------TATGT-------------------------------------------------------------------------------------------------T-----GGACGG---C---------ACAAAGTGAACCAGGCA-GTTTAA-GAAACGA-AGAAAGCCGATACATCTGTTCGAACAATAGACCCCTTGGAG--------CTGAACTAGAGCTTAGCAGCGCCATCGTA--------CCCGCTTATGAGGTGT-TCGAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATTTATCTCAGGGGCCTAG +>I36 +AGCCGGTT-GGGGCCATAGCAGGATCACAGAT---CTTGCGATGATCTC-TGGGAAATCTGGCTA--------GCAAG------------------------------------------------------------------TGTATCGTGAAATTACCTGACCAGTA----GATGA-CCG-------------CGGATAGTAACTGACTCTCTCGCGCCGGGAAGTCTTCTCGCAGCA-----------C-AGGCGTCTC-CAGCCATAGCATCGCT--TA-------AGTC-GGCTTGCCTAACCCTGCACGCTCACGAAAAT---CGCTACTTATTTTCTTTTGCCCCGTAGGAGACGCAATAAGGTCGCAGCGACCCAAATGGTAACCCTGGGCCCCTCTATGCTA---GGCGTGCAGAGAGACATCCTATCGATACGATTCGACTCTGT--------------TATGA-------------------------------------------------------------------------------------------------T-----CGACGG---C---------AAAAAGTGCACCAGGCT-GTTTAA-GAAACGA-AGAATGCCGATACATCTGTTCGAACAAAAGTGTCCTTGGAG--------CTGAAATAGAGCTTAGCAGCCTCATCGAA--------ACCGCTTGTGAGGTGT-TCTAC--------T----TAAGTCC------CC-TT-AAAGTCAGATATCTATCTAAGGGGGCTTG +>I40 +AGCCCTTT-GGGCGCATAACAGGATCACGGTT---CCTGTGAAGCGATC-TGGCAAAGCTGGCCA--------GCAAG------------------------------------------------------------------TGGATCGTGCACCTGCCTGAACAGTG----GATAT-CCG-------------AGGAGAGTTGCTGACTCTTTCGCGCCGGGAAGTCCTTTCCAAGCA-----------C-AGGCTTCTC-TATGCATGGCATCGCT--TA-------TATT-TACCTGCCTAGCCATGCACCCTGACGAACAT---CGACACTTTTTTTATTCTGCCTTGTGAGTGAGGCTCTAAGATCGCAGCGACCGGACTTGTAACCCGGGGCCCCTCTATGATA---GACATGCAGAGAGACATACTATGGTTACGTATCCACAATGT--------------TATGT-------------------------------------------------------------------------------------------------T-----GTACGG---C---------ACAAAGTCAACCAGGCAAGTTTTA-GAAGCGA-AGAAGGCCGATACATCTGTTCCCACAACAGACCCCTGGGAG--------GTGAACTAGAGCTTAGCAGCGCCATCGTA--------CTCGCTTATGAGGTGT-TCGAC--------T----CAAGCCC------CC-TT-AAAGTCAGATATTTATCTCAGGGTTCTTG +>I39 +AGCCCTAT-TCGCGCCTCACCGTTTCACTGTA---CCTTTGACGGGATC-CGCCAAAGCTGCCCA--------GCAAG------------------------------------------------------------------TAGATCGTGCTTCTGCCATAACGGTG----GACAT-CCT-------------AGGAGAGTCGCTGAATCTTTCGCCCAAGGAGGTTCTTCCATACCA-----------A-TTTCTTCTA-TATGC-TGGGATCGCT--TA-------TCTA-GCCCTGCACTTTCCTGCGGCCTGAGGAAAAT---CGACACGTAATTTATTCTTCCTTGTGGCAGAGGGTTTAAGCTCGTAAGAACCGGCCTTTTACATAGGGGCCCCTCCATGTTG---GTCATGCTAAGAGTCAT-----GATTGCGTATACACAATGT--------------TGATT-------------------------------------------------------------------------------------------------T-----GTTCGG---A---------ACAATGTTCACATGGCAAGTTTCA-GAAGCGA-AGTAGGCAGATACATCGATCCCGACGATGGACCCCGGGAAG--------GTGAACTAGAGCATAGTAGCTCTTTGGTA--------CTCGCTTACGCCTTAT-TAGAC--------T----CAAGCCA------CC-TC-AACCTCCGAGAACTATCTAAGCGTTCTTG +>I38 +AGCCCTAT-TCGCGCCTCACCGTTTCACTGTA---CCTTTGACGGGATC-CGCCAAAGCTGCCCA--------GCAAG------------------------------------------------------------------TAG-TCGTGCTACTGCCATAACGGTG----GACAT-CCT-------------AGGAGAGTCGCTGAATCTTTCGCCCAAGGAGGTTCTTCCCGACCA-----------A-TTTCTTCTA-TATGC-TGGGATCGCT--TA-------TCTA-GGCCTGCAGTTCCCTGCGGCCTGAGGAAAAT---CGCCACGTAATTTATTCTTCCTTGTGGCAGAGGGTTTAAGCTCGTAAGAACCGGGCTTTTACATAGGGGCCCCTCCATGTTG---GTCATGATAAGAGTCAT-----GATTGCGTATACACAATGT--------------TGATT-------------------------------------------------------------------------------------------------T-----GTTCGG---A---------ACAATGTTCACATGGCAAGTTTCA-GAAGCGA-AGTAGGCAGATACATCGATCCCGACGATGGACCCCGGGAAG--------GTGAACTAGAGCATAGTAGCTCTTTGGTA--------CTCTCTTACGCCTTAT-TAGAC--------T----CAAGCCA------CC-TC-AATCTCCGAGAACTATCTAAGCGGTCTTG +>I37 +GGCCCTAT-TCGCGCCTCACCGTTTCACTGTA---CCTGTGACCGGATC-CGCCCAAGGTGGCCT--------ACAAG-------------------------------------------------------------------------------------------------------------------------------------------GCCCAAGGAGGTTCTTCCCGATCC-----------A-TCTCTCCTA-TATGC-TGGGATCGCT--TA-------TCTA-GGCCTGCAGTTCCCTGCGGCCTGAGAAAAAT---CGCCATGGAATTTAATCTACGTAGTGGCAGAGGGCTTAAGCTCGTAAGATCCGAGCTATTACATATGGGCCCCTATTTGTTG---GTCATGATAATACTCAT-----TATTTCGTATACACAATGT--------------TGATT---------------------------------------------------------------------------------------------GTACT-----GTTCGG---A---------ACAATGTTCACATGGAAAGTTGCA-GAAGCGA-AGTAGGAAAATACCTCGGTCCCGACGATGGAGCCCGGAAAG--------GTGAACTCGAGCATAGTAGCTGTTTGGTA--------CTCTCTTACGCCCTAT-TAGAC--------T----GAACCCA------CC-TC-TATCTCCGAGAACCATCTAAGCGGTCTTG +>I61 +TGTCCGTACACGTGGATGGCACGATCCCTGATATTTCTGAAAAAGGTGG-TAGCTGCCATGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGGCTAACAGACT---------------------TCGCAGAA-----------C-GGTTTATTG-TATGCATAGGGGCTGT--CC-------TGTA-CGACCGTCTAATATCATGCACCCACGACCAC---GAACAACTAGTTGCATCTCAGTGGTAAGTCTGACACAACAAATGCGGCACCAGGGGGTGGAAGGTGGGTGCCTTCCCCGAAA---GGTGCGCAA-GAGCAAGCCTGCCGTGACGTATACGTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATACAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATGGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCTA--------CCCGCCTATTGCCCGTATTGAT--------T----AAAATAA------CG-TC-GGAGCCAGACTCTTATATCGATGTACTAT +>I49 +TGTCCGTACACGTGGATGGCACGATCCCTGATATTTCTGAAAGAGGTGG-TAGCTGCCAAGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGTCTAACAGACT---------------------TTGCAGAA-----------C-GGTTTATTG-TATGCATAGGGGCTGT--CC-------TGTA-CGACCGTCTAATATCATGCACCCACGACCAC---GAACAATTAGTTGCATCTCAGTGGTAAGTATGGCACAACAAATGCGGCAGCAGGGAGTGGAAGGTGGGTGCTTTCCCCGAAA---GGTGCGCAA-GAGCAAGCCTGCCGTGACGTATACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATACAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATGGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCTA--------CCCGCCTATTGCCCGTATTGAT--------T----AAAATTA------CG-TC-GGAGCCAGACTCTTATATCGATGTACTAT +>I46 +TGTCCGTACACGTGAATGGCACGATCCCTGATATTTCTGAAAGAGGCGG-TAGCTGCCAAGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGTCTAACAGACT---------------------TGGCAGAA-----------C-GGTTTAGTG-TATGCATAGGGGCTGT--CC-------TGTA-CGACCGTCTAATATCATGCACCCACGACCAT---GAAGAATTAGTTGCATCTCAGTGGTTAGTATGGCACAACAAATGCGGCAGCAGGGAGTGGAAGGTGGGTGCTTTCGCCGAAA---GGTGCGCAA-GAGCAAGCCTGCCGTGATGTATACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATACAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATAGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGTTGCATTCTA--------CCCGCCTATTGCCCCTATTGAT--------T----AAAATTA------CG-TC-GGAGCCAGACTCGTATATCGATGTACTAT +>I45 +TGACCGTACACGTGAATGGCACGATCCCTGATATTTCTGAAAGAGGCGG-TAGCTGCCAAGACCG--------GCCTG------------------------------------------------------------------ATGAACATACGACAAGCTGTTCAGTTCACTGAAAG-CGA-------------ACTAGTCTAACAGACT---------------------TGGCAGAA-----------C-GGTTTAGTG-TATGCATAGAGGCTGT--CC-------TGTA-CGACCGTCTAATATCATGCACCCACGTCCAT---GAAGAATTAGTTGCATCTCCGTGGTTAGTATGGCACAACAGATGCGGCAGCAGGGAGTGGAAGGTGGGTGCTTTCCCCGAAA---GGTGCGCAA-GAGCAAGCCTGCCGTGATGTATACGTTCTAT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATAAAATCAC---------ACGAACTGAGCCAACAT-GTTTAG-TGAATGC-GAACCGCCCCTTCATAGCTTGCGGTAGGAGAACACTTCGAT--------CTAATCAACGT-GTTCGAGTTGCATTCTA--------CCCGCCTATTGCCCCTATTGAT--------T----AAAATTA------CG-TC-GGAGCCAGACTCGTGTATCGATGTACTAT +>I48 +GGTCCGTAATCGTGGAAGGCTCGCTCCCTTATATTTCTTGGATAGGTGG-AACCTGCCAAGACCA--------GTCTG------------------------------------------------------------------ATAAACATACGACAAACTGTTCAGGTCAATGAGAG-CGA-------------AATAATCTAACAGACG---------------------TTGCAGAA-----------C-GGCTACGTG-TTTGGATCGTGGTTGT--CG-------TGTA-CGACCGTATAATATCATGCACCCATGACCCC---GACGAATTAGTTGAAGCTCGCTGGTAACTATAGCACATCGATTGACCCGGCAGGGAGTGGAAGATCGGTGCTTTTACTGAGA---CGTACGCAA-GGGCCAGCCTTCCGTAACGTTTACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----AAACAATCAC---------ATGAATCGAGCCAACAT-GTTTAG-GGAGTGC-GAACCACCCCTTCCTGACTTACGGTAGGATAACACTTCGCT--------CTGGTCAGCTT-GATCGAGCAGGATTCCA--------CCCGCCGATTGCCGGTATTGAT--------T----AAAATTA------AG-AC-TGAGTCAGAGTGTGAAATCGATCTACTAT +>I47 +GGTCCGTAATCGTGGAAGGCTCGCTCCCTTATATTTCTTGGATAGGTGG-AACCTGCCAAGACCA--------GTCTG------------------------------------------------------------------ATAAACATACGACAAACTGTTCAGGTCAATGAGAG-CGA-------------AATAATCTAACAGACG---------------------TTGCAGAA-----------C-GGCTACGTG-ATTGGATCGTGGTTGT--CG-------TGTA-CGACCGTATAATATCATGCACCCATGACCCC---GACGAATTAGTTGAAGCTCGCTGGTAACTATAGCACATCGATTGACCCGGCAGGGAGTGGAACATCGGTGCTTTTACTGAGA---CGTACGCAA-GGGCCAGCCTTCCGTAACGTTTACGTTCTCT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----AAACAATCAC---------ATGAATCGAGCCAACAT-GTTTAG-GGAGTGC-GAACCACCCCTTCCTGACTTACGGTAGGATAACACTTCGCT--------CTGGTCAGCTT-GATCGAGCAGGATTCCA--------CCCGCCGATTGCCGGTATTGAT--------T----AAAATTA------AG-AC-TGAGTCAGAGTGTGAAATCGATCTACTAT +>I60 +TGGCCGTTCACGTGGATCGTATTATCCACGAGATTTCTAAAAAAGATGG-TAGCTGCCATGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCACTGAACG-CGA-------------ACAAGGCTAACAGACT---------------------TCGGAGAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGACTGTCTAATATC-TGCACCTTCGGCCAC---GCACAACTAGTTGCATCTCAGTGGTGATTCTGACACAAAAGATGCGGCATCAGGGGGTGGAAGGAGGCTGCTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACGTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTTAGACAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATGGCGAGCGGTAGGAGAACACTGCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCCA--------CCCGCCTACTGCCCGCATTGAT--------T----AGAAGAC------AT-TC-GGAGCAAGACTCTTATATCGATCGACTAT +>I51 +TGGCCGTTCACGTGGATCGTATTATCCACGAGATTTCTAAAAAAGATGG-TAGCTGCCATGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCACTGAACG-CGA-------------ACAAGGCTAACAGACC---------------------TCGGAGAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGAATGTCTAATATC-TGCACATTCGGGCAC---GCACAACTAGTTGCATCTCAGTGGTGATTCTGACACAAAAGATGCGGCATCAGGGGGTGGAAGGAGGCTGCTTTCGGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACGTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAT---------ACGAACATAGACAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATGGCGAGCGGTAGGAGAACACTGCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCCA--------CCCGCCTACTGCCCGCATTGAT--------T----AGAAGAC------AT-TC-GGAGCAAAACTCTTATATCGATCGACTAT +>I50 +TGGCCGTTCACGTGGATCGTAATCTCCACGAGATTTCTAAAAAAGATGT-TAGCTGCCCTAACCG--------GCCTT------------------------------------------------------------------ATGAACTTACAACAAGCTTATCAGTCCACTGAACG-CGA-------------ACAAGGCTAACAGGCC---------------------TCGGAGAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGAATGTCTAATATC-TGCACATTCGGGCAA---GCACAACGAGTTGCATCTCAGTGGTGATTCTCGCACAAAAGATGCGGCATTAGGGGGTGGAAGGAGGCTGCTTTCGGCGAAA---GGTGCGC-A-GAGCATGCCTGCCGTGACGTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAT---------ACGAACATAGACAACAT-GTTTAG-GGAGTGC-GAACCGCCCCTTCATGGCGAACGGTAGGAGATCACTGCGAT--------CTAATCAACGT-GTTCGAGCTGCATTCCA--------CCCGCCTACTGCACGCATTGAT--------T----AGAAGAC------AT-TC-GGAGCAAAACTCTTATATCGATCGACTAT +>I59 +TGGCCGTTCACGTGGATCGTATTATCCACGAGCTTTCTAAAAAAGATAG-TAGCTGCCACGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCACTGAACG-CGA-------------ACAAGCCTAACAGACT---------------------TCGGAGAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGACTGTCTAATATC-TGCACCTTCGGCCAC---GCACAACTAGTTGCATCTCAGTGGTGATTCTGACACAAAAGATGCGGCATCAGGGGGTGGAAGGAGGTTGCTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTTAGACAACAT-GTTTAG-GGAGCGC-GAACCGCCCCTTCATGTCGAGCGGTAGGAGAACACTGCAAT--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGAGCAAGACTCTTATATCGGTCGCCTAT +>I55 +TGGCCGTTCACGTTGATCGTATTATCCACGAGCTTTCTAAAAAAAATAG-TAGCTGCCACGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCACTGACCA-CGA-------------CCAGGCCTAACAGACT---------------------TCGGAGCA-----------C-GATTTATGG-TATACATGGGTGCTGA--CC-------TGCA-CGACTGTCTAATATC-TGCACCTTCGGCCAC---GCACAACTCGTTGCATCTCAGTGGTGATTCTGACACAAACGATGCGGCATCAGGCGGTGGAAGGAGGTTGCTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTTAGACAACAT-GTTTAG-TGCGCGC-GAACGGCCCCTTCCTGTCGAGCGTTAGGAGAACACTGCAAT--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-AC-GGAGCAAAACTCTTATATCGGTCGCCTAG +>I53 +TGGCCTGTCCCGCCCGTCGTAACATCCACCAACTTTCTACAAAAATTAG-TAGCTGTCTCCACCG--------GTCAG-------------------------------------------------------------------TGCACA-AAAACGAGGTGATCAGTAGAGTCACTA-CGA-------------CCAGGCCTAGCGGACT---------------------TCTGCGAA-----------C-GCTTTATTG-TAGACGGCGAGGCGGA--CC-------TGTA-CGACTGTCGAATGTC-GTCAGCTTGGGCCAC---GCACAACTCGTCGCTTGTCAGTTGTTGTTCTGACACATAAGATACAGAACCAGGCGGTAGGCGGAGGTGGGTGTCCGCGAAT---GGTGTGC-A-GAGCAAGCCTATCGTGACCCATTCCTTCTCT--------------AGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCCTTAAA---------ATGCACTAAGACCAATT-GTTTAG-TGCGCTC-GAACGGCCCCTTCCTGTCGAGCGAGAGGAGAACACAGTC-------------ATCAACGG-GTTCGTGCG---------------------------------------------------AGATGGC------AT-AC-AGAGCAAACCCTTATTATAGGTATCCTAG +>I52 +TGGGTTGACCTACCCGTCTTCATACGGTGCATCTTCCAACTAGTATTAG-TGGTTTACTCCC--G--------GGCAG--------------------------------------------------------------------CCACT-GGAAAGCTCAGATGAGTAGATGCACTA-TCA-------------CCAGGCCTAGCGGGCA---------------------TCTCCGACACC--------C-ACTTTGTTA-TGAACAACGAGGCTGG--GC-------GCTT-CGGCGGTAGAATAAC-GTCACCTTGGGCCTA---TCACCAATCGTGGCTTGGATTTTTCTCTGCTGACTCAAAATACACAGAAACGGGCGGTAGG---AGACCGATGTCTCGGTAT---ATTGAGC-A-GGGGCTCATTAACCTGAACCGTTCGTGCTGC--------------AATAG-------------------------------------------------------------------------------------------G-----C-----TTCACATAGA------------CACAAAGACGAACT-GATTAT-TGCCCGT-CAACTGCCCCTTCATGGCGAGAGCGAGGAGGATACATTG-------------ATCAT--C-GTGCGCTCG---------------------------------------------------AGGTCGC------AC-AGTACGTGCCACCCTTAAGTTAGGCTTCCTAG +>I54 +TGGCCGTTCACGTTGATCGTCTTATCCTCGAGCTTTCTAGAAAAGCTAG-TAGCTGCCACGAGCG--------GCCTG------------------------------------------------------------------ATGAACATCCAACAATCTCATCAGTCCACTGACTC-CGA-------------GTAGGGCTAACAGACT---------------------TCGGAGCA-----------C-GATTTATTG-TATACATGGGTGCTGA--CC-------TGCA-CGATTGCCTAAGATC-TGCACCCACGGCCAC---GCACAACTCGTTGCATCTCAGTGGTGATTCTGACACAAAGGATGCGGCATCAGGCGGTGGAAAGAGGTTACTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTACTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTCAAACGCCAT-GTTTAG-TGCGCGC-GAACGGCCCCTTGCCGTCGAGCGTTAGGAGAACACTGCAAT--------CTAATCGACGT-GTTCGATCT---------------------------------------------------AGAAGAC------AT-AC-GGAGCAATACCCTTAGATGGGTCGCCTAG +>I58 +TGGCCGTTCACGTGGATCGTATTATCCACGAGCTTTCTAAAAAAGATAG-TAGCTGCCACGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCACTGAACG-CGA-------------ACAAGCCTAACAGACT---------------------TCGGAGAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGACTGTCTAATATC-TGCACCTTCGGCCAC---GCACAACTAGTTGCATCTCAGTGGTGATTCTGACACAAAAGATGCGGCATCAGGGGGTGGAAGGAGGTTGCTTTCCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCAATCAC---------ACGAACTTAGACAACAT-GTTTAG-GGAGCGC-GAACCGCCCCTTCATGTCGAGCGGTAGGAGAACACTGCAAT--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGAGCAAGACTCTTATATCGGTCGCCTAT +>I56 +TGACCGTTCACGTGGATCGTAGGATTCACGAGCTTTCTACAAAAGAGAG-TAGCTGCCACGACCG--------GCCTG-------------------------------------------------------------------TGAACATACCACAAGCTGAGCAGTCCCCTGAAAG-CGA-------------ACAATCCTAACAGACT---------------------TCGGAGAA---GCCGGTAGC-GATTTATGG-AATGCATGGGGGCTGA--CC-------TGAA-CGACTGTCTGATATC-TTCACCTTCGGCCAC---GCACAACGAGTTCCATCTCAGTGGTGAGTTTGACACAAAAGATGCGGCATCAGGGGGTG---GGAGGTTACTTTCCGCGATA---GGTTAGC-A-GAGCAAGCCTGTCGTGTTCCATCCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------------T-----ATCCAATCAC---------ACGAACTTAGACACCAT-GTTTAG-TGAGCGC-GAACCGTCCCCCCAGGTCGATCGGTAGCAGACCACTACAAT--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGGGAAAGACTCTTATATCGGTCGCCTCT +>I57 +TGGCCGTTCACGTGGATCGTAGCATCCACGAGCTTTCTAAAAAAGATAG-TAGCTGCCAGGACCG--------GCCTG------------------------------------------------------------------ATGAACATACAACAAGCTGATCAGTCCCCTGAACG-CGA-------------ACAAGCCTAACAGACT---------------------TCGGAAAA-----------C-GATTTATGG-TATGCATGGGTGCTGA--CC-------TGAA-CGACTGTCTAATATC-TGCCCCTTCGGGCAC---GCACAACTAGTTGCATCTCAGTGGTGATTCTGACACAAAGGATGCGGCATCAGGGGGTGGAAGGAGGTTGCTTTTCGCGAAA---GGTGCGC-A-GAGCAAGCCTGCCGTGACCTATTCCTTCACT--------------TGTAG-------------------------------------------------------------------------------------------G-----T-----ATCCACTAAC---------ACGAACTTAGACAACAT-GTTTAG-GGACCGC-GTACCGCCCCTTCATGTCGAGCGGTAGGAGAACACTGCAAC--------CTAATCAACGT-GTTCGAGCT---------------------------------------------------AGAAGAC------AT-TC-GGAGCAAGACTCTTATATCGGTCGCCTAT diff --git a/phyne/data/tkf/params/tree.nwk b/phyne/data/tkf/params/tree.nwk new file mode 100644 index 00000000..077d35c1 --- /dev/null +++ b/phyne/data/tkf/params/tree.nwk @@ -0,0 +1 @@ +((((((T19:0.18253564794288973,T38:0.11367942155605133)I0:0.16425146544541422,(((T9:0.22583938333180698,T64:0.22097885669157666)I1:0.054183821944129976,(T26:0.3340918865732189,T45:0.020802408576665286)I2:0.09063885549294241)I3:0.0047460824223079555,T53:0.32209213609679777)I4:0.1839292401884232)I5:0.07314241512109995,((T58:0.08394683607877096,(T27:0.0491489104689503,T20:0.1362122314698043)I6:0.01650547329500701)I7:0.48267922529669793,T28:0.3044049327795188)I8:0.008367001972542086)I9:0.010196170811873594,((((T55:0.12428046929819993,T60:0.050085343119415704)I10:0.11820270164771268,(T49:0.01848690698958188,T63:0.1774423979803579)I11:0.10838789088802715)I12:0.1376004245586173,(T11:0.1642276343063575,T16:0.23394546782410375)I13:0.24476506567139106)I14:0.004237620391101004,((T17:0.44997272536442356,T7:0.01393710522409032)I15:0.026726619505505467,(T43:0.010296842961598943,(T12:0.013285992238172823,(T61:0.04854232125532925,(T15:0.08008380273115935,(T2:0.1884842957547165,T62:0.008318333310841837)I16:0.04285393226529239)I17:0.059542193031121564)I18:0.22895126415363304)I19:0.06569601149001386)I20:0.056123869230948274)I21:0.21302543812529912)I22:0.04004580133067057)I23:0.012445655095494477,(T59:0.05057180362001219,((((((T34:0.009011379268137306,(T35:0.15092936735000784,T18:0.14421500117374628)I24:0.2195911543125767)I25:0.19482360545148897,(T3:0.13497597675010692,((T47:0.004173032352879809,T8:0.03698040578185024)I26:0.10354138339568661,(T33:0.07284176715819353,T4:0.006926480297948681)I27:0.003565525870755711)I28:0.07343030357240612)I29:0.12164092793060345)I30:0.015668637774335442,T6:0.023394260001709032)I31:0.04663931056287233,((T57:0.07536371686048421,T46:0.015844617929883884)I32:0.034610637639086665,T23:0.06056260240953739)I33:0.18284343535272585)I34:0.008756284636401615,T30:0.020400223534816533)I35:0.09026894762141248,((T24:0.01938171467253015,T41:0.0709037427939508)I36:0.10666773456595566,(T32:0.008738851914929117,(T52:0.014363125985653624,(T14:0.004537605881616293,(T10:0.1617130077530088,T42:0.027224014132103343)I37:0.11072710045883706)I38:0.02496444232084729)I39:0.2935687797149779)I40:0.14578540314359661)I41:0.02541892067309106)I42:0.16648929748105917)I43:0.07465518152493947)I44:0.2633651531166342,((((T5:0.0015642318939174669,T21:0.010929646468554873)I45:0.030608235578181143,T39:0.1676954040171613)I46:0.03168557455717944,((T48:0.060082658555633786,T36:0.04719974765845245)I47:0.008110770252420133,T54:0.08040962503719787)I48:0.23073639547292413)I49:0.029391830019327318,(((T13:0.12467911675223041,T31:0.07625879488099768)I50:0.04684824250198325,T37:0.07923244627200617)I51:0.016833079514477097,(((T22:0.10325235201624201,(T40:0.23917218416733682,T44:0.04287095335037999)I52:0.44144184404921705)I53:0.31356753624075046,(T25:0.16354471117510178,T51:0.03041603974875198)I54:0.07837158237710552)I55:0.03727556659982548,((T50:0.014078595495307702,T29:0.2852890552055393)I56:0.11119058202350303,(T56:0.035043706834188504,T1:0.1719311210997305)I57:0.021675960026272287)I58:0.00011928968311421468)I59:0.021631984247609603)I60:0.13365688129311915)I61:0.16167251826699555)I62:0.1924133290575341; \ No newline at end of file diff --git a/phyne/src/likelihood/mod.rs b/phyne/src/likelihood/mod.rs index 1ccc76e3..cfbd7e8f 100644 --- a/phyne/src/likelihood/mod.rs +++ b/phyne/src/likelihood/mod.rs @@ -41,6 +41,32 @@ pub trait ModelSearchCost { fn set_param(&mut self, param: usize, value: f64); /// Returns the valid range for a model parameter [min, max], inclusive. fn param_range(&self, param: usize) -> ParamRange; + + // These are not atomic parameters of the model but rather combinations of parameters like their + // product or ratio. They are useful during parameters optimisation for the case that the atomic + // parameters are not independent and the optimisation is more efficient in the space of the + // helper parameters. + fn helper_param_count(&self) -> usize { + 0 + } + fn helper_param(&self, _param: usize) -> f64 { + panic!( + "helper_param not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } + fn set_helper_param(&mut self, _param: usize, _value: f64) { + panic!( + "set_helper_param not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } + fn helper_param_range(&self, _param: usize) -> ParamRange { + panic!( + "helper_param_range not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } fn set_freqs(&mut self, freqs: FreqVector) -> Result<()>; fn empirical_freqs(&self) -> FreqVector; fn freqs(&self) -> &FreqVector; diff --git a/phyne/src/optimisers/model_optimiser.rs b/phyne/src/optimisers/model_optimiser.rs index 639cb91d..1d462d5a 100644 --- a/phyne/src/optimisers/model_optimiser.rs +++ b/phyne/src/optimisers/model_optimiser.rs @@ -6,7 +6,7 @@ use argmin::solver::brent::BrentOpt; use log::{debug, info, warn}; use crate::evolutionary_models::FrequencyOptimisation; -use crate::likelihood::ModelSearchCost; +use crate::likelihood::{ModelSearchCost, ParamRange}; use crate::optimisers::{ModelOptimisationResult, SingleValOptResult, StopCondition}; use crate::Result; @@ -102,19 +102,64 @@ impl ModelOptimiser { Ok(cost) } + fn num_params(&self) -> usize { + self.c.param_count() + self.c.helper_param_count() + } + + fn get_param_or_helper(&self, idx: usize) -> f64 { + if idx < self.c.param_count() { + self.c.param(idx) + } else { + self.c.helper_param(idx - self.c.param_count()) + } + } + + fn set_param_or_helper(&mut self, idx: usize, value: f64) { + if idx < self.c.param_count() { + self.c.set_param(idx, value); + } else { + self.c.set_helper_param(idx - self.c.param_count(), value); + } + } + + fn param_range_or_helper(&self, idx: usize) -> ParamRange { + if idx < self.c.param_count() { + self.c.param_range(idx) + } else { + self.c.helper_param_range(idx - self.c.param_count()) + } + } + fn single_optimisation_iteration(&mut self) -> Result { let mut curr_cost = self.c.cost(); - for param in 0..self.c.param_count() { - let start_value = self.c.param(param); + for param in 0..self.num_params() { + println!("Optimising parameter {param:?} with current cost {curr_cost}"); + let start_value = self.get_param_or_helper(param); debug!("Optimising parameter {param:?} from value {start_value} with cost {curr_cost}"); let param_opt = self.opt_parameter(param, start_value)?; if param_opt.final_cost < curr_cost { // Parameter will have been reset by the optimiser, set it back to start value - self.c.set_param(param, start_value); + // self.c.set_param(param, start_value); continue; } - self.c.set_param(param, param_opt.value); + self.set_param_or_helper(param, param_opt.value); + + // The helper might set the lambda and mu correctly when setting a helper param, + // but not the before_helper_lambda and before_helper_mu. + // Alternatively we could have two trait methods one for setting helpers in the + // optimiser and one for setting helpers. Could be also called set_helper_repeatedly + // (for the optimiser that stays correct when calling the idx for the setter), and set + // helper that is not intended for repeated calls. Dont know what good names would be + for set_param_idx in 0..self.c.param_count() { + self.c.set_param(set_param_idx, self.c.param(set_param_idx)); + } curr_cost = param_opt.final_cost; + debug_assert_eq!( + curr_cost, + self.c.cost(), + "Cost mismatch after optimisation of parameter {param:?}" + ); + debug!( "Optimised parameter {param:?} to value {} with cost {curr_cost}", param_opt.value @@ -134,7 +179,8 @@ impl ModelOptimiser { cost: RefCell::new(self.c.clone()), param, }; - let range = self.c.param_range(param); + //let range = self.c.param_range(param); + let range = self.param_range_or_helper(param); let min = range.0; let max = range.1.min(start_value * 100.0); debug_assert!(min <= start_value && start_value <= max); @@ -155,6 +201,16 @@ pub(crate) struct ParamOptimiser { pub(crate) param: usize, } +impl ParamOptimiser { + fn set_param_or_helper(&self, cost: &mut C, param: usize, value: f64) { + if param < cost.param_count() { + cost.set_param(param, value); + } else { + cost.set_helper_param(param - cost.param_count(), value); + } + } +} + impl CostFunction for ParamOptimiser { type Param = f64; type Output = f64; @@ -165,7 +221,8 @@ impl CostFunction for ParamOptimiser { } else { *value }; - self.cost.borrow_mut().set_param(self.param, value); + // self.cost.borrow_mut().set_param(self.param, value); + self.set_param_or_helper(&mut self.cost.borrow_mut(), self.param, value); Ok(-self.cost.borrow().cost()) } diff --git a/phyne/src/optimisers/model_optimiser_tests.rs b/phyne/src/optimisers/model_optimiser_tests.rs index 0bcb28b3..6f698010 100644 --- a/phyne/src/optimisers/model_optimiser_tests.rs +++ b/phyne/src/optimisers/model_optimiser_tests.rs @@ -602,7 +602,9 @@ fn tkf_model_opti_template(c: C) { .unwrap(); assert_eq!(intermediate_cost, o.initial_cost); assert_eq!(o.final_cost, o.cost.cost()); - assert!(o.initial_cost < o.final_cost); + // The following assertion is not guaranteed to be true as the empirical optimisation may not + // improve the likelihood. + // assert!(o.initial_cost < o.final_cost); assert_eq!(o.cost.freqs(), &o.cost.empirical_freqs()); for param in 0..c.param_count() { assert_ne!(c.param(param), o.cost.param(param)); @@ -640,3 +642,19 @@ fn tkf92_model_opti() { .unwrap(); tkf_model_opti_template(tkf92); } + +#[cfg_attr(feature = "ci_coverage", ignore)] +#[test] +fn tkf92_model_opti_params() { + let fldr = Path::new("./data/tkf/params/"); + let info = PIB::with_attrs(fldr.join("masa.fasta"), fldr.join("tree.nwk")) + .build_with_ancestors() + .unwrap(); + // let subst_model = SubstModel::::new(&[], &[2.0]); + let subst_model = SubstModel::::default(); + let tkf_model = TKF92IndelModel::new(&[0.01, 0.01025, 0.65]).unwrap(); + let tkf92 = TKFCostBuilder::new(subst_model, tkf_model, info) + .build() + .unwrap(); + tkf_model_opti_template(tkf92); +} diff --git a/phyne/src/tkf_model/mod.rs b/phyne/src/tkf_model/mod.rs index 49c9d9cc..62004184 100644 --- a/phyne/src/tkf_model/mod.rs +++ b/phyne/src/tkf_model/mod.rs @@ -90,6 +90,22 @@ impl ModelSearchCost for TKFCos fn freqs(&self) -> &FreqVector { self.subst_cost.freqs() } + + fn helper_param_count(&self) -> usize { + self.indel_cost.helper_param_count() + } + + fn helper_param(&self, param: usize) -> f64 { + self.indel_cost.helper_param(param) + } + + fn set_helper_param(&mut self, param: usize, value: f64) { + self.indel_cost.set_helper_param(param, value); + } + + fn helper_param_range(&self, param: usize) -> ParamRange { + self.indel_cost.helper_param_range(param) + } } impl TreeSearchCost for TKFCost { diff --git a/phyne/src/tkf_model/tkf92.rs b/phyne/src/tkf_model/tkf92.rs index df3f18ac..26b5e89d 100644 --- a/phyne/src/tkf_model/tkf92.rs +++ b/phyne/src/tkf_model/tkf92.rs @@ -27,6 +27,17 @@ pub(crate) enum TKF92Parameters { Invalid(usize), } +pub(crate) const TKF92_HELPER_PARAM_COUNT: usize = 2; + +#[derive(Debug, Eq, PartialEq, FromPrimitive, IntoPrimitive)] +#[repr(usize)] +pub(crate) enum TKF92HelperParameters { + Rate = 0, + Ratio = 1, + #[num_enum(catch_all)] + Invalid(usize), +} + #[derive(Clone, Debug, PartialEq)] pub struct TKF92IndelModel { params: Vec, @@ -34,6 +45,10 @@ pub struct TKF92IndelModel { ln_r: f64, /// precomputed ((1 - r)/r).ln() ln_one_minus_r_over_r: f64, + /// original lambda before ratio or rate was set + before_helper_lambda: f64, + /// original mu before ratio or rate was set + before_helper_mu: f64, } impl TKF92IndelModel { @@ -42,16 +57,17 @@ impl TKF92IndelModel { /// 0 < lambda < mu < infinity, 0 < r < 1. pub fn new(params: &[f64]) -> Result { valid_number_tkf_params("TKF92", "lambda, mu, r", params, TKF92_PARAM_COUNT)?; - valid_lambda_mu( - params[usize::from(TKF92Parameters::Lambda)], - params[usize::from(TKF92Parameters::Mu)], - )?; - valid_r(params[usize::from(TKF92Parameters::R)])?; + let lambda = params[usize::from(TKF92Parameters::Lambda)]; + let mu = params[usize::from(TKF92Parameters::Mu)]; + valid_lambda_mu(lambda, mu)?; let r = params[usize::from(TKF92Parameters::R)]; + valid_r(r)?; Ok(Self { params: params.to_vec(), ln_r: r.ln(), ln_one_minus_r_over_r: ((1.0 - r) / r).ln(), + before_helper_lambda: lambda, + before_helper_mu: mu, }) } @@ -67,6 +83,8 @@ impl Default for TKF92IndelModel { params: vec![DEFAULT_LAMBDA, DEFAULT_MU, r], ln_r: r.ln(), ln_one_minus_r_over_r: (-r).ln_1p() - r.ln(), + before_helper_lambda: DEFAULT_LAMBDA, + before_helper_mu: DEFAULT_MU, } } } @@ -92,10 +110,17 @@ impl TKFModel for TKF92IndelModel { self.ln_r = value.ln(); self.ln_one_minus_r_over_r = (-value).ln_1p() - value.ln(); } - _ => { + TKF92Parameters::Lambda => { self.params[idx] = value; + self.before_helper_lambda = value; } + TKF92Parameters::Mu => { + self.params[idx] = value; + self.before_helper_mu = value; + } + _ => panic!("Invalid parameter index for TKF model: {param:?}"), }; + println!("set param: lambda = {}, mu = {}", self.lambda(), self.mu(),); } fn param_range(&self, idx: usize) -> ParamRange { @@ -134,6 +159,53 @@ impl TKFModel for TKF92IndelModel { fn get_blocks(&self, msa: &AA) -> Vec { blocks_of_alignment(msa) } + + fn helper_param_count(&self) -> usize { + // with that instead of zero the optimisation converges much faster and achieves better logl + TKF92_HELPER_PARAM_COUNT + // 0 + } + + fn helper_param(&self, idx: usize) -> f64 { + let helper_param = TKF92HelperParameters::from_primitive(idx); + match helper_param { + TKF92HelperParameters::Rate => (self.lambda() + self.mu()) / 2.0, + TKF92HelperParameters::Ratio => self.lambda() / self.mu(), + _ => panic!("Invalid helper parameter index for TKF model: {helper_param:?}"), + } + } + + fn set_helper_param(&mut self, idx: usize, value: f64) { + let param = TKF92HelperParameters::from_primitive(idx); + match param { + TKF92HelperParameters::Rate => { + let lambda = self.before_helper_lambda; + let mu = self.before_helper_mu; + self.params[usize::from(TKF92Parameters::Lambda)] = + 2.0 * value / (mu / lambda + 1.0); + self.params[usize::from(TKF92Parameters::Mu)] = 2.0 * value / (lambda / mu + 1.0); + } + TKF92HelperParameters::Ratio => { + let lambda = self.before_helper_lambda; + let mu = self.before_helper_mu; + let x = (lambda + mu) / 2.0; + let y = (x * (1.0 - value)) / (1.0 + value); + self.params[usize::from(TKF92Parameters::Lambda)] = x - y; + self.params[usize::from(TKF92Parameters::Mu)] = x + y; + } + _ => panic!("Invalid helper parameter index for TKF model: {param:?}"), + } + println!("set helper: lambda = {}, mu = {}", self.lambda(), self.mu()); + } + + fn helper_param_range(&self, _idx: usize) -> ParamRange { + let param = TKF92HelperParameters::from_primitive(_idx); + match param { + TKF92HelperParameters::Rate => (f64::EPSILON, f64::MAX), + TKF92HelperParameters::Ratio => (f64::EPSILON, 1.0 - f64::EPSILON), + _ => panic!("Invalid helper parameter index for TKF model: {param:?}"), + } + } } impl Display for TKF92IndelModel { @@ -297,6 +369,8 @@ mod private_tests { params: vec![0.5, 1.0, 0.3], ln_r: 0.0, // cache filled with dummy, since it is not needed here ln_one_minus_r_over_r: 0.0, // cache filled with dummy since it is not needed here + before_helper_lambda: 0.0, // cache filled with dummy since it is not needed here + before_helper_mu: 0.0, // cache filled with dummy since it is not }; // Use an invalid index model.param_range(3); @@ -308,6 +382,8 @@ mod private_tests { params: vec![1.1, 2.0, 0.3], ln_r: 0.0, // cache filled with dummy since it is not printed ln_one_minus_r_over_r: 0.0, // cache filled with dummy since it is not printed + before_helper_lambda: 0.0, // cache filled with dummy since it is not needed here + before_helper_mu: 0.0, // cache filled with dummy since it is not }; let fmt = format!("{}", tkf_indel_model); @@ -321,6 +397,8 @@ mod private_tests { params: vec![1.0, 2.0, 0.3], ln_r: 0.0, // dummy ln_one_minus_r_over_r: 0.0, // dummy + before_helper_lambda: 0.0, // cache filled with dummy since it is not needed here + before_helper_mu: 0.0, // cache filled with dummy since it is not }; let new_lambda = 1.1; model.set_param(usize::from(TKF92Parameters::Lambda), new_lambda); diff --git a/phyne/src/tkf_model/tkf_indel.rs b/phyne/src/tkf_model/tkf_indel.rs index 2a7503ed..4f391e8a 100644 --- a/phyne/src/tkf_model/tkf_indel.rs +++ b/phyne/src/tkf_model/tkf_indel.rs @@ -109,6 +109,28 @@ pub trait TKFModel: Clone + Display { fn params(&self) -> &[f64]; fn set_param(&mut self, idx: usize, value: f64); fn param_range(&self, idx: usize) -> ParamRange; + fn helper_param_count(&self) -> usize { + 0 + } + + fn helper_param(&self, _idx: usize) -> f64 { + panic!( + "helper_param not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } + fn set_helper_param(&mut self, _idx: usize, _value: f64) { + panic!( + "set_helper_param not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } + fn helper_param_range(&self, _idx: usize) -> ParamRange { + panic!( + "helper_param_range not implemented for this model. \ + Make sure to implement this method if helper_param_count() > 0" + ); + } /// Returns the factor corresponding to an insertion event at the root. fn ln_insertion_factor_at_root(&self) -> f64; /// Returns the factor corresponding to an insertion event at a non-root node. @@ -536,6 +558,24 @@ impl ModelSearchCost for TKFIndelCost &FreqVector { &DUMMY_FREQS } + + fn helper_param_count(&self) -> usize { + self.model.helper_param_count() + } + + fn helper_param(&self, param: usize) -> f64 { + self.model.helper_param(param) + } + + fn set_helper_param(&mut self, param: usize, value: f64) { + self.model.set_helper_param(param, value); + self.model_info.borrow_mut().valid.clear(); + self.model_info.borrow_mut().valid_for_reestimation.clear(); + } + + fn helper_param_range(&self, param: usize) -> ParamRange { + self.model.helper_param_range(param) + } } impl TreeSearchCost for TKFIndelCost {