add UCSC's axt to MAF converter (#6199)

* add UCSC's axt to maf converter * allow selecting fasta indexes from data tables match folder name and tool id * add data table test, cleanup * Apply suggestions from code review --------- Co-authored-by: Björn Grüning <[email protected]>
galaxyproject · Aug 27, 2024 · 329ba03 · 329ba03
1 parent 8f49573
commit 329ba03
Show file tree

Hide file tree

Showing 10 changed files with 244 additions and 0 deletions.
diff --git a/tools/ucsc_tools/ucsc_axttomaf/.shed.yml b/tools/ucsc_tools/ucsc_axttomaf/.shed.yml
@@ -0,0 +1,9 @@
+name: ucsc_axttomaf
+owner: iuc
+description: Convert dataset from axt to MAF format.
+long_description: This is a tool from the UCSC suite to convert dataset from axt to MAF format.
+homepage_url: http://hgdownload.cse.ucsc.edu/admin/exe/
+remote_repository_url: https://github.com/galaxyproject/tools-iuc/tree/master/tools/ucsc-tools/axttomaf
+type: unrestricted
+categories:
+- Convert Formats
diff --git a/tools/ucsc_tools/ucsc_axttomaf/axttomaf.xml b/tools/ucsc_tools/ucsc_axttomaf/axttomaf.xml
@@ -0,0 +1,105 @@
+<tool id="ucsc_axtomaf" name="axtToMaf" version="@TOOL_VERSION@+galaxy0" profile="21.05" license="MIT">
+    <description>
+        Convert dataset from axt to MAF format
+    </description>
+    <macros>
+        <token name="@TOOL_VERSION@">377</token>
+    </macros>
+    <xrefs>
+        <xref type="bio.tools">UCSC_Genome_Browser_Utilities</xref>
+    </xrefs>
+    <requirements>
+        <requirement type="package" version="@TOOL_VERSION@">ucsc-axttomaf</requirement>
+    </requirements>
+    <version_command><![CDATA[ echo "@TOOL_VERSION@" ]]></version_command>
+    <command detect_errors="exit_code"><![CDATA[
+    axtToMaf
+        '$in_axt'
+        #if str( $target_reference_index_source.target_reference_index_source_selector ) == "history":
+            '${target_reference_index_source.in_tar_ref_index}'
+        #else:
+            '$target_reference_index_source.tar_ref_index_file.fields.path'
+        #end if
+        #if str( $query_reference_index_source.query_reference_index_source_selector ) == "history":
+            '${query_reference_index_source.in_que_ref_index}'
+        #else:
+            '$query_reference_index_source.que_ref_index_file.fields.path'
+        #end if
+        '$out'
+    ]]></command>
+    <inputs>
+        <param name="in_axt" format="axt" type="data" label="axt file" />
+
+        <conditional name="target_reference_index_source">
+            <param name="target_reference_index_source_selector" type="select" label="Choose the source for the target reference genome index">
+                <option value="cached">Built-in</option>
+                <option value="history">My history</option>
+            </param>
+            <when value="cached">
+                <param name="tar_ref_index_file" type="select" label="Using reference genome index">
+                  <options from_data_table="fasta_indexes" />
+                </param>
+            </when>
+            <when value="history">
+                <param name="in_tar_ref_index" format="tabular" type="data" label="Dataset with chrom sizes for target sequence. Typically a '.fai' file." />
+            </when>
+        </conditional>
+
+        <conditional name="query_reference_index_source">
+            <param name="query_reference_index_source_selector" type="select" label="Choose the source for the target reference genome index">
+                <option value="cached">Built-in</option>
+                <option value="history">My history</option>
+            </param>
+            <when value="cached">
+                <param name="que_ref_index_file" type="select" label="Using reference genome index">
+                  <options from_data_table="fasta_indexes" />
+                </param>
+            </when>
+            <when value="history">
+                <param name="in_que_ref_index" format="tabular" type="data" label="Dataset with chrom sizes for query sequence. Typically a '.fai' file." />
+            </when>
+        </conditional>
+
+    </inputs>
+    <outputs>
+        <data name="out" format="maf" />
+    </outputs>
+    <tests>
+        <test>
+            <param name="in_axt" value="1.axt"/>
+            <conditional name="target_reference_index_source">
+                <param name="target_reference_index_source_selector" value="history" />
+                <param name="in_tar_ref_index" value="hg38.fasta.fai"/>
+            </conditional>
+            <conditional name="query_reference_index_source">
+                <param name="query_reference_index_source_selector" value="history" />
+                <param name="in_que_ref_index" value="mm39.fasta.fai"/>
+            </conditional>
+            <output name="out" file="1.maf"/>
+        </test>
+        <test>
+            <param name="in_axt" value="1.axt"/>
+            <conditional name="target_reference_index_source">
+                <param name="target_reference_index_source_selector" value="cached" />
+                <param name="tar_ref_index_file" value="hg38"/>
+            </conditional>
+            <conditional name="query_reference_index_source">
+                <param name="query_reference_index_source_selector" value="cached" />
+                <param name="que_ref_index_file" value="mm39"/>
+            </conditional>
+            <output name="out" file="1.maf"/>
+        </test>
+    </tests>
+    <help><![CDATA[
+**What it does**
+
+`axtToMaf`_ is a tool to convert dataset from `axt format`_ to `MAF format`_.
+
+.. _axtToMaf: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/FOOTER.txt
+.. _axt format: https://genome.ucsc.edu/goldenPath/help/axt.html
+.. _MAF format: https://genome.ucsc.edu/FAQ/FAQformat.html#format5
+    ]]>    </help>
+    <citations>
+        <citation type="doi">10.1093/bib/bbs038</citation>
+    </citations>
+</tool>
diff --git a/tools/ucsc_tools/ucsc_axttomaf/test-data/1.axt b/tools/ucsc_tools/ucsc_axttomaf/test-data/1.axt
@@ -0,0 +1,19 @@
+0 hg38.chr1 52448 53340 mm39.chr2 111209343 111210259 + 33103
+GCCAATGGATAGAGTAAATAATTCTGCGGTATCTAAATTTGTATTGATTGGACTTTCAAGCTCTTGGGAGATGCATCTTTTTCTTTTTTGGTTCTTCTCTGTGTTCTACATGGGAATTATCCTGGAAAATCTCTTCATTGTGTTCACAGTAATTATTGACTCTCATTTA-AATTCCCCAGGTACTGCCTACTGGCCAACATTTATCTTCTTGATCTGGG----TCTTCTCCTACAGT-------------TCTGACTTTTTCACTAACTGCAGCATCATTTCTTTTCCAAGATGCATCATACAGATATTTTTCATTTGTGTCATGCGTAAAAATTGAGATGGTGCTGCTCATAACCATGGCATAGAGCAGGTACACTGCCAATCTGTAAGCCTCCCCATTACCTGACCACAATGAACCCCAAAATGTGTGTTTCCTTTGT---TGGAGGCATCCTGGATAGTCAGGATAATCCATGCTGTATCTCAGTTTGTTTTTGCCATAAACTTGCCTTTTTGTGGCCCTAATAGAGTAGGTAGTTTTCACTGTGATTTTCCTTATGTCATGAAACTTGCTTGTGTAGACACTTACAAACTAGAGGTTGTAGTCACTGCTAACAGTGGGCTTATATCCATAGCTACCTGTTTCTTATTAATAATATCCTATATTTTCATTTCGGTAACCGTCTAGAATCCTTCTTCAGGAGACTTATCTAAAGCATTTGT-GTCATGTTAGATCACATCACAGTAGGGATTTTGTTTTTTATGCCATGTATATTTCTGTATGTGTAGCCTTTGCCTAAAACAACACATGATTAATATTTGTTCATTGTTCCTTTTGCTATCACCCCTG--------TCTAGGATCTACACATTAAGAAACAAAGACATGAACGTCTCCATGGAAAGACTGGGAAAATGGA
+GTCAATGGATCAAATAAATGAAACTGTAGCTTTTGAATTTGTGTTGTTGGGACTCTCCAGTTCTTGGAAAAATACAATTTTTCTCATGTCCACCTTCTCTTTGCTCTATGTAAGCATCATCGTGGGAAACCTTTTCATTGTCTTTTTGGTAATTAATGACTCCCGTTTACAGTCTCCTATGTATTTTCTTTTGGCCAATCTGTCCCTCATTGATGTTGGACTTTCCTCTACCACAGTGCCCAAGATGATCTCAGATCTTCTGAAAGAGCACAAAGTAATATCTTTCCACAGTTGCATGACTCAGATCTGCTCCATCCACATTATG-GGAGGAGTGGAGATGGTGCTGCTCATAGCCATGGCGTTTGACAGGTACACGGCCA-TCTGTAAGCCTCTGCACTACATGAGCATCATGAGCCCTAGAATATGCATTTCATTTGTAATTGCTGGC----TGGGTCACTGGAGTGGTCCATGCCATGTCACAGTTCTCTTTTGTTGTGAACCTGCCTTTTTGTGGTCCTAACAAAGTGGACAGCTTTTACTGCGACTTCCCTCGAATCATACAACTGGCATGCACGGATAGAGACACATTTGAGTTTGTTGTTGCTGCCAACAGTGGCTTCATGACCTTGGGGACCTTCTTCTTGCTTCTCCTCTCCTATGTCTTCATTTTGGTCACAGTCTGGCAAAGGTCTTCAGGGGACTTGTCAAAGGCGCTTGTCACTCTGTCAGCTCACATCACTGTGGTTGTTCTGTTTTTTACTCCATGCATGTTTCTCTATGTGTGGCCTTTCCCCACATCATCGACTGACAAATACCTGTTCATTGTTGACTTTGCTGTCACCCCTGCTCTTAATCCTGTCATCTATACTTTAAGGAACAAAGATATGAAGGAAGCCATAAAAAAGCTGAGCAAGCAGA
+
+1 hg38.chr1 56555 58924 mm39.chr2 111303515 111305845 + 43395
+CACATTATTCTGCAGTTTTTATCTTGCTTACGTATTTTTGGAATGTTACTATTTTTCAAAAATTAATTTGGGATCAACCAACACTTCTTATTCTGCTGCTGTTCTAGAGAAAATCATTTTCCTCATTTCTGAACAAGAGAAAATGAAAT-----------------ACAGCTCTAAACAAATGCCACTGTAAACCAAGGTGGAGCCTTTGCACTTTCAGGCCACCATGATAACCT---GGAGATTAGATTTTTCTGTGTCTTTATATCAATAATAAAGCCAAGCTTCTCCAGGGGTATCCACTAGGCTTGTCTCAATGGCTCAATACAGGTCCTTTTGTGAATGATTACCTCACCCTCATGGAAACACACTCTT-GTTACAGAAACTCAGAATGATTCTATTTTTTCTTTTATATTTGTATATGTTTTTCCAATACCTCTGAAAAAACTGATCCAAAAAAAATACAAATTTTAATTGTAGCCAGTCAATTCAGGAAGGATAAAGGTCAAAAACTTTCAAAGAAACCTTCAGCCCCAACACACTAAACTTTGGGAGCACAGGTTGGCATCCAGAGGTAAACATTTGCTATAACTGATAACAGGAGAAGGATCCATTTATTCACCTGTTATCAATTACAGGCATTGTATTTAAAGATCAGATGTTTTATATTTATTTCTTCAAATTTCATTCATGGTGCCATAAGTGAAGGTATCTCTGTCCACCCTGAATATATTTTCACTCCCTCATCTCAGTCATTCCGAACAATTCACACACTAAGATTACCCATGCTAAATGGGGATTCTTTTTTACTAGCCAATGTAGTACCTCAAATCCTTCCTTCCCTCCCCCTATTTCATCAGCAGGCAATTCTTTTGATACTTTTGTCAAGGGGAAATTGTGTGACTCAGAGATCTAGTCCCCAAGAGAAACTAATAATGGGCTGGGTATTGTCTGTCTCAGCAGCATCAGTGGGTCCCTCTCCTGTGCAGCTAATTAGCTTCCTTTCCAATATGAAGAATCTTATATATAGCTTTGTCTTT-GGGGTATTAC----ATAAATGAAGATTAAGCTATCTGAATTTCTCCTTCTCCTAAAAATGCACATCCTATGACTGAAAAGACAGGTAAAAGAGATGCTTTTAATTACAAAACTTTCCCTGTCGTGGTTGCTTCTCTCTATCCTTCTAAACTCCCTTTCAATTTCTTCTCTTCTGTAACATATTTGTGCCCAAAATCTTCTGCTTTCTGAAATATTTT--ATCTTTTTCTTCCACACTATCTC----TTATTTTCCAATTTTAATCATTAAATTATATTATGTCTTATAAAACTAATCCCACATATAAACCCCTATGATA--ATTTCAGTTTG--TCCCTAGTATGAAG-TTCTTTAAAGATGTGTAGTTTTCTAACTTTCATGCTCTCCAATTAATTATAAACTTCATTTTCCACTCTGAAAAGGAGATGTCTGATCTCAGCTATTTCCATCCTATTTGAAAACCAGATTTAGTTTTAAACCAGAGGAAGGGAATCTCAAGTCTTTACCTCCCACAGTCTGGTGTGATTCTCTCTCTTTTGGTATTACCTTCCTCCACATTGGAACACTCCAGCCAATGCATAGGCTGAGAGGCTATCTC-AGATTCAGAAAGATTTGGCCTCATCCCAGGGGAGGGTACAGAGGAGCTGATGACTATGAATTC--TGAAATGGAACTGTTCCAG-------------GTTGAAGAAATAAGAAAGGGAATTGGGAAGAGCAAT-----GCCCAGTGAAAAAGAAGAAATAA-TATTTTAGGAAGT-GAATGC-TAATTTTATTTTAAACAAAATAAGAACTCAAGGAATAAGAGGGTTCTTCCAATAGGTTAGAGTGATCCTGTCAAACATATATGCTTCT-AGATTTTTTTAAAGACTGTTTCTACTAAGAAAGCATAGACCGCTATTGAGAAAGATCATTAAACTGGAATTTAGGAGGTCTGCCTTCTGATTCTGACTTCTTGAATGTATTGTTAGCCAT----TTAACCACACTGTGTTGTTTC--------TCA--TTCTACCTGTAGAATCTCAAAGTTCTTTCCCACTTCTATACAAAACTATAATTCTGAACATCCTTTTTGTTTAATATAAGTCTGCATTTCCTGTTTGAA--GATATGTGTCCC-AGACCCTAAATGACTGACAAATTTTAAATCTCC---AATAGGAAAGATGACAAACTCTATGGAAACTTGGCTTCTGAAGAACTCCTAGAAGCTTTCCAAAGTCATCAGTGTTTCCTAAGAAGGCAGAGAAATCAAACAC-----------------------------------ATGG----------------------------TCTTTTCCTCCAGACAAGCTCCTTTGGGTCATCAGGATTTCTTCAACAATA---AAATGTAATAATTCCAAATGTTTGTAACAGAATGGGTAGGACTTTCTTCACTTATTTAAATACTCCCTTTTT
+cacattgtattgta--ctctattatgttt-tattgtattatattatTGTTGTTTTTCAAAAATT----TTGGTTAAACTAGCATCCCTTCTTTTTCTGTTGACCTAAGAAGAAACCCTTTCCTTCTTT-TGAACAAGATAAAATGAAGTGCAGATTTCTAAACAACACAGCTTTAAAGAAGCACACTTATAAACAAAGGTAGAGCTTTCACTGTTTCA-------ATGA-AACCTCAGTGAGATCAGATT-----ATGACTCTACATAAATGACAAAG------------------ACCCACCGAGACTGGTACCATGTTTTCATTCATATGGTGTTA--AATGATGGTCTC---CTCGGGTATCCAGACCCTTCATTTAAGTAACTCAAAACAACT----------------ATCTATAT-----CTTCCAAAATTTCTAAACAAATTGACCCGATGACT-TACAAGTTTAACTTATAATAGATCAACTCA--AATGATAAAA--TAGTGATTTTCAGATATACCTATATCCTCTGG---------TTTTGTAGAGCATTCAGGCTTTCTTGGACAAACATTTGTCA-AAATGA-AACAATAAGAACTTTTATCTTTTTCCATACTAACAAGTGAAG----TGTATTTTCAGATAATTAGTTTTAGACT-ATTTTT---AATTTGATATATGAATCCA-GAGTAAAAATACAGCTATCCACTTTTGGAATATTTTTTATTCT-----------ATTTCAAGCAATTCATAAAAAAAGATTATAACTCCTCAA-GGAGGTTATTTTAAATT---TAATACATTACACAGAATATTTTTATTTTTTCCCTT-TTATATAACCAGACAATTCTTTTGATACTTCTCTCTTGGGGAAACCA-ACAGCTCAT-GATCCAATTCCCAATAGAAACTATTAATGAACTGGGTATTGGACACCTTCACAGTGCCAATGGGTCACTCTCCTGAGCAAATAATTGGCTTCTTTGCCACTATAAACAGTCATACATATATCTTTGTCTTCAGGGGCATTGCTTATAGAAATGAATGATAACCTACTTAAATTT-TCTTTTTCTTCAAAGCGTACCTGCTGTTACTTTAGAGA-AGGT--AAGGGGTGGTTTACATTTAAAACCTG--CCGATTATAGTTG----TTTCAATCTTGCTAA------CTCTAAT---CTCCCATCTGTAAAATAATTATGTCAGACAT-GTCACCTTTCAAGAGTATTTTTACCCTTTTAATTTTAAACTATCACAGAATCGTTTTGTGATTTTATTTA--AAAATAT----TGCTTTATATAATTAATCTCGTGCTCACAACCTTTTGATATTTTTTTAATTTGCTTCTTTGGTATGAAGATTTTTCAAGAACGTGTG--TTTCTGCTTTTTGTAGCCCTCAGACAATTGTGCAATTCATTTCCAACTATGCCAAGCAGATGTCTGGTCTCAGATATTTCCTGCATA-TTAAAGATCATATTACAATCTCAAT-----------AATCTCAAGCCTTGTCTGCCAATACTCTGGAGTAATTCTTACCTTT----TATT-------TCAATATTTGAATACTGTAAGCAGTGAATACACTGAGAGACTATTCCAAGACTTAGATATATT-----CTAGCCCACAGGATTCTAAAAAAGTGC-ATTAAATAGAAATTCTTTTTATTGCAACTGTTCAAGTAGTGGTCTAGATGTTTGAAAAA-AAGAAAAAAAATTGGCAAAGTCTATTAAGGGCACAGGGTAGGGAAAGGACTAGTTTTTTTACGAAGTAGGATGTGTAGTTTTATGCCAAACTGAGTATGAGTTCAAAGACTAACACAATCCTTTCAGTAAAGCAAAGTGATCATATCAAATATATAGTATTCTGAGATTTT----ACAACTGTTTCTCCTAATCATTTGTAGAGTAATACTGACGAA---CAGTACACTAGAGTTTAGGACACCTACCCTCTGA-TCTGTAGTTCTGTATGTTTTGCTCTTGATGCTGTTAA--ACACTAGGCTGTTTTATGCAGAATCAAGTGTTAAATGCGTTTTCTCACAGTTGTTCCCAGTTTTCATGTAAATCTACGATTCCATTAGTATGCTTTGTTTAACACATGTCTATGTATCTTGTTTGAAGTGATACTTATCCCTAGACCTTGAATAATAGACAAATTTTAAGTCACCACAACCAGGAAAACCATCAAACACTATGGGGACTTGGTTTCTGATGCACTCCTTGAAGCTCCATGCTGTTATCAACAGCTCT---GAAGGGAAAGGAATAAAGCACTATattttgcttttgttttgattttatattttgaGATGGAGCTGAGAAGTTGTTTTGTTGTATTTTTACTATTTTTACACATAAATTCTCTAGAAACATGTGGATTtaat-aataatatataaaaataataaCCCAAAATGCTTGTAA-------GATAAAAAGACCTTCATTTATTTAAGTTCTCATTCTTT
+
+2 hg38.chr1 60917 61125 mm39.chr2 111306441 111306644 + 8501
+TTAATATTTCCCACTTCATCATTACTAACCTCACTCACTCTGAACTTGATGAGAGCACCTGAGCATTAATTTTTCTTATAATTATTTAATGATTACCAGAATTCGTTCAGTATGGCCAGCTCTGGTCAAAGTGAGGCAGGCAAGATGCTTTGTCAACTGCCTGGATGGAATGTCTCAAAAGGTTTCCATTTCATGGTAGCATTATGCAA
+ttAATAT--CCCAGTGCACC-TCACTACTCTGGCTCACGCTGAGCCTGATGAGGTCATA-GGGCATCCATTTTTCTTGTAGTTATTTTATTATTACCAGAATTCATTCAGAATGGCCACCTCTAGTCAACACAAAATAAGCAGCGTGCTTTGTC-ACTATTTGAGTAGACTTCCTTAGGAAGTTTTCATCTTAGGGTAGTTTTACATAA
+
+3 hg38.chr1 61126 61836 mm39.chr2 111307411 111308031 + 9500
+AGTTCAAGACGTTTAATCAAGACCCTTCACTTACTTAACTATACCTCCTTGAGA-ATCCCATCTATGAAAAAATTCTAGTCATTATAAAAATGATTGATTAAATGAGGGAAGTAGTAGAGTTCTTCATTTCTTTAGTTGGTTTAGTCTCCTATGAGTCAATCCTA-TTTTCAAA-ATTCTTAATAAACCATTTATTCCTTCAACTTTCTATGCCATTTGATGTTTTGTAAAAAAAAAAATATAA------TATGTATACAAAAAGATATTTCAAAATCTAGAAAGAGAGCTTTAGAGCTTTGTAAAGCTCTTTTAAAAATCAAAAACAACTACTGTTAATTAACATGTTGTACTATGC-AATTTGTTTACCATTATTACTCTTGGTATTTTTAAGAAAAGTCTT----TCCATTGTTATTATAAATGCTTCTATTGATATTTATTTTAATAACTGTTATTACAGTCCGTCATGTACATACACTATACTTAAAcctaatgtttggtatttaaatcgtttcaagattttatc----actgtcaacaaagtatgatgaatattttta--tgctgaaaacttctgtaaaaatagaattccaagagtattattgcaccaaaaggcatggacttaaaattcttgatacatgatttcaaaatattttctttaaggtttgaatcagtctatattccctccagcagcgtataaaagtgccaatttctctgatccttagcc
+aGCCCAAGAAAATTAAACAAATCCCTTTACTAATTTAAATGTATCTTGCTACTAGATATCATTTGTGA---------------CATAAATATTTGCAATAAAAT--------CACTAACATTATTTATTTCATCAGTG-------------ATGGGTTAACCCTAGTTTTCAAACAGCTTTGATAAATCTGTCATCTCTTTAGATCACTAGATAATTTGATATTTGAGGGGAAAATAGATATGAATGCACTATGTGCATTAAAATATCCCTTAAA-TCTAAAAATAGAA-TTAGAAGCTCTG------------AAACCTGAAATGCTATCATTGTAACCGGGAATGATGGGATATGCTAATTACTTTACAACAAT----------GTTTTAAGCAAGAATCTTAAAGTCCA------------------------ACATCTAGTTT------TGTTATTATAGCCAACTGTGTATATTTGATATATTTATACCTAATTGCTTTCATTTAAACTGCCTTAAAGCTTTACTTAAAACTATCAATAAAGTCTTGCTAGTATTTGTGACTGGAGATATTATCTGTACAAGTAGTTATTTCAGAGAAGTACTTTACAACAAGGTATGGAGATAAAATTATTGACATAT---------------TGTTTAAGGGTTAAA-----ATATCTTCCCTCTTGTGAATTATAAATAAACCAATTTCATTGAGCTTTTGCT
+
+4 hg38.chr1 61944 64435 mm39.chr2 111308032 111310602 + 69061
+ctccgcctcccggggtcaagctattctcctgcctcagcctcccaagtagctgggactacaggcatgcaccaccatgcccagctaatttttgttatttttagtagagatggagtttccccatgttggacaggatggtctcgatctcttgacctcgtgatccaccctcctcggcctcccaaagtgctgggataacaggcgtgaacaaccatgcccggcctgtaaaactttttcctaatttaacagaaaaataatagtattatattttatcatatttctttgatttctaAGacacacatacacacaca---cacacacatatctgtatatacaaatacacgtatagc---------------------------------------------ttacaTTTT--AATTCTTCATTTCATTTGTTCATTTATTAGGTCTTGGAGATTTTGTGAAAC---------TGTTTAAATTCTTTTTTATACTATGAAGATATCAACCTT-TTGTCTCTACAGCATTTCAAATTCAAGTATGAT----TCACGTGTTGGTTTGGGGTAGATCATTATAGGCACATGTAGGAAACAGCTTT-----------------------------------------------------------------------CAGAGA----------------TGC----------CTTAACCGT-----------------------------------AATTATGCATTTGTATTC--------------TAATTTTTATTTAATGTTATTATT----GATTGCATTTTTAAAG--ATTCTGT--ATTTTTTAAACCATTTATTTGTATATGTTGGTATACAATCTTGCCATTTTCTGGGATTTCATATTTCCTTATTTTTGTTTTT--TACCTTTTTTGGCTTGAATT---TTTTGAGTTTTTATGCATTCTTTTCCAGTTTCTTAAGATGCTAATAAGTTCATGTATTTGAGCAATTGAGAACATTTAAAGCAATAGACTGCCTCTGAGCACAGCTTTGTCCATATTACATTAACCTTTTATACCCTGGGTTCCCACTAGTTTTTAAATAATCTACTATCAAATAAAAG-ATTTGTTAATAATAA-------------------------------------------------------ATTTTAAATCATTAACACTTAACGCATTATTTTCAGTCACACTAAGTTGATTCCTTCGTTTCTTTCAGGTTGCTTCAGAGTCTTCCCTTCTATCTGATTCAGTGGACCAAGTAAATGACTCTCTGGTAACAGAATTTGTATTACTTGGACTTGCACAATCCTTGGAAATGCAGTTTTTCCTTTTTCTCTTCTTCTCTTTATTCTATGTGGGAATTATCCTGGGAAAACTCTTCATTGTGTTCACAGTGATCTTTGATCCTCACTTACACTCCCCCATGTATATTCTGCTGGCCAACCTATCGCTCATTGACTTGAGCCTTTCATCTACCACAGTTCCTAGGTTGATCTACGATCTTTTTACTGATTGTAAAGTTATTTCCTTCCATAATTGCATGATACAAAAGTTCTTTATCCATGTTATGGGAGGAGTTGAAATGGTGCTGCTGATAGTCATGGCATATGATAGGTACACTGCGATCTGCAAGCCTCTCCACTATCCAACTATTATGAATCCCAAAATGTGCATGTTTTTGGTAGCAGCAGCTTGGGTCATTGGGGTGATTCATGCTATGTCTCAGTTTGTTTTTGTCATAAATTTACCCTTCTGTGGCCCTAATAATGTGGGGAGCTTTTATTGTGATTTTCCTCGGGTTATTAAACTTGCATGCATGGACACTTATGGGCTAGAATTTGTGGTCACTGCCAACAGTGGATTCATATCGATGGGCACCTTCTTTTTCTTAATTGTATCATACATTTTTATTCTGGTCACTGTCCAACGACATTCCTCAAATGATTTATCCAAAGCATTCTTCACTTCGTCGGCTCACATCACCGTAGTGGTTTTGTTTTTTGCTCCATGCATGTTTCTCTACGTGTGGCCTTTCCCTACTAAGTCATTGGATAAATTTTTTGCCATCATGAACTTTGTTGTCACCCCTGTCGTAAATCCTGCCATCTATACTTTAAGGAACAAAGATATGAAGTTTGCAATGAGAAGGCTGAATCAACATATTTTAAATTCTATGGAGACGACATAACACATTTGGTTGATGAGAGCACAGGATAAATGCCATGGACCATCAAGACTCCTGTGATCACCATGATCACTATGGAACGCGCACATTTTTAGTATTGCCTGAAAAAACTGAAAAATCTGCAAAAAGGATGCATTAAATCTAAGAATTGTATTTCAGATAAAGTTGCAACATTTTTTGTTAATCATA----AAAAGTA-----------------------------------------TATATTTCTATCT--------AATGTGTGTATCTAATTAA-----------CAGCAATGACTAC-------------CTTT----------------AATTTTGATGTAGTTAT-TTTATATCTGTATATAAGCACATACACATATATATGACCTAGGTTTATTTATCAGTATTTTTATGCTGATAATAAGCATCACTGGAAATTAAT-----TTTCTTATGGAAATTATGTGGATCCAATGGATAAAATATGAGTTTATATAAATTAGTAAATGCCAAAATCAAGGAAGAAACAATTTTTATTTTAATTGTACTTTAAGTTAGATAAATGGTAAGGTCAACAGCTTGTTACAACCCTTAAGTATTATTTTCAGGCTGATTGTCAATATGTTTTGTACAatgttctcac
+CTCTGCTTTGTGTGGAGAAG-TATTGTCC--------ACTCTGAATGATTTGGTAATA-------ATAATAATATATTTTGCTAATT--------------------------------ATATT---------------------------------------------------------------------------------------------TGTAA----TTTTTCTTATTTAACA-AGAAATGAGA-------ATTTTCTTATATTTCACTGATTTCTAacacacacacatatacacatcacatacacacaaagagagagacagagacacagagagacacagagagacagacagacagacagacagacagacagatacacagtcaCATTTCTCAATTCTCTATATCATT-GTTCACTT----TGCCTGATGGGTTTTTAGTAGCAATCATGTATATTTTAAT-CATTTTAATATTATAAAGATATCAGCTTTATTCTCTCTA-AACATTG-AAATTCAAGTACAAGATAAGCATGAGTTAGTTTTGTG---------ACAAGCATGTGTAAGAAATAATTCTTtctctctctctctctctctctctctctctctctctctctctctctccatttttctttctttctttTGGTACAGAGGATGAAAGTTCCACATGTGCTCTTAAAGCACTTTATCATTGAGTCTATTTTTTCAGTTGGCATAAGTAAAAATAAGTTAAGTGTTTTGACTCACCATTCTAATGCTTGATATTCACTTAATTTTATTAGTTTTAAATTTCAGTTTTAAAATAATTTTTTTAAGTTTTGAAAACAAT-----GTTAATGTCCTTACAGAGTCTTTTCATTTTCTGGACTACTGCATTCATTTAGGATTGTTCCCAATATCTCTTTTTATTTCACCTGCACTTGAGAATTTAATACATTTTATCCTACTACAT--AAATCGTGTCATGTTTATGT-TTTTATCAATTGCTATAATTTAATGCAGTATATCAT-TCTGAGTAGAGCTTTTTTCAAACTCCATAAATCTTATATACACTGGATTCCTCACATTTTGTAAATTGTCTAC-ACCAATTAAAAATACTTGTTAATAGAAATTGTTTTGATAAACAAATTGTCTAGTGTTAAACAATATCTCTGGTACAGTGACTGATTTTAAATGACTAACATTTAAT--CTTACTTT--GATGCCTTAAAATAATTCAT-CTTGTCTTTCAGGATTTTTATGTATTTTCCCAACTTTCTGATCCAATGGACCAAGTTAATGCTTCTGCTCTACCTGAATTTGTTTTGCTGGGACTTGCACAATCGTTTGGAACACAGATTTTCTTTGGCCTCTTCTTCTCCTTATTTTATGTGGGGATTCTTTTTGGAAACCTCTTCATTGTGTTCATAGTGATTGTTGATTCTCATTTACACTTCCCCATGTATATTCTATTGGCCAACCTTTCGCTCATTGACTTGGGCCTTTCATCTACAACAATTCCTAGGACAATATCTGATCTTTTTACTGGCTGTAAAGTCATTTCTTTCCACAGCTGCATGACACAAATGTTCTTCATTCATGTGATGGGTGGAGTTGAGATGGTACTGCTCATAGCCATGGCATATGACAGATATATAGCAATCTGCAAGCCTCTTCACTACCTGATGATCATGAATCCCAAAAAATGTATAATTTTGGTAATAGCAGCTTGGGTCATAGGGATGATTCATGCAGTGTCTCAGTTTCTGTTTGTTGTAAATTTACCCTTCTGTGGTCCCTACAATGTAGGAAGTTTTTATTGTGATTTTCCAAGGGTCATTAAACTTGCATGCATGGACACTTATAAACTAGAATTTGTGGTCTCTGCCAACAGTGGCTTTATTTCTATGTGTACCTTCTTTTTTCTGATTACATCATACATTTTTATCCTGGCCAGTGTACGACAACATTCTTCAACTGATTTATCCAAAGCATTTGTCACCTTGTCAGCTCACATCACTGTAGTGGTTTTGTTTTTCATACCATGCATGTTTCTGTATGTATGGCCATTTCCGACCAAGTCACTTGATAATTTCTTTGCTATTGTTGACTTTGTTCTCACTCCTGTCTTAAATCCTACTATCTATACTTTAAGGAATAAAGATATGAGGTTGGCCATCAGAAGGCTAAGTAGACAGGTTTTAAGTTCTAGGGAATTTATATAGTAAATTTTATTAGTATCAACACTGA--AAATGTGATAAACTATCAA---------------TATTAACACTATAGAATAAAAACATTATT-TTATTGGGTTTAGAAACATAAGAATGTAC--------------CTATCAAAGAGTTTCAATTAATAAGGAGTTACAATAACTTCTGCTGGTTACAAAATAAAAGTGTTATATGTTTTCAAATCAACATAACTAAATCTTGTTGGTTTTATATTTTTGTATTGAATTGAAGTATTCATATTTATTTAAGTTATTATAACCAGCAATTATTGCACAAATGTTCATGCTTTTTGGTTTGTATAAATAAATGTGCatatatatatatatatatatatatatatatatatatatatatatataAAATGAACC--TTTATCAGCATCTATATGCAGATAATA--CAACACTGGTTACTAATATATATTTCTTGTGGAAATTATGTGGATCCAATGGATAAA-----GGTTTATGTAAACCAATA------------------------------ATTTTAA------TTTAAGACTTATAAGTAGGA------------------------TAACCATCACTATTAAGAT-------------TTTTGTGTCCTTCTCTTAC