| Previous changeset 1:271eb7f4b8bc (2022-05-08) |
|
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc commit 908fb9148a3a116dbf7ccd16e0992e7882e748c2 |
|
modified:
hapcut2.xml |
|
added:
macros.xml test-data/ref.fasta |
| b |
| diff -r 271eb7f4b8bc -r 800f8086da7d hapcut2.xml --- a/hapcut2.xml Sun May 08 10:44:03 2022 +0000 +++ b/hapcut2.xml Tue Apr 30 08:38:15 2024 +0000 |
| [ |
| b'@@ -1,17 +1,16 @@\n-<tool id="hapcut2" name="Hapcut2" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@+ga@GA_VERSION_SUFFIX@">\n- <description> - haplotype assembly for diploid organisms</description>\n+<tool id="hapcut2" name="Hapcut2" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="23.2">\n+ <description>haplotype assembly for diploid organisms</description>\n+ <macros>\n+ <token name="@TOOL_VERSION@">1.3.4</token>\n+ <token name="@VERSION_SUFFIX@">1</token>\n+ <import>macros.xml</import> \n+ </macros>\n <xrefs>\n <xref type="bio.tools">hapcut2</xref>\n </xrefs>\n- <macros>\n- <token name="@TOOL_VERSION@">1.3.3</token>\n- <token name="@VERSION_SUFFIX@">0</token>\n- <token name="@GA_VERSION_SUFFIX@">1</token>\n- </macros>\n <requirements>\n- <requirement type="package" version="1.3.3">hapcut2</requirement>\n+ <requirement type="package" version="@TOOL_VERSION@">hapcut2</requirement>\n </requirements>\n-\n <command detect_errors="exit_code"><![CDATA[\n \n ## Prep inputs\n@@ -19,23 +18,31 @@\n ln -s \'$input_bam\' input.bam\n && ln -s \'$input_vcf\' input.vcf\n \n+## Set reference genome if required\n+#if $optimization.choice in [\'pacbio\', \'ont\']:\n+ #if $optimization.reference_genome.source == \'index\':\n+ #set $ref_genome_path = $optimization.reference_genome.index.fields.path\n+ #else\n+ #set $ref_genome_path = $optimization.reference_genome.fasta\n+ #end if\n+#end if\n \n ## Run program\n ## =====================================================================\n ## Extract variant fragments from alignment\n-&& extractHAIRS --bam input.bam --VCF input.vcf --out frags.dat\n+&& extractHAIRS --bam \'input.bam\' --VCF \'input.vcf\' --out frags.dat\n #if $optimization.choice == \'pacbio\':\n --pacbio 1\n- --ref \'$optimization.reference_fasta\'\n+ --ref \'$ref_genome_path\'\n #elif $optimization.choice == \'ont\':\n --ont 1\n- --ref \'$optimization.reference_fasta\'\n+ --ref \'$ref_genome_path\'\n #elif $optimization.choice == \'hic\':\n --HiC 1\n #end if\n \n #if $advanced.minIS\n---minIS $advanced.minIS\n+ --minIS $advanced.minIS\n #end if\n \n #if $advanced.maxIS\n@@ -47,13 +54,11 @@\n #if $optimization.choice == \'hic\':\n --HiC 1\n #end if\n-\n ]]></command>\n \n <inputs>\n- <param name="input_bam" type="data" format="bam" label="Input BAM file"/>\n- <param name="input_vcf" type="data" format="vcf" label="Input VCF file"/>\n-\n+ <param name="input_bam" argument="--bam" type="data" format="bam" label="Input BAM file" help="Coordinate-sorted BAM file"/>\n+ <param name="input_vcf" argument="--VCF" type="data" format="vcf" label="Input VCF file" help="Variant file with genotypes for a single individual"/>\n <conditional name="optimization">\n <!-- TODO: include 10X (requires extra processing step) -->\n <param name="choice" type="select" display="radio" label="Optimization">\n@@ -62,71 +67,111 @@\n <option value="ont">Oxford Nanopore</option>\n <option value="hic">Hi-C</option>\n </param>\n-\n+ <when value="default"></when>\n <when value="pacbio">\n- <param name="reference_fasta" type="data" format="fasta"\n- label="Reference genome fasta file"\n- help="The reference genome is required for long-read optimization."\n- />\n+ <expand macro="reference_genome_input"/>\n </when>\n-\n <when value="ont">\n- <param name="reference_fasta" type="data" format="fasta"\n- label="Reference genome fasta file"\n- help="The reference genome is required for long-read optimization."\n- />\n+ <expand macro="reference_genome_input"/>\n </when>\n+ <when value="hic"></when>\n </conditional>\n \n- <param name="output_phased" type="boolean" label="Output phased VCF file?"\n- checked="true"\n- help="Output variant calls on the haplotype assembly"\n- />\n- <param name="output_fragments" type="boolean" label="Output fragments file?"\n- help="Output fragments collected by extractHAIRS"\n- />\n-\n+ <param name="output_ph'..b'- >\n+ <data name="frags" format="txt" from_work_dir="frags.dat" label="${tool.name} on ${on_string}: Fragments">\n <filter>output_fragments</filter>\n </data>\n </outputs>\n \n <tests>\n+ <!-- Defaults -->\n+ <test expect_num_outputs="1">\n+ <param name="input_bam" ftype="bam" value="input.bam"/>\n+ <param name="input_vcf" ftype="vcf" value="input.vcf"/>\n+ <param name="output_fragments" value="0"/>\n+ <param name="output_phased" value="0"/>\n+ <conditional name="optimization">\n+ <param name="choice" value="default"/>\n+ </conditional>\n+ <output name="haplotype" ftype="txt" file="output_haplotype.out"/>\n+ </test>\n+\n+ <!-- Defaults with all outputs -->\n+ <test expect_num_outputs="3">\n+ <param name="input_bam" ftype="bam" value="input.bam"/>\n+ <param name="input_vcf" ftype="vcf" value="input.vcf"/>\n+ <param name="output_fragments" value="1"/>\n+ <param name="output_phased" value="1"/>\n+ <conditional name="optimization">\n+ <param name="choice" value="default"/>\n+ </conditional>\n+ <output name="frags" ftype="txt" file="output_frag.dat"/>\n+ <output name="haplotype" ftype="txt" file="output_haplotype.out"/>\n+ <output name="haplotype_phased" ftype="vcf" file="output_haplotype.out.phased.vcf"/>\n+ </test>\n+\n+ <!-- Hi-C optimization -->\n <test expect_num_outputs="3">\n <param name="input_bam" ftype="bam" value="input.bam"/>\n <param name="input_vcf" ftype="vcf" value="input.vcf"/>\n <param name="output_fragments" value="1"/>\n <param name="output_phased" value="1"/>\n- <param name="optimization" value="default"/>\n+ <conditional name="optimization">\n+ <param name="choice" value="default"/>\n+ </conditional>\n <output name="frags" ftype="txt" file="output_frag.dat"/>\n <output name="haplotype" ftype="txt" file="output_haplotype.out"/>\n <output name="haplotype_phased" ftype="vcf" file="output_haplotype.out.phased.vcf"/>\n </test>\n+\n+ <!-- Pacbio optimization with ref genome -->\n+ <test expect_num_outputs="3">\n+ <param name="input_bam" ftype="bam" value="input.bam"/>\n+ <param name="input_vcf" ftype="vcf" value="input.vcf"/>\n+ <param name="output_fragments" value="1"/>\n+ <param name="output_phased" value="1"/>\n+ <conditional name="optimization">\n+ <param name="choice" value="pacbio"/>\n+ </conditional>\n+ <conditional name="reference_genome">\n+ <param name="source" value="history"/>\n+ </conditional>\n+ <param name="fasta" ftype="fasta" value="ref.fasta"/>\n+ <output name="frags" ftype="txt" file="output_frag.dat"/>\n+ <output name="haplotype" ftype="txt" file="output_haplotype.out"/>\n+ <output name="haplotype_phased" ftype="vcf" file="output_haplotype.out.phased.vcf"/>\n+ </test>\n+\n+ <!-- Oxford nanopore optimization with ref genome -->\n+ <test expect_num_outputs="3">\n+ <param name="input_bam" ftype="bam" value="input.bam"/>\n+ <param name="input_vcf" ftype="vcf" value="input.vcf"/>\n+ <param name="output_fragments" value="1"/>\n+ <param name="output_phased" value="1"/>\n+ <output name="frags" ftype="txt" file="output_frag.dat"/>\n+ <output name="haplotype" ftype="txt" file="output_haplotype.out"/>\n+ <output name="haplotype_phased" ftype="vcf" file="output_haplotype.out.phased.vcf"/>\n+ <conditional name="optimization">\n+ <param name="optimization" value="ont"/>\n+ </conditional>\n+ <conditional name="reference_genome">\n+ <param name="source" value="history"/>\n+ </conditional>\n+ <param name="fasta" ftype="fasta" value="ref.fasta"/>\n+ </test>\n </tests>\n \n <help><![CDATA[\n@@ -172,6 +217,7 @@\n \n ]]></help>\n <citations>\n- <citation type="doi">https://doi.org/10.1101/gr.213462.116</citation>\n+ <citation type="doi">10.1101/gr.213462.116</citation>\n </citations>\n+ <expand macro="creator"/>\n </tool>\n' |
| b |
| diff -r 271eb7f4b8bc -r 800f8086da7d macros.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros.xml Tue Apr 30 08:38:15 2024 +0000 |
| b |
| @@ -0,0 +1,30 @@ +<macro> + <xml name="reference_genome_input"> + <conditional name="reference_genome"> + <param name="source" type="select" label="Source for the reference genome" help="Built-in references were created using default options."> + <option value="indexed" selected="true">Use a built-in genome</option> + <option value="history">Use a genome from history</option> + </param> + <when value="indexed"> + <param name="index" type="select" label="Select a reference genome" help="If your genome of interest is not listed, contact the Galaxy team."> + <options from_data_table="fasta_indexes"> + <filter type="sort_by" column="2" /> + <validator type="no_options" message="No genomes are available for the selected input dataset" /> + </options> + </param> + </when> + <when value="history"> + <param name="fasta" type="data" format="fasta,fasta.gz" + label="Reference genome fasta file" + help="The reference genome is required for long-read optimization." + /> + </when> + </conditional> + </xml> + + <xml name="creator"> + <creator> + <organization name="Galaxy Australia" url="https://site.usegalaxy.org.au"/> + </creator> + </xml> +</macro> \ No newline at end of file |
| b |
| diff -r 271eb7f4b8bc -r 800f8086da7d test-data/ref.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/ref.fasta Tue Apr 30 08:38:15 2024 +0000 |
| b |
| b'@@ -0,0 +1,428 @@\n+>NC_045512.2 Severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1, complete genome\n+ATTAAAGGTTTATACCTTCCCAGGTAACAAACCAACCAACTTTCGATCTCTTGTAGATCTGTTCTCTAAA\n+CGAACTTTAAAATCTGTGTGGCTGTCACTCGGCTGCATGCTTAGTGCACTCACGCAGTATAATTAATAAC\n+TAATTACTGTCGTTGACAGGACACGAGTAACTCGTCTATCTTCTGCAGGCTGCTTACGGTTTCGTCCGTG\n+TTGCAGCCGATCATCAGCACATCTAGGTTTCGTCCGGGTGTGACCGAAAGGTAAGATGGAGAGCCTTGTC\n+CCTGGTTTCAACGAGAAAACACACGTCCAACTCAGTTTGCCTGTTTTACAGGTTCGCGACGTGCTCGTAC\n+GTGGCTTTGGAGACTCCGTGGAGGAGGTCTTATCAGAGGCACGTCAACATCTTAAAGATGGCACTTGTGG\n+CTTAGTAGAAGTTGAAAAAGGCGTTTTGCCTCAACTTGAACAGCCCTATGTGTTCATCAAACGTTCGGAT\n+GCTCGAACTGCACCTCATGGTCATGTTATGGTTGAGCTGGTAGCAGAACTCGAAGGCATTCAGTACGGTC\n+GTAGTGGTGAGACACTTGGTGTCCTTGTCCCTCATGTGGGCGAAATACCAGTGGCTTACCGCAAGGTTCT\n+TCTTCGTAAGAACGGTAATAAAGGAGCTGGTGGCCATAGTTACGGCGCCGATCTAAAGTCATTTGACTTA\n+GGCGACGAGCTTGGCACTGATCCTTATGAAGATTTTCAAGAAAACTGGAACACTAAACATAGCAGTGGTG\n+TTACCCGTGAACTCATGCGTGAGCTTAACGGAGGGGCATACACTCGCTATGTCGATAACAACTTCTGTGG\n+CCCTGATGGCTACCCTCTTGAGTGCATTAAAGACCTTCTAGCACGTGCTGGTAAAGCTTCATGCACTTTG\n+TCCGAACAACTGGACTTTATTGACACTAAGAGGGGTGTATACTGCTGCCGTGAACATGAGCATGAAATTG\n+CTTGGTACACGGAACGTTCTGAAAAGAGCTATGAATTGCAGACACCTTTTGAAATTAAATTGGCAAAGAA\n+ATTTGACACCTTCAATGGGGAATGTCCAAATTTTGTATTTCCCTTAAATTCCATAATCAAGACTATTCAA\n+CCAAGGGTTGAAAAGAAAAAGCTTGATGGCTTTATGGGTAGAATTCGATCTGTCTATCCAGTTGCGTCAC\n+CAAATGAATGCAACCAAATGTGCCTTTCAACTCTCATGAAGTGTGATCATTGTGGTGAAACTTCATGGCA\n+GACGGGCGATTTTGTTAAAGCCACTTGCGAATTTTGTGGCACTGAGAATTTGACTAAAGAAGGTGCCACT\n+ACTTGTGGTTACTTACCCCAAAATGCTGTTGTTAAAATTTATTGTCCAGCATGTCACAATTCAGAAGTAG\n+GACCTGAGCATAGTCTTGCCGAATACCATAATGAATCTGGCTTGAAAACCATTCTTCGTAAGGGTGGTCG\n+CACTATTGCCTTTGGAGGCTGTGTGTTCTCTTATGTTGGTTGCCATAACAAGTGTGCCTATTGGGTTCCA\n+CGTGCTAGCGCTAACATAGGTTGTAACCATACAGGTGTTGTTGGAGAAGGTTCCGAAGGTCTTAATGACA\n+ACCTTCTTGAAATACTCCAAAAAGAGAAAGTCAACATCAATATTGTTGGTGACTTTAAACTTAATGAAGA\n+GATCGCCATTATTTTGGCATCTTTTTCTGCTTCCACAAGTGCTTTTGTGGAAACTGTGAAAGGTTTGGAT\n+TATAAAGCATTCAAACAAATTGTTGAATCCTGTGGTAATTTTAAAGTTACAAAAGGAAAAGCTAAAAAAG\n+GTGCCTGGAATATTGGTGAACAGAAATCAATACTGAGTCCTCTTTATGCATTTGCATCAGAGGCTGCTCG\n+TGTTGTACGATCAATTTTCTCCCGCACTCTTGAAACTGCTCAAAATTCTGTGCGTGTTTTACAGAAGGCC\n+GCTATAACAATACTAGATGGAATTTCACAGTATTCACTGAGACTCATTGATGCTATGATGTTCACATCTG\n+ATTTGGCTACTAACAATCTAGTTGTAATGGCCTACATTACAGGTGGTGTTGTTCAGTTGACTTCGCAGTG\n+GCTAACTAACATCTTTGGCACTGTTTATGAAAAACTCAAACCCGTCCTTGATTGGCTTGAAGAGAAGTTT\n+AAGGAAGGTGTAGAGTTTCTTAGAGACGGTTGGGAAATTGTTAAATTTATCTCAACCTGTGCTTGTGAAA\n+TTGTCGGTGGACAAATTGTCACCTGTGCAAAGGAAATTAAGGAGAGTGTTCAGACATTCTTTAAGCTTGT\n+AAATAAATTTTTGGCTTTGTGTGCTGACTCTATCATTATTGGTGGAGCTAAACTTAAAGCCTTGAATTTA\n+GGTGAAACATTTGTCACGCACTCAAAGGGATTGTACAGAAAGTGTGTTAAATCCAGAGAAGAAACTGGCC\n+TACTCATGCCTCTAAAAGCCCCAAAAGAAATTATCTTCTTAGAGGGAGAAACACTTCCCACAGAAGTGTT\n+AACAGAGGAAGTTGTCTTGAAAACTGGTGATTTACAACCATTAGAACAACCTACTAGTGAAGCTGTTGAA\n+GCTCCATTGGTTGGTACACCAGTTTGTATTAACGGGCTTATGTTGCTCGAAATCAAAGACACAGAAAAGT\n+ACTGTGCCCTTGCACCTAATATGATGGTAACAAACAATACCTTCACACTCAAAGGCGGTGCACCAACAAA\n+GGTTACTTTTGGTGATGACACTGTGATAGAAGTGCAAGGTTACAAGAGTGTGAATATCACTTTTGAACTT\n+GATGAAAGGATTGATAAAGTACTTAATGAGAAGTGCTCTGCCTATACAGTTGAACTCGGTACAGAAGTAA\n+ATGAGTTCGCCTGTGTTGTGGCAGATGCTGTCATAAAAACTTTGCAACCAGTATCTGAATTACTTACACC\n+ACTGGGCATTGATTTAGATGAGTGGAGTATGGCTACATACTACTTATTTGATGAGTCTGGTGAGTTTAAA\n+TTGGCTTCACATATGTATTGTTCTTTCTACCCTCCAGATGAGGATGAAGAAGAAGGTGATTGTGAAGAAG\n+AAGAGTTTGAGCCATCAACTCAATATGAGTATGGTACTGAAGATGATTACCAAGGTAAACCTTTGGAATT\n+TGGTGCCACTTCTGCTGCTCTTCAACCTGAAGAAGAGCAAGAAGAAGATTGGTTAGATGATGATAGTCAA\n+CAAACTGTTGGTCAACAAGACGGCAGTGAGGACAATCAGACAACTACTATTCAAACAATTGTTGAGGTTC\n+AACCTCAATTAGAGATGGAACTTACACCAGTTGTTCAGACTATTGAAGTGAATAGTTTTAGTGGTTATTT\n+AAAACTTACTGACAATGTATACATTAAAAATGCAGACATTGTGGAAGAAGCTAAAAAGGTAAAACCAACA\n+GTGGTTGTTAATGCAGCCAATGTTTACCTTAAACATGGAGGAGGTGTTGCAGGAGCCTTAAATAAGGCTA\n+CTAACAATGCCATGCAAGTTGAATCTGATGATTACATAGCTACTAATGGACCACTTAAAGTGGGTGGTAG\n+TTGTGTTTTAAGCGGACACAATCTTGCTAAACACTGTCTTCATGTTGTCGGCCCAAATGTTAACAAAGGT\n+GAAGACATTCAACTTCTTAAGAGTGCTTATGAAAATTTTAATCAGCACGAAGTTCTACTTGCACCATTAT\n+TATCAGCTGGTATTTTTGGTGCTGACCCTATACATTCTTTAAGAGTTTGTGTAGATACTGTTCGCA'..b'AAAAGACTGTGTTGTATTACACAGTTACTTCACTTCAGACTATTACCAGCTGTACTCA\n+ACTCAATTGAGTACAGACACTGGTGTTGAACATGTTACCTTCTTCATCTACAATAAAATTGTTGATGAGC\n+CTGAAGAACATGTCCAAATTCACACAATCGACGGTTCATCCGGAGTTGTTAATCCAGTAATGGAACCAAT\n+TTATGATGAACCGACGACGACTACTAGCGTGCCTTTGTAAGCACAAGCTGATGAGTACGAACTTATGTAC\n+TCATTCGTTTCGGAAGAGACAGGTACGTTAATAGTTAATAGCGTACTTCTTTTTCTTGCTTTCGTGGTAT\n+TCTTGCTAGTTACACTAGCCATCCTTACTGCGCTTCGATTGTGTGCGTACTGCTGCAATATTGTTAACGT\n+GAGTCTTGTAAAACCTTCTTTTTACGTTTACTCTCGTGTTAAAAATCTGAATTCTTCTAGAGTTCCTGAT\n+CTTCTGGTCTAAACGAACTAAATATTATATTAGTTTTTCTGTTTGGAACTTTAATTTTAGCCATGGCAGA\n+TTCCAACGGTACTATTACCGTTGAAGAGCTTAAAAAGCTCCTTGAACAATGGAACCTAGTAATAGGTTTC\n+CTATTCCTTACATGGATTTGTCTTCTACAATTTGCCTATGCCAACAGGAATAGGTTTTTGTATATAATTA\n+AGTTAATTTTCCTCTGGCTGTTATGGCCAGTAACTTTAGCTTGTTTTGTGCTTGCTGCTGTTTACAGAAT\n+AAATTGGATCACCGGTGGAATTGCTATCGCAATGGCTTGTCTTGTAGGCTTGATGTGGCTCAGCTACTTC\n+ATTGCTTCTTTCAGACTGTTTGCGCGTACGCGTTCCATGTGGTCATTCAATCCAGAAACTAACATTCTTC\n+TCAACGTGCCACTCCATGGCACTATTCTGACCAGACCGCTTCTAGAAAGTGAACTCGTAATCGGAGCTGT\n+GATCCTTCGTGGACATCTTCGTATTGCTGGACACCATCTAGGACGCTGTGACATCAAGGACCTGCCTAAA\n+GAAATCACTGTTGCTACATCACGAACGCTTTCTTATTACAAATTGGGAGCTTCGCAGCGTGTAGCAGGTG\n+ACTCAGGTTTTGCTGCATACAGTCGCTACAGGATTGGCAACTATAAATTAAACACAGACCATTCCAGTAG\n+CAGTGACAATATTGCTTTGCTTGTACAGTAAGTGACAACAGATGTTTCATCTCGTTGACTTTCAGGTTAC\n+TATAGCAGAGATATTACTAATTATTATGAGGACTTTTAAAGTTTCCATTTGGAATCTTGATTACATCATA\n+AACCTCATAATTAAAAATTTATCTAAGTCACTAACTGAGAATAAATATTCTCAATTAGATGAAGAGCAAC\n+CAATGGAGATTGATTAAACGAACATGAAAATTATTCTTTTCTTGGCACTGATAACACTCGCTACTTGTGA\n+GCTTTATCACTACCAAGAGTGTGTTAGAGGTACAACAGTACTTTTAAAAGAACCTTGCTCTTCTGGAACA\n+TACGAGGGCAATTCACCATTTCATCCTCTAGCTGATAACAAATTTGCACTGACTTGCTTTAGCACTCAAT\n+TTGCTTTTGCTTGTCCTGACGGCGTAAAACACGTCTATCAGTTACGTGCCAGATCAGTTTCACCTAAACT\n+GTTCATCAGACAAGAGGAAGTTCAAGAACTTTACTCTCCAATTTTTCTTATTGTTGCGGCAATAGTGTTT\n+ATAACACTTTGCTTCACACTCAAAAGAAAGACAGAATGATTGAACTTTCATTAATTGACTTCTATTTGTG\n+CTTTTTAGCCTTTCTGCTATTCCTTGTTTTAATTATGCTTATTATCTTTTGGTTCTCACTTGAACTGCAA\n+GATCATAATGAAACTTGTCACGCCTAAACGAACATGAAATTTCTTGTTTTCTTAGGAATCATCACAACTG\n+TAGCTGCATTTCACCAAGAATGTAGTTTACAGTCATGTACTCAACATCAACCATATGTAGTTGATGACCC\n+GTGTCCTATTCACTTCTATTCTAAATGGTATATTAGAGTAGGAGCTAGAAAATCAGCACCTTTAATTGAA\n+TTGTGCGTGGATGAGGCTGGTTCTAAATCACCCATTCAGTACATCGATATCGGTAATTATACAGTTTCCT\n+GTTTACCTTTTACAATTAATTGCCAGGAACCTAAATTGGGTAGTCTTGTAGTGCGTTGTTCGTTCTATGA\n+AGACTTTTTAGAGTATCATGACGTTCGTGTTGTTTTAGATTTCATCTAAACGAACAAACTAAAATGTCTG\n+ATAATGGACCCCAAAATCAGCGAAATGCACCCCGCATTACGTTTGGTGGACCCTCAGATTCAACTGGCAG\n+TAACCAGAATGGAGAACGCAGTGGGGCGCGATCAAAACAACGTCGGCCCCAAGGTTTACCCAATAATACT\n+GCGTCTTGGTTCACCGCTCTCACTCAACATGGCAAGGAAGACCTTAAATTCCCTCGAGGACAAGGCGTTC\n+CAATTAACACCAATAGCAGTCCAGATGACCAAATTGGCTACTACCGAAGAGCTACCAGACGAATTCGTGG\n+TGGTGACGGTAAAATGAAAGATCTCAGTCCAAGATGGTATTTCTACTACCTAGGAACTGGGCCAGAAGCT\n+GGACTTCCCTATGGTGCTAACAAAGACGGCATCATATGGGTTGCAACTGAGGGAGCCTTGAATACACCAA\n+AAGATCACATTGGCACCCGCAATCCTGCTAACAATGCTGCAATCGTGCTACAACTTCCTCAAGGAACAAC\n+ATTGCCAAAAGGCTTCTACGCAGAAGGGAGCAGAGGCGGCAGTCAAGCCTCTTCTCGTTCCTCATCACGT\n+AGTCGCAACAGTTCAAGAAATTCAACTCCAGGCAGCAGTAGGGGAACTTCTCCTGCTAGAATGGCTGGCA\n+ATGGCGGTGATGCTGCTCTTGCTTTGCTGCTGCTTGACAGATTGAACCAGCTTGAGAGCAAAATGTCTGG\n+TAAAGGCCAACAACAACAAGGCCAAACTGTCACTAAGAAATCTGCTGCTGAGGCTTCTAAGAAGCCTCGG\n+CAAAAACGTACTGCCACTAAAGCATACAATGTAACACAAGCTTTCGGCAGACGTGGTCCAGAACAAACCC\n+AAGGAAATTTTGGGGACCAGGAACTAATCAGACAAGGAACTGATTACAAACATTGGCCGCAAATTGCACA\n+ATTTGCCCCCAGCGCTTCAGCGTTCTTCGGAATGTCGCGCATTGGCATGGAAGTCACACCTTCGGGAACG\n+TGGTTGACCTACACAGGTGCCATCAAATTGGATGACAAAGATCCAAATTTCAAAGATCAAGTCATTTTGC\n+TGAATAAGCATATTGACGCATACAAAACATTCCCACCAACAGAGCCTAAAAAGGACAAAAAGAAGAAGGC\n+TGATGAAACTCAAGCCTTACCGCAGAGACAGAAGAAACAGCAAACTGTGACTCTTCTTCCTGCTGCAGAT\n+TTGGATGATTTCTCCAAACAATTGCAACAATCCATGAGCAGTGCTGACTCAACTCAGGCCTAAACTCATG\n+CAGACCACACAAGGCAGATGGGCTATATAAACGTTTTCGCTTTTCCGTTTACGATATATAGTCTACTCTT\n+GTGCAGAATGAATTCTCGTAACTACATAGCACAAGTAGATGTAGTTAACTTTAATCTCACATAGCAATCT\n+TTAATCAGTGTGTAACATTAGGGAGGACTTGAAAGAGCCACCACATTTTCACCGAGGCCACGCGGAGTAC\n+GATCGAGTGTACAGTGAACAATGCTAGGGAGAGCTGCCTATATGGAAGAGCCCTAATGTGTAAAATTAAT\n+TTTAGTAGTGCTATCCCCATGTGATTTTAATAGCTTCTTAGGAGAATGAC\n' |