Repository 'maker'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/maker

Changeset 7:7bac6f5fcfa2 (2021-06-22)
Previous changeset 6:d46d803ca6cc (2021-01-05) Next changeset 8:3c1a01127868 (2021-06-28)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/maker commit 38b61f14db4a6e474c55e6107e31eaba6c677ad9"
modified:
macros.xml
maker.xml
added:
test-data/annot_fix_nucleotides.gff3
test-data/evidences_fix_nucletides.gff3
test-data/non_IUPAC.fasta
b
diff -r d46d803ca6cc -r 7bac6f5fcfa2 macros.xml
--- a/macros.xml Tue Jan 05 09:00:59 2021 +0000
+++ b/macros.xml Tue Jun 22 10:15:59 2021 +0000
b
@@ -7,6 +7,7 @@
         </requirements>
     </xml>
     <token name="@VERSION@">2.31.11</token>
+    <token name="@SUFFIX_VERSION@">1</token>
 
     <xml name="citations">
         <citations>
b
diff -r d46d803ca6cc -r 7bac6f5fcfa2 maker.xml
--- a/maker.xml Tue Jan 05 09:00:59 2021 +0000
+++ b/maker.xml Tue Jun 22 10:15:59 2021 +0000
[
@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<tool id="maker" name="Maker" profile="16.04" version="@VERSION@">
+<tool id="maker" name="Maker" profile="16.04" version="@VERSION@+galaxy@SUFFIX_VERSION@">
     <description>genome annotation pipeline</description>
     <macros>
         <import>macros.xml</import>
@@ -44,7 +44,7 @@
             MPI_CMD="mpiexec -n \${GALAXY_SLOTS:-4}";
         fi &&
 
-        \${MPI_CMD} maker --ignore_nfs_tmp maker_opts.ctl maker_bopts.ctl maker_exe.ctl < /dev/null
+        \${MPI_CMD} maker --ignore_nfs_tmp $advanced.fix_nucleotides maker_opts.ctl maker_bopts.ctl maker_exe.ctl < /dev/null
 
         &&
 
@@ -453,6 +453,7 @@
         </section>
 
         <section name="advanced" title="Advanced settings" expanded="False">
+            <param name="fix_nucleotides" argument="-fix_nucleotides" type="boolean" truevalue="-fix_nucleotides" falsevalue="" checked="false" label="Fix nucleotides" help="This will replace non-IUPAC characters with 'N's."/>
             <param name="other_gff" type="data" format="gff" label="Extra features to pass-through to final Maker generated GFF3 file" optional="True"/>
             <param name="alt_peptide" type="text" value="C" size="1" label="Amino acid used to replace non-standard amino acids in BLAST databases">
                 <validator type="regex" message="This must be a single uppercase letter">^[A-Z]$</validator>
@@ -563,6 +564,15 @@
             <output name="output_gff" file="annot_dfam_up.gff3"/>
             <output name="output_evidences" file="evidences_norm_dfam_up.gff3" compare="sim_size"/>
         </test>
+        <test>
+            <param name="genome" value="genome.fasta"/>
+            <param name="est_evidences|est" value="non_IUPAC.fasta"/>
+            <param name="est_evidences|est2genome" value="1"/>
+            <param name="repeat_masking|repeat_source|source_type" value="no"/>
+            <param name="advanced|fix_nucleotides" value="true"/>
+            <output name="output_gff" file="annot_fix_nucleotides.gff3"/>
+            <output name="output_evidences" file="evidences_fix_nucletides.gff3" compare="sim_size"/>
+        </test>
     </tests>
     <help><![CDATA[
         MAKER is a portable and easily configurable genome annotation pipeline. Its purpose is to allow smaller eukaryotic and prokaryotic genome projects to independently annotate their genomes and to create genome databases. MAKER identifies repeats, aligns ESTs and proteins to a genome, produces ab-initio gene predictions and automatically synthesizes these data into gene annotations having evidence-based quality values. MAKER is also easily trainable: outputs of preliminary runs can be used to automatically retrain its gene prediction algorithm, producing higher quality gene-models on seusequent runs. MAKER's inputs are minimal and its ouputs can be directly loaded into a GMOD database. They can also be viewed in the Apollo genome browser; this feature of MAKER provides an easy means to annotate, view and edit individual contigs and BACs without the overhead of a database. MAKER should prove especially useful for emerging model organism projects with minimal bioinformatics expertise and computer resources.
b
diff -r d46d803ca6cc -r 7bac6f5fcfa2 test-data/annot_fix_nucleotides.gff3
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/annot_fix_nucleotides.gff3 Tue Jun 22 10:15:59 2021 +0000
b
@@ -0,0 +1,41 @@
+##gff-version 3
+HS08198 maker gene 352 1848 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0;Name=maker-HS08198-exonerate_est2genome-gene-0.0
+HS08198 maker mRNA 352 1848 2869 + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1;Parent=maker-HS08198-exonerate_est2genome-gene-0.0;Name=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1;_AED=0.00;_eAED=0.00;_QI=70|1|1|1|0|0|7|0|192
+HS08198 maker exon 352 397 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:9;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 421 582 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:10;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 812 894 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:11;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 1053 1123 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:12;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 1208 1315 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:13;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 1587 1688 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:14;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker exon 1772 1848 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:exon:15;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker five_prime_UTR 352 397 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:five_prime_utr;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker five_prime_UTR 421 444 . + . ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:five_prime_utr;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 445 582 . + 0 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 812 894 . + 0 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 1053 1123 . + 1 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 1208 1315 . + 2 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 1587 1688 . + 2 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+HS08198 maker CDS 1772 1848 . + 2 ID=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS08198-exonerate_est2genome-gene-0.0-mRNA-1
+###
+HS04636 maker gene 1813 6903 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0;Name=maker-HS04636-exonerate_est2genome-gene-0.0
+HS04636 maker mRNA 1813 6903 8728 + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1;Parent=maker-HS04636-exonerate_est2genome-gene-0.0;Name=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1;_AED=0.00;_eAED=0.00;_QI=49|1|1|1|0|0|9|0|572
+HS04636 maker exon 1813 1934 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:0;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 2055 2198 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:1;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 2852 2995 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:2;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 3426 3607 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:3;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 4340 4423 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:4;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 4543 4789 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:5;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 5072 5358 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:6;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 5860 6007 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:7;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker exon 6494 6903 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:exon:8;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker five_prime_UTR 1813 1861 . + . ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:five_prime_utr;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 1862 1934 . + 0 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 2055 2198 . + 2 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 2852 2995 . + 2 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 3426 3607 . + 2 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 4340 4423 . + 0 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 4543 4789 . + 0 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 5072 5358 . + 2 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 5860 6007 . + 0 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+HS04636 maker CDS 6494 6903 . + 2 ID=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1:cds;Parent=maker-HS04636-exonerate_est2genome-gene-0.0-mRNA-1
+###
b
diff -r d46d803ca6cc -r 7bac6f5fcfa2 test-data/evidences_fix_nucletides.gff3
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/evidences_fix_nucletides.gff3 Tue Jun 22 10:15:59 2021 +0000
b
@@ -0,0 +1,42 @@
+##gff-version 3
+HS08198 . contig 1 2344 . . . ID=HS08198;Name=HS08198
+###
+HS08198 blastn expressed_sequence_match 444 1848 140 + . ID=HS08198:hit:2:3.2.0.0;Name=HS08198.g2
+HS08198 blastn match_part 444 583 140 + . ID=HS08198:hsp:19:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 79 218 +;Gap=M140
+HS08198 blastn match_part 445 523 79 + . ID=HS08198:hsp:20:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 1 79 +;Gap=M79
+HS08198 blastn match_part 808 895 88 + . ID=HS08198:hsp:21:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 214 301 +;Gap=M88
+HS08198 blastn match_part 1053 1123 71 + . ID=HS08198:hsp:22:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 301 371 +;Gap=M71
+HS08198 blastn match_part 1206 1315 110 + . ID=HS08198:hsp:23:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 370 479 +;Gap=M110
+HS08198 blastn match_part 1586 1689 104 + . ID=HS08198:hsp:24:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 479 582 +;Gap=M104
+HS08198 blastn match_part 1771 1848 78 + . ID=HS08198:hsp:25:3.2.0.0;Parent=HS08198:hit:2:3.2.0.0;Target=HS08198.g2 581 658 +;Gap=M78
+HS08198 est2genome expressed_sequence_match 352 1848 2869 + . ID=HS08198:hit:3:3.2.0.0;Name=HS08198.g2
+HS08198 est2genome match_part 352 397 2869 + . ID=HS08198:hsp:26:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 3 52 +;Gap=M18 D1 M2 I4 M10 I1 M15
+HS08198 est2genome match_part 421 582 2869 + . ID=HS08198:hsp:27:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 53 217 +;Gap=M4 I2 M3 I3 M5 D1 M9 D1 M139
+HS08198 est2genome match_part 812 894 2869 + . ID=HS08198:hsp:28:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 218 300 +;Gap=M83
+HS08198 est2genome match_part 1053 1123 2869 + . ID=HS08198:hsp:29:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 301 371 +;Gap=M71
+HS08198 est2genome match_part 1208 1315 2869 + . ID=HS08198:hsp:30:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 372 479 +;Gap=M108
+HS08198 est2genome match_part 1587 1688 2869 + . ID=HS08198:hsp:31:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 480 581 +;Gap=M102
+HS08198 est2genome match_part 1772 1848 2869 + . ID=HS08198:hsp:32:3.2.0.0;Parent=HS08198:hit:3:3.2.0.0;Target=HS08198.g2 582 658 +;Gap=M77
+HS04636 . contig 1 9453 . . . ID=HS04636;Name=HS04636
+###
+HS04636 blastn expressed_sequence_match 1815 6903 120 + . ID=HS04636:hit:0:3.2.0.0;Name=HS04636.g1
+HS04636 blastn match_part 1815 1934 120 + . ID=HS04636:hsp:0:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 129 248 +;Gap=M120
+HS04636 blastn match_part 1815 1845 31 + . ID=HS04636:hsp:1:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 50 80 +;Gap=M31
+HS04636 blastn match_part 2055 2198 144 + . ID=HS04636:hsp:2:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 249 392 +;Gap=M144
+HS04636 blastn match_part 2852 2997 146 + . ID=HS04636:hsp:3:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 393 538 +;Gap=M146
+HS04636 blastn match_part 3424 3609 186 + . ID=HS04636:hsp:4:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 535 720 +;Gap=M186
+HS04636 blastn match_part 4339 4423 85 + . ID=HS04636:hsp:5:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 718 802 +;Gap=M85
+HS04636 blastn match_part 4542 4790 249 + . ID=HS04636:hsp:6:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 802 1050 +;Gap=M249
+HS04636 blastn match_part 5068 5360 293 + . ID=HS04636:hsp:7:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 1046 1338 +;Gap=M293
+HS04636 blastn match_part 5859 6008 150 + . ID=HS04636:hsp:8:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 1336 1485 +;Gap=M150
+HS04636 blastn match_part 6490 6903 414 + . ID=HS04636:hsp:9:3.2.0.0;Parent=HS04636:hit:0:3.2.0.0;Target=HS04636.g1 1481 1894 +;Gap=M414
+HS04636 est2genome expressed_sequence_match 1813 6903 8728 + . ID=HS04636:hit:1:3.2.0.0;Name=HS04636.g1
+HS04636 est2genome match_part 1813 1934 8728 + . ID=HS04636:hsp:10:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 127 248 +;Gap=M122
+HS04636 est2genome match_part 2055 2198 8728 + . ID=HS04636:hsp:11:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 249 392 +;Gap=M144
+HS04636 est2genome match_part 2852 2995 8728 + . ID=HS04636:hsp:12:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 393 536 +;Gap=M144
+HS04636 est2genome match_part 3426 3607 8728 + . ID=HS04636:hsp:13:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 537 718 +;Gap=M182
+HS04636 est2genome match_part 4340 4423 8728 + . ID=HS04636:hsp:14:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 719 802 +;Gap=M84
+HS04636 est2genome match_part 4543 4789 8728 + . ID=HS04636:hsp:15:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 803 1049 +;Gap=M247
+HS04636 est2genome match_part 5072 5358 8728 + . ID=HS04636:hsp:16:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 1050 1336 +;Gap=M287
+HS04636 est2genome match_part 5860 6007 8728 + . ID=HS04636:hsp:17:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 1337 1484 +;Gap=M148
+HS04636 est2genome match_part 6494 6903 8728 + . ID=HS04636:hsp:18:3.2.0.0;Parent=HS04636:hit:1:3.2.0.0;Target=HS04636.g1 1485 1894 +;Gap=M410
b
diff -r d46d803ca6cc -r 7bac6f5fcfa2 test-data/non_IUPAC.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/non_IUPAC.fasta Tue Jun 22 10:15:59 2021 +0000
b
@@ -0,0 +1,35 @@
+>HS04636.g1
+dddrrrrrrrrrgcccrgctgctgtgcgcggtcctggcgctcagccatacagcaaatccttgctgttcccacccatgtca
+tgctcgcccgcgccctgctgctgtgcgcggtcctggcgctcagccatacagcaaatccttgctgttcccacccatgtcaa
+aaccgaggtgtatgtatgagtgtgggatttgaccagtataagtgcgattgtacccggacaggattctatggagaaaactg
+ctcaacaccggaatttttgacaagaataaaattatttctgaaacccactccaaacacagtgcactacatacttacccact
+tcaagggattttggaacgttgtgaataacattcccttccttcgaaatgcaattatgagttatgtcttgacatccagatca
+catttgattgacagtccaccaacttacaatgctgactatggctacaaaagctgggaagccttctctaacctctcctatta
+tactagagcccttcctcctgtgcctgatgattgcccgactcccttgggtgtcaaaggtaaaaagcagcttcctgattcaa
+atgagattgtggaaaaattgcttctaagaagaaagttcatccctgatccccagggctcaaacatgatgtttgcattcttt
+gcccagcacttcacgcatcagtttttcaagacagatcataagcgagggccagctttcaccaacgggctgggccatggggt
+ggacttaaatcatatttacggtgaaactctggctagacagcgtaaactgcgccttttcaaggatggaaaaatgaaatatc
+agataattgatggagagatgtatcctcccacagtcaaagatactcaggcagagatgatctaccctcctcaagtccctgag
+catctacggtttgctgtggggcaggaggtctttggtctggtgcctggtctgatgatgtatgccacaatctggctgcggga
+acacaacagagtatgcgatgtgcttaaacaggagcatcctgaatggggtgatgagcagttgttccagacaagcaggctaa
+tactgataggagagactattaagattgtgattgaagattatgtgcaacacttgagtggctatcacttcaaactgaaattt
+gacccagaactacttttcaacaaacaattccagtaccaaaatcgtattgctgctgaatttaacaccctctatcactggca
+tccccttctgcctgacacctttcaaattcatgaccagaaatacaactatcaacagtttatctacaacaactctatattgc
+tggaacatggaattacccagtttgttgaatcattcaccaggcaaattgctggcagggttgctggtggtaggaatgttcca
+cccgcagtacagaaagtatcacaggcttccattgaccagagcaggcagatgaaataccagtcttttaatgagtaccgcaa
+acgctttatgctgaagccctatgaatcatttgaagaacttacaggagaaaaggaaatgtctgcagagttggaagcactct
+atggtgacatcgatgctgtggagctgtatcctgcccttctggtagaaaagcctcggccagatgccatctttggtgaaacc
+atggtagaagttggagcaccattctccttgaaaggacttatgggtaatgttatatgttctcctgcctactggaagccaag
+cacttttggtggagaagtgggttttcaaatcatcaacactgcctcaattcagtctctcatctgcaataacgtgaagggct
+gtccctttacttcattcagtgttccagatccagagctcattaaaacagtcaccatcaatgcaagttcttcccgctccgga
+ctagatgatatcaatcccacagtactactaaaagaacgttcgactgaactgtag
+>HS08198.g2
+atgctgccccctgggactgcgaccctcttgactctgctcctggcagctggctcgctgggccagaagcctcagaggccaca
+tgctgccccctgggactgcgaccctcttgactctgctcctggcagctggctcgctgggccagaagcctcagaggccacgc
+cggcccgcatcccccatcagcaccatccagcccaaggccaattttgatgcgcagcaggagcagggccaccgggccgaggc
+caccacactgcatgtggctccccagggcacagccatggctgtcagtaccttccgaaagctggatgggatctgctggcagg
+tgcgccagctctatggagacacaggggtcctcggccgcttcctgcttcaagcccgaggcgcccgaggggctgtgcacgtg
+gttgtcgctgagaccgactaccagagtttcgctgtcctgtacctggagcgggcggggcagctgtcagtgaagctctacgc
+ccgctcgctccctgtgagcgactcggtcctgagtgggtttgagcagcgggtccaggaggcccacctgactgaggaccaga
+tcttctacttccccaagtacggcttctgcgaggctgcagaccagttccacgtcctggacggtgagtgcacagcgggggca
+agcatggcggcgtggtga