2f032b790bd3d0354f1d61de5dcadb2e469c1817 braney Tue Mar 10 17:42:42 2026 -0700 Add unit tests for mafFrag tool Add test suite with 5 tests covering basic extraction, reverse strand, custom outName, small region, and larger multi-block region using hg38 multiz100way. Replace broken old test target that referenced a missing expected.maf file. Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com> diff --git src/hg/ratStuff/mafFrag/tests/expected/larger.maf src/hg/ratStuff/mafFrag/tests/expected/larger.maf new file mode 100644 index 00000000000..dfa1fe173c6 --- /dev/null +++ src/hg/ratStuff/mafFrag/tests/expected/larger.maf @@ -0,0 +1,29 @@ +##maf version=1 scoring=zero +a score=0.000000 +s hg38.chr1 100000 500 + 248956422 actaagcacacagagaataatgtctagaatctgagtgccatgtt--atcaaa-ttgtactgagactcttgcagtcacacaggc-tgacatg-taagcatcgccatgcctagtacagactctccctgca-gatgaaattatatgggatgctaaattataatgagaacaatgtttggtgagcca-aaactacaacaaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaattcactctcccttatcataagg------a---taatcttgccta-aagggacagtaatattAAAGAC-ACTAGGAATAACCTCTGTACTTTGGACAGT----AGACCTGCATAGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttgccattt-aaggtgattctcca-tgtttctgctcta-aCTGTGCTTCACAA-TACTCAAAACACTA----------AATCAGGATGTTTCCTGGAGTTCAGGGAGCTGTCCGTGTTACTGAGCAGTTC +s panTro4 0 384 + 384 ..........................................................................................................................cctgca-gatgaaattatatgggatgccaaattataatcagaacaatgtttggtgagcca-aaactacaacaaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaattcagtctcccttatcataagg------a---taatcttgccta-aagggacagtaatattAAAGAC-ACTAGGAATAACCTCTGTACTTTGGACAGT----AGACCTGCATAGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttgccgttttaaggtgattctcca-tgtttctgctcta-aCTGTGCTTCACAA-TACTCAAGACACTG----------AATCAGGATGTTTCCTGGAGTGCAGGGAGCTGTCCGTGTTACTGAGCAGTTC +s ponAbe2 0 499 + 499 actaaacacacagagaataatgtatagaatctgagtgccatatt--atcaaa-ttatactgagactcttgcagtcacacgggc-tgatatg-taagcatcgtcatgcctagtacagactctctctgca-gatgaaattatatgggttgctaaattataatcagaacaatgtttggtgagcca-aaactacaacaaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaatgcagtctcccttaccataagg------a---taa--ttgccta-aagggacagtaatattaaagac-actaggaataacctctctactttggacagt----AGACCTGCATAGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttaccattttaaggtgattctcca-tgtttctgctcta-aCTGTGATTCACAA-TACTCAAGACACCA----------AATCAGGATGTTCCCTGGAGTGCAGGGAGCTGTCCGTGGTCCTGAGCAGCTC +s nomLeu3 0 501 + 501 actaagcacacagagaataatgtacagaatctgagtgccatgtt--atcaaa-ttatactgagacgcttgcagtcacacaggc-tgatatg-taagcatcgtcatgcctagtacagactctccctgca-gatgaaattatatgggatgctaaattataatcagaacaatgtttggtgagcca-gaactacaacaaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaatgcagtctcccttatcataagg------a---taatcttgccta-aagggacagtaatattaaagac-actaggaataacctctctactttggacagt----AGACCTGCATAGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttaccattttaaggtgattctcca-tgtttctgctcta-aCTGTGCTTCACAA-TACTTAAAACACTG----------AATCAGGATGTTCCCTGGAGTTCAGGGAGCTGTCCGTGGTCTTGAGCAGCTC +s rheMac3 0 496 + 496 actaagcacacagagaataatgtatggaatctgagtgccatgtt--atcaaa-ttatactgagactcttgcagtcacacgggc-tgatatg-caaccattgtcacgcctagtacagactctccctgca-gatgaaattatatggcatgctaaattacaatcagaacaatgtttggtgagcca-aaactacaatgaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaatgcagtctcccttatcaaaaga------a---taatcttgccta-aagggatagtaat-----aGAC-ACTAGGAATGACTTCTGTACTTTGGACAGT----AGACCTGCAGGGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttaccattttaaggtgattctcca-tgtttctgctcta-aCTGTGCTTTACAA-TACTCAAAACACTG----------AATCAGGATGTTCCCTGGAGGTCAGGGAGCTGTCCGTGGTCCTGAGCAGCTC +s macFas5 0 496 + 496 actaagcacacagagaataatgtatggaatctgagtgccatgtt--atcaaa-ttatactgagactcttgcagtcacacgggc-tcatatg-taagcattgtcatgcctagtacagactctccctgca-gatgaaattatatggcatgctaaattacaatcagaacaatgtttggtgagcca-aaactacaatgaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcttaatgcagtctcccttatcaaaaga------a---taatcttgccta-aagggatagtaat-----aGAC-ACTAGGAATGACTTCTGTACTTTGGACAGT----AGACCTGCAGGGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatttaccattttaaggtgattctcca-tgtttctgctcta-aCTGTGCTTTACAA-TACTCAAAACACTG----------AATCAGGATGTTCCCTGGAGATCAGGGAGCTGTCCGTGGTCCTGAGCAGCTC +s chlSab2 0 488 + 488 actaagcacacagagaataatgtatggaatctgagtgtcatgtt--atcaaa-ttatactgagactcttgcagtcacacgggc-tgatata-taagcattgtcatgcctagtacagactctccttgca-gatgaaattatatggcatgctaaattacaatcagaacaatgtttggtgagc---------caatgaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagtttaatgcagtctcccttatcaaaaga------a---taatcttgccta-aagggatagtaat-----aGAC-ACTAGGAATGACTTCTGTACTTTGGACAGT----AGACCTGCAGGGCCCa-ttaggcctcaatgaagtcttatgcaagaccagaagccaatgtaccattttaaggtgattctcca-tgtttctgctcta-aCTGTGCTTCACAA-TACTCAAAACACTG----------AATCAGGATGTTCCCTGGAGGTCAGGGAGCTGTCCATGGTCCTGAGCAGCTC +s saiBol1 0 499 + 499 actaagcacacagagaataatgtgtagaatctgggtgtcatgtt--atcaaa-atacctggagattcttatagtcacatgggc-tgatatg-taagcattgtcatgcctagtacattctctctctgca-gatgaaattatacgggatgctaaattataatcagaacaatgtatggtgagcaagaaagtacaacaaggg-------aagctaattggatgaatttat-aaaaatatgcctcagccaaaatagcctaatgcaatctccctcatcataagg------g---taatcttgccta-aagggacagtaatgttgaagac-actaggaaaaacttctctactttggacaAC----AGACCTGCATAGCCTA-TTAAGCCTCAATTCAGTCTTATGT-ATATCAGAAGCCAATTTACCATTTGAAGGTGATTCTCCA-TGTTTCTGCTCTA-AGTGTGCTTCACAG-TACTCGAAACA-.G----------AATCACGATGTTCCCTGGAGTTCAGGGAGCTGTCCGTGGTCCTGAGCAGCTC +s otoGar3 0 374 + 374 actaaacaaacagaaaataatacatagaaccaagttgacaagtt--atcaaa-ttatactgagacttctgctgccacacatgc-taaaatg-ttagcattatcatgcgtagtacatcctctccctaca-ggtggaactatata-------acatcataatcagaataatgtctgatgatcaaggaaccaaaacaatgg-------aaattaattggataaatct-t-agaaatatacctcagccacaatagcttaaggcagtcttcctcatcagaagg------g---taatcttgcccc-atgagacagtaatattgaagAC-ATTAGGAATAACTTCTCTACTTTAGATAGT-----------------CCA-TTAGCCTTTAATTAAGTCCTATACAAGACCAAAAGCCggctcggcgccta-....................................................................................................................... +s tupChi1 0 495 + 495 ACTAAGCTTACCAAGAAAAATGTATGGAAACAGGGTGTCAAGTT--CTCAAA-TTATACAGAGACTTTTGCAATCACATGGGC-TGAAATG-TTAGTAT---CATACC---CACATTTTCTCCCTGCA-GATGGAATTATATGAGATGCTGAATTATAATAAGAATCATGTTTCATGAGCACGGATACACAAC------------AAGTTATTTGAATAAATTTGT-AGACATATGCTTCAGCCAAAATAGCATAAGGCAGTGTCCCTTGT---AAGG------A---TAATCTTGCCCA-ATGGGACAGTAATACTGaagac-acaagaaaaaaatcctcaattTTGGACAGTTC--TGACTTGCATAGTCCAGTTAGCCTTCAGTTCAGTCTTATACAAAACC-TAAGCCAGTTTACCATTTGAATGTGTTTTTCCA-GATCTGTGCTCTA-A--GTAACTCACAG-GGCTGGAGACCCTGA-AGTAGGGTAACCAAGATCTTCCCTAGAGTG-ACAGAGCTGTCCGTGGTGCTGAACCAA-A +s speTri2 0 450 + 450 actaagcatacagaaaataatatatagaa--------ctagTCA---------------TGGACTGGTTGCAGTCATATGAGC-TGAAAT--TTAATACTGTTatgtctagtaaactctcc-tctgca-gatggaatcatatgggattatg-attataataaga--------ttgtgggtaagtaaccacaactttgg-------aaataaagtagataaatttgt-agaagtatgcccctgtctaaatagcatacggcagtctccctcatcataaggataagga---taatcttgctca-atgggacagtaatattgagca--------agctactcctctgctttagacagttc------ttgcctggCTCA-TTGAGCT--CATAAAGTTTCATAGAAAATCAA--ATTAATTCATCATTTTAAGGTGTTTCACTATTATCTCTGTTTTA-ACTCTGGTTCAAAT-GACTCAAGACATCA----------AAATAG---------TGG------GGGAGTTGTCCATGGTCCTGAACCATTA +s octDeg1 0 372 + 372 attaagcccacagaaaataatgtatagca--------ccaagta--agcaaa-tcatactggatcctttgcagcaccatgggt-tgaaata-tttgcattctcatgcctagtatattctct-cctgcgtgatgaaattacatgaggtgttgaattataaccaaa-taatgcctggtgagcgaggaa--aaaacaatag-------aagttaactggacaatttttt-agaaatataccacagccaaaacaacataaagcagcccctctcctcataggga-----a---t-atcttaacca-atggagcagtaatattAAAGAC-AAAAAGAATAACTCATTTACTTCGGACAGTTCCTGGATTTGTGTGGTCCA-TTAACCTGCAATTAAGTCATGTGCAAAACCAAAA........................................................................................................................................ +s orcOrc1 0 510 + 510 attaagcacacagaaaataatgtatagaacccagatgccaagtt--accaaa-ttatactgaaatgtttgcattgacatgggc-tgaaatg-ttagcattgtcatacctaatatatattctctctaca-gctagaattatatgggatgttgaattataatcagaataatgtttggtgagcaagaaaccacaacaggggaagttaaaagttaattggatgcattttt-agaaacatacctcagccaaaacagcttaggacagtctcccttacaatagga------a---taatcttgccca-atgaggtagtaagattgaagac-actaggagtgacttccctactttgaacagtttctggatttgcacagtcca-gcagccttctgtgaAGTCTTACGCAAGACCAGAAATCAGCTTATCATTTTAAGGTGACTCTCCG-TATCTCTTATCTA-ACTCTgcccCACAA-TGTTTAAGATGCT----------TGAGGAAGACGGTTCCTAGAG---GGATAGCTGTCCGTGGTCCTGAACCACTA +s panHod1 0 503 + 503 gttaagcacac--aaaataatgtatagaaccaaggtgccaagtt--accaaa-ttatcttgcaatgtttgtagtgacatgggt-tgaaatgtttagcattgtcatgcctaatatattttctccctgca-gatggaattttatgggatgttgaatcataatcagaataatgtttggtgagcaagaaaccacaacaaagg-------aagttaattggatgcattttt-agaaacatgcctcagccaaaatagtgtgaggcagtctccctcacgaaggga------a---taatctagccCA-GTGAAGTGGTAAGATTGAATAC-ACTAGGAATGGCATCCCTACTTTGAATGGCTTCTGAATTTGCACAGTCCA-TTAGCCTTCTCTTAGGTTTCATGGAAGACTACATATCAACTTATCATTTTAAGGTGACTTTCTG-TATCTCTTCTCTA-ACTCTGACTTATAAAGGTTTAAGACCCT----------TGAGGAAGATGGGTCCTAGAG---GGATAGCTGTCCGTGGTTCTGAACCACTA +s bosTau8 0 505 + 505 actaagcacatggaaaataatgtatagaaccaaggtgccaagtt--accaaa-ttatgttgaaatgtttgtagtgacatgggt-tgaaatgtttagcattgtcatgcctaatatattttctccctgca-gatggaattttatgggatgttgaatcataatcagaataatgtttggtgagcaagaaaccacaacaaagg-------aagttaattggatgcattttt-agaagcatgcctcagccaaaatagcttgaggcagactccctcacgatggaa------a---taatctagccca-atgagGTGGTAAGATTGAATAC-ACTAGGAATGACATCCCTACTTCGAATGGTTTCTGAATTTGCACAGTCCA-TTAGCCTTCTCTTAGGTTTCATGGAAGACTAGAAATCAACTTATCATTTTAGGGTGACTCTCCG-TATCTCTTCTCTA-ACTCTGACTCACAAAGGTTTAAGACCCT----------TGAGGAGGATGGTTCCTAGAG---GGATAGCTGTCCATGGTTCTGAACCACTA +s oviAri3 0 503 + 503 gttaagcacac--aaaataatgtatagaaccaaggtgccaagtt--accaaa-ttatgttgcaatgtttgtagtgacatgggt-tggagtgtttagcattgtcatgcctaatatattttctccctgca-gatggaattttatgggatgttgaatcataatcagaataatgtttggtgagcaagaaaccacaacaaagg-------aagttaattggatgcattttt-agaaacatgcctcagccaaaatagcttgaggcagtctccctcacgatggga------a---taatctagccCA-GTGAAGTGGTAAGATTGAATAC-ACTAGGAATGGCGACCCTACTTTGAATGGTTTCTGAATTTGCACAGTCCA-TTAGCCTTCTCTTAGGTTTCATGGAAGACTAGAAATCAACTTATCATTTTAAGGTGACTTTCTG-TATCTCTTCTCTA-ACTCTGACTCACAAGGGTTTAAGACCCT----------TGAGGAAGATGGTTCCTAGAG---GGATAGCTGTCCGTGGTTCTGAACCACTA +s capHir1 0 503 + 503 gttaagcacac--aaaataatgtatagaaccaaggtgccaagtt--accaaa-ttatgttgcaatgtttgtagtgacatgggt-tgaaatgtttagcattgtcatgcctaatatattttctccctgca-gatggaattttatgggatgttgaatcataatcagaataatgtttggtgagcaagaaaccacaacaaagg-------aagttaattggatgcattttt-agaaacatgcctcagccaaaatagcttgaggcagtctccctcacgatggaa------a---taatctagccCA-GTGAAGTGGTAAGATTGAATAC-ACTAGGAATGGCGTCCCTACTTTGAATGGTTTCTGAATTTGCACAGTCCA-TTAGCCTTCTCTTAGGTTTCATGGAAGACTAGAAATCAACTTATCATTTTAAGGTGACTTTCTG-TATCTCTTCTCTA-ACTCTGACTCACAAAGGTTTAAGACCCT----------TGAGGAAGATGGTTCCTAGAG---GGATAGCTGTCCATGGTTCTGAACCACTA +s equCab2 0 479 + 479 tctaagcacacagagaataatgtatagaaccaaggggccaagtc--atcaaa-ttatgctgagacttctgcattgacatggac-tgaaatt-ttagcgttgtcatgcctagtacattctctccctgca-gatggaattatatgggatgctgaactgaaatcagaataatgtttggtgagcaaggaagcataacaagag-------aagttaattggataattTTTT-AGAAACATATCTCAGACAGAACAGCTTAAGTCAGTCTCTCTCACAACACAG------G---TAATCTTGCCCA------------------AAGAT-GCTAAGAATGACCTTTCTACTTTGGACAGTTTCTGGACTTGCACAGTCCA-TTAGCCTTCTGTTAAGTCTTATACAAGACCAGAAATCAATTTCTCATTTTAAAACGACTCTC------CCATTTTCTA-ACTCTAATTCACAA-T-TTCAAGACGCT----------CAACTCAGACGGTGCCTGCTG----GGGAGTTGTCCATGGTCCTGAACCACTA +s cerSim1 0 487 + 487 actaagcacacagagaataatgtatagaaccatggtgccaagtt--atcaaa-ttatgctgagacttttgcagtgatatggacttgaaatg-ttagcattgtcatgcctcatacattctctccctgca-gatggaattaaacgggatgctgaactggaatcagaatgatgtttgatgagcaaggaaccacatcgaggg-------aagctaactggataaattttt-agaaacatgtctcagccaaaacagcttaaggtagtctccctcacaacagag------a---cagtcttgccCA------------------ACAACGGCCAAGGATGACCTCCCTACTTCGGAAAGTTTCTCGACTTGCACAGTCCA-TTAGTCTTCTATTAAGTCTTATACAAGACCAGAAACTGATAAATCATTTTAAGATGACTCTCCA-TGTCTCCTCTCTA-ACTCTGATTCGCAA-TGTTCAGGACACT----------CAACTAAGATACTGCCTAGGG----GGGAGCTGTCCGCAGTCCTGAACCACTG +s felCat8 0 501 + 501 accaagaacatagagaacattgtatagaactaagatgccaagtc--atcaaa-ttatactgacatttttgcagt-acgtggac-tgaaatg-ttagcattgtcaagcctagcatattctctccctcta-gatggaattatatgaaatgctgaattataatcagaatagtgttttgtgagcaaagaatcacaacaaggg-------aagttaattggataaattttt-agaaacatgcctcagccaaaacagcttaagacagtcaccctcagagtaggg------a---taatcttgctcaagggggatagtaatattgaagAT-ACTAGAAATGACTCCTCAACTTTGGATAGTTTCTGGGCTTGAACAGTCTA-TTGGCCTTCTTATAAGTCTTAGGGAAGATCAGAAATCAGCTTGTCATCTTAAGGGGA--TGACA-TATCTGTTCTCCA-ACTCTGATTCACAC-TGTTCAAGACATT----------TAACCAAGAC-TTTTCTGGAGT--GGGGAACTGTCCGTGGTTCTGAACAACTA +s canFam3 0 494 + 494 actaagaacatagagaatattgtataaaactaagatgccaagtt--atcaaa-ttatactgcgattcttgcagt-acaggggc-tgaaatg-ttagcattgtcaagcctagtacattctctccctgta-catggagttatatgagttgctgaattataatcagaatagtgtttggtgagcaaagaatcacaacaaggg-------aagtt-attggataaATTTT--AGAAATGTGCCTTGGTCAAAATAGCTTAAGGCGGTCC------CAGGGGGT------ATATTAATCTTGCCCAGGGGGTATATTAA--TTGAAGAT-ACTAGGAATGGCCCTTCGACTTTGGACAGTTTCTGGGCTGGAACAGTCTA-TTGGCGGTTTTATAAGTCTTATGCAAGACCAGAAATCAATGTGTCATCCTAAGGGGACTCTCCA-TATCTGTTCTTTA-ACTCTGATT--CAG-TGTTCAAGACAT.----------TAACCAAAATGTTTCCTGGAGC--AGGGAGCTGTCCCTGGTTCTGAACAACTA +s musFur1 0 499 + 499 actaagaacatagagaatattgtatagaactaagatgccaagtt--atcaaa-ttatactgagatttttgcagt-acatgggc-tgaaatg-ttagcattgtcgagcctagtacagtctctccc-ata-gctggaat--tatgagatgctgaattataatcagaatcgtgtttggtgagcaaagaatcacaacaaggg-------aagtcaattggataaattttt-agaaacgtgcctcagccaaaacagcttaagacaatccgcctcacagtaggg------a---taatcttgccCAGAGGGGATAGCAAGGTTGAAGAT-ACTAGGAATGTCCCCCGGACTTTGGACAGTTTCTGGGCTTGAACAGTCTG-ATGGCCTTTTACAAAGTCTTATGCAAGACCAGAAATCAATTTGTCATCCTAAGAGGACCCTCCA-TATCTGTTCTCTA-ACTCTGGTTCACAG-GGTTCAACATTC.----------AAGCAAGAGC-TGTCCTGGAGT--AGGGAGCTGTCCCTGGTTCTGAACAACCA +s odoRosDiv1 0 495 + 495 actaagaacatagagaatattgtatagaactaagatgccaagtt--atcaaa-ttatactgagatttttgcagt-acatgggc-tgaaatg-ttagcattgtcaagcctagcacattctttccctgta-gatggaattatatgaaatgctgaattataatcagaatagtgtttggtgagcaaagaatc---acaaggg-------aagttaattggataaatgtttaagaaacgtgcctcagccaaaacagctcaatacagtacccctcacaataggg------a---caatcttgcccaggggGGATATTAATGTTGAAGAT-ACTGAGAATGG-CCCTGGACTTTGGACAGTTTCTGGGCTTGAACAGTCTA-TTGGCCTTCTTATAAATCTTTTGCAAGACCAGAAATCAATTTGTCATCCTAAGGGAACTCTTCA-TATCTATTTTCTA-ATTCTGATTCACAG-TGTTCAAGATTC.----------AAACTGG-----TTCCTGGAGC--AGGGAGCTGTCCCTGGTTCTGAACAACTA +s pteAle1 0 498 + 498 actaagcacacagagaataaagtatagaaccaaggtgccaagtt--gtcaaatttatactgagatctttgcagtgacatgagc-tgaaatg-ttagcactgtcatgcctactacattttctccctgca-gatggaattatatgggacgctgaattataattagaataatg-ttggtgaacaaggaaccac-acaatgg-------aagttaattggataaaatttt-agaaacacatttcagcccaaacagagtaagacagcctcacttacaataggg------a---caatcttgctta-ccgggacAGACAGATTGA--AC-ACTAGGAATGACCGCTCTACATTGGACAGTTTCT-GACTTGCACAGTCCA-TTAGGCCTTTATTAAGTCTTATACAAGACCAGGAATCAACTTCTCATTTTAAGGAGACTCTCCA-TATCTCTTCTCTA-ACTCTACTTAACCA-TGTTCAAGATACT----------TGACCAAGATAGTTCCTGCAG---GGAGCACTCTCCA-GGTCCTGGACCACTA +s loxAfr3 0 510 + 510 actaagctcacagataggaatgtatagaactaggttaccagggtacatcaaa--tatagtgagacttctgtagtcacatgggc-tgaaatg-tcagcattgtcatacctagcacattctctccctgaa-gatgaaattatatggggtgctgaattataataagaatagtatttggtaagtaaggaaccacaacaagga-------acgttaattggataatttttt-agaaaaatgtcttagtcaaaatagcttaaggctttctccttcaccataagg------a---taatcttgcc-a-acgaaacagtaatattgaagAT-GCTAGGAATGATCCCTCTAC-TAGGACTAGTCCTGGGCTTGCACAGCCCA-GAAGCCTTCAGTTAAGCCTTATACAAGACCAGAAGTCAGTTTACCATCTTAAATTCACTCTCTA-TATCTCTGTTCTA----CTATTGCACAG-TGCTCAAAACACT-GAAGCAGAATAACCACGATGTTTCCTGGACCA-TGTGTGCTGTCCATGGTGCTGAACCACTG +s triMan1 0 508 + 508 actaagctcacagatataaatgtatggaactagggtaccagggt--atcaaa--tatagtaagacttctgtagtcacatggac-tgaaatg-tcagcattgtcatacctagcacattctacccctaaa-gatg-aattatatggaatgttgaattataatcagaataatgtttgataagtaaggaacgacaacaaggg-------a-gtcaattggataaattttt-agaaacataccttagtcaaaatagcttaaggcagtctccttcaccataagg------a---taatcttgcc-a-gtggaacagtaatattgtagac-actaggaatgacccctctac-ttggacaattcctggGCTTG--GAGTCCA-TTAGCCTTCAGTTAAGCCTTATACAAGATCAGAAGTCAGTTTATTATCTTAAGTTGACTCTCCA-TATCTCTCTTCTAGTCTGTGTTTCACAA-TACTCAAAACACT-GAAACAGCATAACCAAGATGTTTCCTGGAGCA-TGGGTGCTGTCCGTGGTTCTGAACCACTA +