Mercurial > repos > iuc > maker

changeset 7:7bac6f5fcfa2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/maker commit 38b61f14db4a6e474c55e6107e31eaba6c677ad9"
author iuc
date Tue, 22 Jun 2021 10:15:59 +0000
parents d46d803ca6cc
children 3c1a01127868
files macros.xml maker.xml test-data/annot_fix_nucleotides.gff3 test-data/evidences_fix_nucletides.gff3 test-data/non_IUPAC.fasta
diffstat 5 files changed, 131 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/macros.xml	Tue Jan 05 09:00:59 2021 +0000
+++ b/macros.xml	Tue Jun 22 10:15:59 2021 +0000
@@ -7,6 +7,7 @@
         </requirements>
     </xml>
     <token name="@VERSION@">2.31.11</token>
+    <token name="@SUFFIX_VERSION@">1</token>
 
     <xml name="citations">
         <citations>
--- 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.
--- /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
@@ -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
+###
--- /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
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/non_IUPAC.fasta	Tue Jun 22 10:15:59 2021 +0000
@@ -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