Mercurial > repos > greg > kaks_analysis
changeset 27:f174450ebc44 draft
Uploaded
| author | greg |
|---|---|
| date | Fri, 28 Apr 2017 14:23:46 -0400 |
| parents | 73db26d39092 |
| children | b807167c1e60 |
| files | kaks_analysis.xml macros.xml test-data/species1.faa test-data/species1.fna test-data/species2.faa test-data/species2.fna |
| diffstat | 6 files changed, 117 insertions(+), 128 deletions(-) [+] |
line wrap: on
line diff
--- a/kaks_analysis.xml Tue Apr 11 13:36:06 2017 -0400 +++ b/kaks_analysis.xml Fri Apr 28 14:23:46 2017 -0400 @@ -79,17 +79,17 @@ ]]> </command> <inputs> - <param name="coding_sequences_species_1" format="fasta" type="data" label="Coding sequences (CDS) fasta file for species1" /> - <param name="proteins_species_1" format="fasta" type="data" label="Aamino acids (proteins) sequences fasta file for species1" /> + <param name="coding_sequences_species_1" format="fasta" type="data" label="Coding sequences for the first species" /> + <param name="proteins_species_1" format="fasta" type="data" label="Protein sequences for the first species" /> <conditional name="comparison_cond"> - <param name="comparison" type="select" label="Select method for pairwise sequence comparison to determine homolgous pairs" help="Cross species comparison requires selection of inputs for second species"> - <option value="paralogs" selected="true">Self species comparison</option> - <option value="orthologs">Cross species comparison</option> + <param name="comparison" type="select" label="Type of sequence comparison"> + <option value="paralogs" selected="true">Paralogous</option> + <option value="orthologs">Orthologous</option> </param> <when value="paralogs" /> <when value="orthologs"> - <param name="coding_sequences_species_2" format="fasta" type="data" label="Coding sequences (CDS) fasta file for species2" /> - <param name="proteins_species_2" format="fasta" type="data" label="Aamino acids (proteins) sequences fasta file for species2" /> + <param name="coding_sequences_species_2" format="fasta" type="data" label="Coding sequences for the second species" /> + <param name="proteins_species_2" format="fasta" type="data" label="Protein sequences for the second species" /> </when> </conditional> <conditional name="options_type"> @@ -100,27 +100,27 @@ <when value="basic" /> <when value="advanced"> <conditional name="set_min_coverage_cond"> - <param name="set_min_coverage" type="select" label="Specify minimum sequence pairwise coverage length between homologous pairs?"> + <param name="set_min_coverage" type="select" label="Alignment coverage configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="min_coverage" type="float" value="0.5" min="0.3" max="1.0" label="Minimum sequence pairwise coverage length between homologous pairs" /> + <param name="min_coverage" type="float" value="0.5" min="0.3" max="1.0" label="match score" /> </when> </conditional> <conditional name="recalibrate_cond"> - <param name="recalibrate" type="select" label="Specify evolutionary rate for recalibrating synonymous subsitutions (ks) of species?"> + <param name="recalibrate" type="select" label="Species rates recalibration configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="recalibration_rate" type="float" value="0.0" min="0.0" label="Evolutionary rate for recalibrating synonymous subsitutions (ks) of species" /> + <param name="recalibration_rate" type="float" value="0.0" min="0.0" label="Recalibration rate" /> </when> </conditional> <conditional name="codeml_ctl_file_cond"> - <param name="codeml_ctl_file_select" type="select" label="Select PAML codeml control file?" help="Used for ML analysis of protein-coding DNA sequences using codon substitution models, select No to use the default control file"> + <param name="codeml_ctl_file_select" type="select" label="PAML codeml configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> @@ -130,33 +130,33 @@ </when> </conditional> <conditional name="fit_components_cond"> - <param name="fit_components" type="select" label="Fit a mixture model of multivariate normal components to synonymous (ks) distribution?" help="Used to identify significant duplication events in a genome"> + <param name="fit_components" type="select" label="Rates clustering configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="num_of_components" type="integer" value="1" min="1" label="Number of components to fit to synonymous subsitutions (ks) distribution" /> + <param name="num_of_components" type="integer" value="1" min="1" label="Number of components" /> </when> </conditional> <conditional name="set_lower_ks_limit_cond"> - <param name="set_lower_ks_limit" type="select" label="Set lower limit of synonymous subsitutions (ks)?" help="Reduces background noise from young paralogous pairs due to normal gene births and deaths in a genome"> + <param name="set_lower_ks_limit" type="select" label="Lower limit synonymous subsitution rates configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="min_ks" type="float" value="0.0" min="0.0" label="Lower limit of synonymous subsitutions (ks)" /> + <param name="min_ks" type="float" value="0.0" min="0.0" label="Minimum rate" /> </when> </conditional> <conditional name="set_upper_ks_limit_cond"> - <param name="set_upper_ks_limit" type="select" label="Set upper limit of synonymous subsitutions (ks)?" help="Excludes likey ancient paralogous pairs"> + <param name="set_upper_ks_limit" type="select" label="Upper limit synonymous subsitution rates configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="max_ks" type="float" value="0.0" min="0.0" label="Upper limit of synonymous subsitutions (ks)" /> + <param name="max_ks" type="float" value="0.0" min="0.0" label="Maximum rate" /> </when> </conditional> </when> @@ -195,9 +195,11 @@ </test> </tests> <help> -This tool is one of the PlantTribes collection of automated modular analysis pipelines that utilize objective classifications of -complete protein sequences from sequenced plant genomes to perform comparative evolutionary studies. This tool performs orthologous -or paralogous ks analyses of coding sequences and amino acid sequences. +This tool is one of the PlantTribes collection of automated modular analysis pipelines for comparative and evolutionary analyses +of genome-scale gene families and transcriptomes. This tool estimates paralogous and orthologous pairwise synonymous (Ks) and +non-synonymous (Ka) substitution rates for a set of gene coding sequences either produced by the AssemblyPostProcessor tool or +from an external source. Optionally, the resulting set of estimated Ks values can be clustered into components using a mixture +of multivariate normal distributions to identify significant duplication event(s) in a species or a pair of species. ----- @@ -205,26 +207,91 @@ * **Required** - - **Coding sequences (CDS) fasta file for species1** - Coding sequences (CDS) fasta file for species1. - - **Aamino acids (proteins) sequences fasta file for species1** - Aamino acids (proteins) sequences fasta file for species1 - - **Select method for pairwise sequence comparison to determine homolgous pairs** - Pairwise sequence comparison to determine homolgous pairs (cross species comparison requires selection of inputs for species2). + - **Coding sequences for the first species** - coding sequence fasta file for for the first species either produced by the AssemblyPostProcessor tool or an external source selected from your history. + - **Protein sequences for the first species** - corresponding protein sequence fasta files for the first species either produced by the AssemblyPostProcessor tool or an external source selected from your history. + - **Type of sequence comparison** - pairwise sequence comparison to determine homolgous pairs. This can be either paralogous for self-species comparison or orthologous for cross-species comparison. Cross species comparision requires data selected for the second species. * **Optional** - - **Minimum sequence pairwise coverage length between homologous pairs** - Minimum sequence pairwise coverage length between homologous pairs (e.g., 0.5 results in 50% coverage. Legal values lie between 0.3 and 1.0. - - **Evolutionary rate for recalibrating synonymous subsitutions (ks) of species** - (applies to paralogous ks analysis) Recalibrate synonymous subsitutions (ks) of species using a predetermined evoutionary rate that can be determined from a species tree inferred from a collection single copy genes from taxa of interest (Cui et al., 2006). - - **Select PAML codeml control file?** - Select PAML's codeml control file from your history. This file is used to to perfom ML analysis of protein-coding DNA sequences using codon substitution models. Selecting No uses the default file which does not include input (seqfile, treefile) and output (outfile) parameters of codeml. - - **Fit a mixture model of multivariate normal components to synonymous (ks) distribution?** - Fit a mixture model of multivariate normal components to synonymous (ks) distribution to identify significant duplication event(s) in a genome. - - **Number components to fit to synonymous subsitutions (ks) distribution** - Number components to fit to synonymous subsitutions (ks) distribution. - - **Lower limit of synonymous subsitutions (ks)** - Lower limit of synonymous subsitutions (ks) - necessary if fitting components to the distribution to reduce background noise from young paralogous pairs due to normal gene births and deaths in a genome. - - **Upper limit of synonymous subsitutions (ks)** - Upper limit of synonymous subsitutions (ks) - necessary if fitting components to the distribution to exclude likey ancient paralogous pairs. + - **Coding sequences for the second species** - coding sequence fasta file for for the second species either produced by the AssemblyPostProcessor tool or an external source selected from your history. Required only for orthologous comparison. + - **Protein sequences for the second species** - corresponding protein sequence fasta files for the second species either produced by the AssemblyPostProcessor tool or an external source selected from your history. Required only for orthologous comparison. + - **Alignment coverage configuration** - select 'Yes' to set the minimum allowable alignment coverage length between homologous pairs. PlantTribes uses global codon alignment match score to determine the pairwise alignment coverage. By default, the match score is set to 0.5 if 'No' is selected. + + - **match score** - number of base matches in a pairwise sequence alignment divided by the length of shorter sequence. Positions in the alignment corresponding to gaps are not considered. The score is restricted to the range 0.3 - 1.0. + + - **Species rates recalibration configuration** - select 'Yes' to recalibrate synonymous substitution rates of a species using a predetermined evolutionary rate. Recalibration evolutionary rate can be determined from a species tree inferred from a collection of conserved single copy genes from taxa of interest as described in Cui et al., 2006. Applies only to paralogous comparisons. + + - **recalibration rate** - a predetermined evolutionary recalibration rate. + + - **PAML codeml configuration** - select 'Yes' to enable selection of a PAML codeml control file to carry out maximum likelihood analysis of protein-coding DNA sequences using codon substitution models. Template file "codeml.ctl.args" can be found in the scaffold data installed into Galaxy via the PlantTribes Scaffolds Download Data Manager tool, and are also available at the PlantTribes GitHub `repository`_. Default settings shown in the template are used if 'No' is selected. + +.. _repository: https://github.com/dePamphilis/PlantTribes/blob/master/config/codeml.ctl.args + + - **Rates clustering configuration** - select 'Yes' to estimate clusters of synonymous substitution rates using a mixture of multivariate normal distributions which represent putative duplication event(s). + + - **Number of components** - number of components to include in the normal mixture model. + + - **Lower limit synonymous subsitution rates configuration** - select 'Yes' to set the minimum allowable synonymous substitution rate to use in the normal mixtures cluster analysis to exclude young paralogs that arise from normal gene births and deaths in a genome. + + - **Minimum rate** - minimum allowable synonymous substitution rate. + + - **Upper limit synonymous subsitution rates configuration** - select 'Yes' to set the maximum allowable synonymous substitution rate to use in the normal mixtures cluster analysis to exclude likely ancient paralogs in a genome. + + + - **Maximum rate** - maximum allowable synonymous substitution rate. </help> <citations> <expand macro="citation1" /> - <citation type="doi">10.1093/bioinformatics/btw412</citation> - <citation type="doi">10.1186/1471-2105-10-421</citation> - <citation type="doi">10.1093/molbev/msm088</citation> - <citation type="doi">10.18637/jss.v004.i02</citation> + <citation type="bibtex"> + @article{Wall2008, + journal = {Nucleic Acids Research}, + author = {2. Wall PK, Leebens-Mack J, Muller KF, Field D, Altman NS}, + title = {PlantTribes: a gene and gene family resource for comparative genomics in plants}, + year = {2008}, + volume = {36}, + number = {suppl 1}, + pages = {D970-D976},} + </citation> + <citation type="bibtex"> + @article{Altschul1990, + journal = {Journal of molecular biology} + author = {3. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ}, + title = {Basic local alignment search tool}, + year = {1990}, + volume = {215}, + number = {3}, + pages = {403-410},} + </citation> + <citation type="bibtex"> + @article{Katoh2013, + journal = {Molecular biology and evolution}, + author = {4. Katoh K, Standley DM}, + title = {MAFFT multiple sequence alignment software version 7: improvements in performance and usability}, + year = {2013}, + volume = {30}, + number = {4}, + pages = {772-780},} + </citation> + <citation type="bibtex"> + @article{Yang2007, + journal = {Molecular biology and evolution}, + author = {5. Yang Z}, + title = {PAML 4: phylogenetic analysis by maximum likelihood}, + year = {2007}, + volume = {24}, + number = {8}, + pages = {1586-1591},} + </citation> + <citation type="bibtex"> + @article{McLachlan1999, + journal = {Journal of Statistical Software}, + author = {6. McLachlan GJ, Peel D, Basford KE, Adams P}, + title = {The EMMIX software for the fitting of mixtures of normal and t-components}, + year = {1999}, + volume = {4}, + number = {2}, + pages = {1-14},} + </citation> </citations> </tool>
--- a/macros.xml Tue Apr 11 13:36:06 2017 -0400 +++ b/macros.xml Fri Apr 28 14:23:46 2017 -0400 @@ -3,7 +3,7 @@ <token name="@WRAPPER_VERSION@">0.8</token> <xml name="requirements_assembly_post_processor"> <requirements> - <requirement type="package" version="0.4">plant_tribes_assembly_post_processor</requirement> + <requirement type="package" version="0.8">plant_tribes_assembly_post_processor</requirement> </requirements> </xml> <xml name="requirements_gene_family_aligner"> @@ -40,7 +40,7 @@ </stdio> </xml> <xml name="param_codon_alignments"> - <param name="codon_alignments" type="select" label="Construct orthogroup multiple codon alignments?"> + <param name="codon_alignments" type="select" label="Construct orthogroup multiple codon alignments"> <option value="yes" selected="true">Yes</option> <option value="no">No</option> </param> @@ -59,13 +59,13 @@ </param> </xml> <xml name="param_orthogroup_fna"> - <param name="orthogroup_fna" type="select" label="Process corresponding gene family classification orthogroups CDS fasta files?"> + <param name="orthogroup_fna" type="select" label="Orthogroups coding sequences"> <option value="yes" selected="true">Yes</option> <option value="no">No</option> </param> </xml> <xml name="param_scaffold"> - <param name="scaffold" type="select" label="Orthogroups or gene families proteins scaffold"> + <param name="scaffold" type="select" label="Gene family scaffold"> <options from_data_table="plant_tribes_scaffolds" /> <validator type="no_options" message="No PlantTribes scaffolds are available. Use the PlantTribes Scaffolds Download Data Manager tool in Galaxy to install and populate the PlantTribes scaffolds data table." /> </param> @@ -78,9 +78,9 @@ </xml> <xml name="cond_alignment_method"> <conditional name="alignment_method_cond"> - <param name="alignment_method" type="select" force_select="true" label="Select method for multiple sequence alignments"> - <option value="mafft" selected="true">MAFFT algorithm</option> - <option value="pasta">PASTA algorithm</option> + <param name="alignment_method" type="select" force_select="true" label="Multiple sequence alignment method"> + <option value="mafft" selected="true">MAFFT</option> + <option value="pasta">PASTA</option> </param> <when value="mafft" /> <when value="pasta"> @@ -90,31 +90,31 @@ </xml> <xml name="cond_remove_gappy_sequences"> <conditional name="remove_gappy_sequences_cond"> - <param name="remove_gappy_sequences" type="select" label="Remove gappy sequences in alignments?"> + <param name="remove_gappy_sequences" type="select" label="Alignment post-processing configuration"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> <conditional name="trim_type_cond"> - <param name="trim_type" type="select" label="Select process used for trimming"> - <option value="gap_trimming" selected="true">Remove gappy sites in alignments (gap trimming)</option> - <option value="automated_trimming">Trim alignments using trimAl's ML heuristic trimming approach (automated trimming)</option> + <param name="trim_type" type="select" label="Trimming method"> + <option value="gap_trimming" selected="true">Gap score based trimming</option> + <option value="automated_trimming">Automated heuristic trimming</option> </param> <when value="gap_trimming"> - <param name="gap_trimming" type="float" value="0" min="0" max="1.0" label="Remove sites in alignments with gaps of" help="Zero value has no affect" /> + <param name="gap_trimming" type="float" optional="true" min="0" max="1.0" label="Gap score" /> </when> <when value="automated_trimming" /> </conditional> <conditional name="remove_sequences_with_gaps_cond"> - <param name="remove_sequences_with_gaps" type="select" label="Remove sequences with specified gaps?"> + <param name="remove_sequences_with_gaps" type="select" label="Remove sequences"> <option value="no" selected="true">No</option> <option value="yes">Yes</option> </param> <when value="no" /> <when value="yes"> - <param name="remove_sequences_with_gaps_of" type="float" value="0" min="0" max="1" label="Remove sequences with gaps of" help="Zero value has no affect" /> - <param name="iterative_realignment" type="integer" value="0" min="0" label="Maximum number of iterations" help="Zero value has no affect"/> + <param name="remove_sequences_with_gaps_of" type="float" optional="true" min="0" max="1" label="Coverage score" /> + <param name="iterative_realignment" type="integer" optional="true" min="0" label="Realignment iteration limit" /> </when> </conditional> </when> @@ -130,34 +130,4 @@ url = {https://github.com/dePamphilis/PlantTribes},} </citation> </xml> - <xml name="citations2to4"> - <citation type="bibtex"> - @article{Sasidharan2012, - journal = {Nucleic Acids Research}, - author = {2. Sasidharan R, Nepusz T, Swarbreck D, Huala E, Paccanaro A}, - title = {GFam: a platform for automatic annotation of gene families}, - year = {2012}, - pages = {gks631},} - </citation> - <citation type="bibtex"> - @article{Li2003, - journal = {Genome Research} - author = {3. Li L, Stoeckert CJ, Roos DS}, - title = {OrthoMCL: identification of ortholog groups for eukaryotic genomes}, - year = {2003}, - volume = {13}, - number = {9}, - pages = {2178-2189},} - </citation> - <citation type="bibtex"> - @article{Emms2015, - journal = {Genome Biology} - author = {4. Emms DM, Kelly S}, - title = {OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy}, - year = {2015}, - volume = {16}, - number = {1}, - pages = {157},} - </citation> - </xml> </macros>
--- a/test-data/species1.faa Tue Apr 11 13:36:06 2017 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,16 +0,0 @@ ->species1_11 -MGVSMGQGNPMGMHLLPSGSSSPRTSPSLRDPPLSLPVLPNSDLSVSLPDLHKLRRNEPVTSGILHVRDLSFLRPRSHNGDDDEETEEMTREQEEKYLQWRSSLVEKLAGIELNLERVKFRMSVEIPPSDDFRAMKKSWENFYASELLSSRNPVRKIAKRPDTILVRGVPSRWFAETRISSKASTLVTHTIIESC ->species1_12 -MSAAAAALRPTEPLPLPSGLSLAPRLKLLLTFFRADLSVRPVDEWQLKTALLAFLRDPPLSLPVLPDSDLSVRTLPDLHKRRRDEPVASGVLHVRDLSFLRPRRRNGDDEEEEAEEMTREQEEEKYFQWRSSLVEKLAGIELNLEGVKFRMSVEIPPSDDFRAMKKSWENFYASELLSSRNPVRKIAKRPDTILVRGVPSRWFAETRISSKASTLVTHTIFSALGKIRNLNISSDDEWGAKQDGTNKEIISGLNCKVWVQFENYDDFNSAMQALCGRSLEKEGSRLKVDYEVTWDHEGFFRNAQYEPVRSNLEERNSSAHGRKKHYTSRIESDHRKRFRD ->species1_15 -MKDGLSLSFALISSSPDSKCELLNSRPSCRAARRGESGLLIRRSYLRPCQCPFGDRMSEQQDSTSKSSSSSISSSTQESEEEVSITIGSLLAQAKNNSGHSLGRRLSQLGSIPHTPRVNGKIPNLDNATLDHERLSERLGNYGLAEFQIEGDGNCQFRALADQIFRNPDYHKHVRKLVMKQLKEFRKQYESYVPMEYKVYLKKMKRSGEWGDHLTLQAAADRFGAKICLLTSFRDTCLIEIVPRDVTPTRELWLSFWCEVHYNSLYATDDLLTRKTKKKHWLF ->species1_16 -MSEQQDHASKSSCSSLSTSTQESEEDVTVGTLLTEAKNSGRSLGKRLSHLDSIPHTPRVNGQIPDVNNATIDHETLLERLGTYGLAEFQIEGDGNCQFRALADQIFRNPDYHKHVRKSVVKQLKEFRKHYEGYVPMEYKVYLKKMKRSGEWGDHVTLQAAADRFAAKICLLTSFRDTCLIEIVPRGATPTKELWLSFWSEVHYNSLYATEDLPNRKTRKKHWLF ->species1_21 -MAGAGAGESLDLPVVDLASSDLAAAAKSVRKACVEYGFFYVVNHGAEGLAEKVFGESSKFFEQPLGEKMALLRNRNYLGYTPLGADKLDASSKFKGDLNENYCIGPIRKEGYQNDANQWPSEENFPCWKETMKLYHETALATGKRILSLIALSLNLDVEFFDCPVAFLRLLHYPGEANESDDGNYGASAHSDYGVLTLVATDGTPGLQICREKDRCPQLWEDVHHIEGALIVNIGDLLQRWTNCVFRSTLHRVVAVGKERYSVAFFLHTNPDLVVQCLESCCSEACPPRFPPIRSGDYLEDRLRARYK ->species1_22 -MWGPHIILYLQPFFLLPSSHMSCVLGRPSAPSLDHPQQPNPPPVAPEKPPAVAKKAAEEEEEKKPPKQARRERHAWSSRSAAAEAVGLGLGGSFANRARGEQVAAGWPAWLSAVVGEAIDGWTLRRADSFEKIDKVRTPALALAIVGGGGRELSSSVLSVAQIGQGTYINVYKARDTVTGKIVALKKMGQVCFLLCKPSYRGDTAAGGRGGRRRQQQQTAALAEEESGMAGGGGGGNRLDLPVVDLASSDPRAAAESIRKACVESGFFYVVNHGVEEGLLKRLFAESSKFFELPMEEKIALRRNSNHRGYTPPYAEKLDPSSKFEGDLKESFYIGPIGDEGLQNDANQWPSEERLPSRRETIKMYHASALSTGKRILSLIALSLNLDAEFFENIGAFSCPSAFLRLLHYPGEVDDSDDGNYGASAHSDYGMITLLATDGTPGLQICREKNRNPQLWEDVHHIDGALIVNIGDLLERWTNCIYRSTVHRVVAVGKERYSAAFFLDPNPDLVVQCLESCCSESCPPRFSPIKSGDYLKERLSATYK ->species1_35 -MAAATTSRRGPGAMDDENLTFETSPGVEVISSFDQMGIRDDLLRGIYAYGFEKPSAIQQRAVLPIISGRDVIAQAQSGTGKTSMISLSVCQIVDTAVREVQALILSPTRELAAQTERVMLAIGDFINIQVHACIGGKSIGEDIRKLEHGVHVVSGTPGRVCDMIKRRTLRTRAIKLLILDEADEMLGRGFKDQIYDVYRYLPPELQVCLISATLPHEILEMTSKFMTDPVRILVKRDELTLEGIKQFFVAVEKEEWKFDTLCDLYDTLTITQAVIFCNTKRKVDWLTERMRSNNFTVSAMHGDMPQKERDAIMGEFRSGATRVLITTDVWARGLDVQQVSLVINYDLPNNRELYIHRIGRSGRFGRKGVAINFVKKEDIRILRDIEQYYSTQIDEMPMNVADLI ->species1_36 -MAAATTSRRGPGAMDDENLTFETSPGVEVISSFDQMGIREDLLRGIYAYGFEKPSAIQQRAVLPIISGRDVIAQAQSGTGKTSMISLSVCQIVDTAVREVQALILSPTRELAAQTERVMLAIGDYINIQVHACIGGKSIGEDIRKLEHGVHVVSGTPGRVCDMIKRRTLRTRAIKLLILDEADEMLGRGFKDQIYDVYRYLPPELQVCLISATLPHEILEMTSKFMTDPVRILVKRDELTLEGIKQFFVAVEKEEWKFDTLCDLYDTLTITQAVIFCNTKRKVDWLTERMRSNNFTVSAMHGDMPQKERDAIMGEFRSGATRVLITTDVWARGLDVQQVSLVINYDLPNNRELYIHRIGRSGRFGRKGVAINFVKKEDIRILRDIEQYYSTQIDEMPMNVADLI
--- a/test-data/species1.fna Tue Apr 11 13:36:06 2017 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,16 +0,0 @@ ->species1_11 -ATGGGTGTGAGTATGGGGCAAGGAAACCCAATGGGTATGCACTTGTTGCCATCTGGCAGCTCAAGTCCGCGCACCTCGCCTTCCCTCCGCGACCCGCCCCTCTCCCTCCCCGTCCTTCCCAACTCCGACCTCTCCGTGTCCCTCCCCGACCTGCATAAGCTTCGCCGCAATGAGCCCGTCACTTCGGGCATCCTCCACGTCCGCGACCTCTCATTCCTCCGCCCCCGCAGCCACAACGGGGATGATGATGAGGAGACCGAGGAGATGACCCGTGAGCAGGAGGAGAAGTACTTGCAGTGGAGGAGCTCCCTGGTCGAGAAGCTGGCCGGGATCGAGCTCAACCTCGAGAGGGTTAAGTTTCGGATGAGCGTCGAAATCCCGCCCTCCGATGACTTCAGGGCAATGAAGAAGTCTTGGGAGAATTTCTACGCCTCCGAGCTCCTCAGTAGCAGGAATCCTGTGAGGAAGATAGCGAAAAGGCCAGACACAATTCTTGTCCGTGGTGTGCCATCCAGGTGGTTTGCGGAGACGAGGATATCATCGAAAGCCTCCACACTGGTCACACACACTATCATCGAAAGCTGC ->species1_12 -ATGTCCGCCGCCGCCGCCGCCCTCCGGCCGACCGAGCCGCTCCCCCTCCCGAGCGGCCTCTCCCTCGCGCCGCGCCTCAAGCTGCTCCTCACCTTCTTCCGCGCCGACCTCTCCGTCCGCCCCGTCGACGAGTGGCAGCTCAAGACCGCGCTCCTCGCCTTCCTCCGCGACCCGCCCCTCTCCCTCCCCGTCCTCCCCGACTCCGACCTCTCCGTGCGCACCCTCCCCGACCTGCATAAGCGCCGCCGCGACGAGCCCGTCGCCTCGGGCGTCCTCCACGTCCGCGACCTCTCCTTCCTCCGCCCACGCCGCCGCAACGGGGATGATGAGGAGGAGGAGGCCGAGGAGATGACCCGTGAGCAGGAGGAGGAGAAGTACTTCCAGTGGAGGAGCTCCCTGGTCGAGAAGCTGGCCGGGATCGAGCTCAACCTCGAGGGGGTTAAGTTTCGGATGAGCGTCGAGATCCCGCCCTCCGATGACTTCAGGGCAATGAAGAAGTCTTGGGAGAATTTCTACGCCTCCGAGCTCCTCAGTAGCAGGAATCCTGTGAGGAAGATAGCGAAAAGGCCAGACACCATTCTTGTCCGGGGTGTGCCATCCAGGTGGTTTGCGGAGACGAGGATATCATCGAAAGCCTCCACGCTGGTCACACACACTATTTTCTCGGCACTTGGTAAAATAAGGAACCTTAATATTTCTAGTGATGATGAATGGGGAGCAAAACAAGACGGAACCAATAAGGAGATTATATCTGGACTAAATTGCAAAGTGTGGGTGCAATTTGAGAACTACGACGATTTCAACAGTGCAATGCAGGCATTATGTGGACGTTCATTAGAAAAAGAAGGATCACGGTTGAAGGTAGACTATGAAGTAACTTGGGATCATGAAGGTTTCTTCCGCAATGCACAATACGAGCCTGTTCGCAGCAATTTAGAAGAGAGAAATTCATCGGCTCATGGAAGGAAGAAACATTACACATCGCGAATTGAGTCAGATCATAGAAAGAGATTTAGGGAT ->species1_15 -ATGAAAGATGGCCTTTCTCTCTCCTTCGCTCTCATCAGCTCGAGCCCCGACAGCAAGTGTGAGCTACTGAACTCGAGACCCTCCTGTCGCGCGGCGCGGCGCGGCGAGAGTGGCCTTTTGATCCGACGAAGCTATCTAAGACCCTGCCAATGTCCATTTGGAGATAGGATGTCGGAACAGCAGGATAGTACTAGTAAAAGCTCTAGCTCAAGCATCAGCAGCAGTACACAGGAGAGCGAGGAGGAGGTATCTATAACAATAGGTAGCCTCCTCGCCCAAGCAAAGAACAACAGTGGGCATAGTCTTGGAAGGCGCCTCTCTCAATTGGGTTCAATCCCGCACACTCCTCGAGTTAATGGAAAAATCCCTAATCTTGATAATGCAACTTTGGATCATGAAAGATTGTCGGAAAGGTTGGGAAATTATGGTTTGGCCGAGTTTCAAATAGAGGGTGATGGAAATTGTCAGTTCCGAGCTTTGGCAGACCAGATATTTCGCAACCCCGATTATCACAAACATGTGAGAAAGTTAGTCATGAAACAGCTAAAGGAATTCAGAAAACAGTATGAAAGCTATGTACCTATGGAATATAAAGTCTACTTGAAGAAAATGAAAAGATCTGGGGAATGGGGGGATCATCTGACTTTACAAGCAGCTGCAGACAGGTTTGGTGCCAAAATTTGTTTGCTGACGTCATTCAGAGACACCTGCCTAATTGAGATAGTCCCCAGGGATGTGACTCCCACAAGGGAGTTGTGGCTAAGCTTCTGGTGTGAAGTGCACTACAATTCCTTGTACGCAACTGACGATCTCCTAACCCGCAAAACCAAGAAGAAGCATTGGTTGTTC ->species1_16 -ATGTCTGAACAACAGGATCATGCTAGCAAAAGTTCTTGCTCAAGTCTTAGCACCAGTACTCAGGAGAGTGAGGAGGATGTGACAGTTGGTACCCTTTTAACTGAAGCAAAGAACAGTGGACGGAGTCTTGGAAAACGCCTTTCCCACTTAGATTCTATCCCGCACACTCCTCGAGTTAATGGGCAAATTCCTGATGTTAATAATGCAACAATAGACCATGAAACATTACTGGAAAGATTGGGCACTTATGGCTTAGCTGAATTCCAAATTGAAGGAGACGGAAATTGTCAGTTCCGAGCTTTGGCAGATCAGATATTCCGCAATCCTGACTATCACAAACATGTGAGGAAGTCAGTCGTGAAGCAGCTAAAGGAATTCAGGAAACACTATGAAGGCTATGTACCGATGGAATATAAGGTGTACTTGAAGAAAATGAAAAGATCTGGAGAATGGGGAGATCATGTGACCTTACAAGCGGCTGCAGACCGGTTTGCTGCCAAGATTTGCCTGCTGACATCATTTAGAGACACATGCCTAATCGAGATAGTCCCCAGAGGTGCCACTCCCACAAAAGAGCTTTGGTTAAGCTTCTGGAGTGAGGTGCACTACAATTCCTTGTATGCAACTGAAGATCTTCCAAATCGCAAGACCAGAAAGAAGCACTGGCTGTTC ->species1_21 -ATGGCCGGCGCCGGCGCCGGCGAGAGCCTGGACCTCCCCGTGGTGGACCTAGCGTCCTCCGACCTCGCCGCCGCCGCCAAATCCGTCCGAAAGGCTTGCGTGGAGTACGGATTCTTCTACGTGGTCAACCATGGAGCCGAGGGATTGGCGGAGAAGGTGTTCGGGGAGAGCAGCAAGTTTTTCGAGCAGCCGCTGGGGGAGAAGATGGCGCTGCTGAGGAACAGAAACTACCTGGGGTACACCCCGCTTGGCGCCGATAAGCTCGACGCCTCGTCCAAATTCAAAGGAGATCTCAATGAAAATTACTGTATCGGACCTATCAGAAAAGAAGGTTATCAGAATGATGCTAACCAATGGCCTTCTGAAGAGAATTTCCCATGTTGGAAGGAGACAATGAAGCTATACCATGAAACTGCACTTGCTACTGGTAAAAGGATACTCTCTCTAATTGCTCTGAGTTTGAATCTCGACGTTGAATTCTTTGACTGCCCAGTGGCCTTTCTTCGGTTATTGCACTACCCAGGTGAAGCTAACGAGTCCGATGATGGCAATTATGGTGCATCAGCTCACTCAGACTATGGAGTACTAACACTTGTAGCAACAGATGGCACTCCTGGGCTGCAGATATGCAGGGAGAAGGATAGGTGCCCCCAGCTTTGGGAAGACGTTCATCACATTGAAGGGGCCCTGATTGTTAATATCGGCGATTTGCTACAAAGGTGGACTAATTGTGTTTTCAGGTCTACACTGCATCGCGTTGTTGCAGTTGGTAAAGAGCGATACTCTGTGGCTTTCTTTCTTCACACAAACCCTGATTTAGTGGTTCAATGCTTGGAAAGCTGCTGCAGTGAGGCATGCCCACCGAGGTTCCCACCTATAAGGAGCGGCGACTATTTGGAAGACCGATTGAGGGCTAGATACAAA ->species1_22 -ATGTGGGGCCCACATATCATCCTCTATCTCCAACCCTTCTTCCTCCTCCCTTCCTCTCACATGAGCTGCGTCCTCGGCCGCCCCTCCGCCCCCTCCCTCGACCACCCCCAGCAGCCCAACCCCCCGCCCGTCGCCCCGGAGAAGCCGCCCGCCGTCGCCAAGAAGGCGGCCGAGGAGGAGGAGGAGAAGAAGCCGCCGAAGCAGGCTAGGAGGGAGAGGCACGCATGGTCGTCGCGGTCTGCCGCCGCCGAGGCGGTCGGCCTGGGGCTCGGGGGGAGCTTCGCCAACAGGGCGCGCGGGGAGCAGGTGGCGGCCGGCTGGCCCGCCTGGCTCTCCGCCGTCGTCGGCGAGGCCATCGACGGCTGGACCCTGCGCCGCGCCGACTCCTTCGAGAAGATCGACAAGGTACGTACTCCTGCCCTCGCGCTCGCCATTGTTGGTGGTGGGGGAAGGGAACTGAGCTCATCGGTCTTGTCGGTGGCGCAGATCGGGCAGGGGACGTACATCAACGTGTACAAGGCGCGGGACACGGTGACGGGCAAGATCGTGGCGCTCAAGAAGATGGGCCAAGTTTGCTTCCTTCTCTGTAAGCCCAGTTACCGTGGGGATACAGCCGCCGGCGGACGCGGAGGGCGGCGGCGGCAGCAGCAGCAAACCGCCGCTTTGGCAGAAGAGGAATCCGGGATGGCCGGCGGCGGCGGCGGCGGGAATCGCCTGGACCTCCCCGTGGTGGACCTCGCGTCCTCCGACCCCCGAGCCGCCGCCGAGTCCATCCGAAAGGCGTGCGTGGAGTCCGGATTCTTCTACGTGGTCAACCATGGGGTGGAGGAGGGATTGCTGAAGAGGTTGTTCGCGGAGAGCTCGAAGTTCTTCGAGCTGCCGATGGAGGAGAAGATAGCGCTGCGGAGGAACAGCAACCACCGGGGATACACCCCGCCCTACGCCGAGAAGCTCGATCCCTCGTCCAAATTCGAAGGAGACCTCAAGGAAAGTTTCTATATTGGGCCTATTGGAGATGAAGGTTTGCAGAATGATGCTAACCAGTGGCCTTCTGAAGAGCGCTTACCAAGTCGGAGGGAGACAATTAAGATGTACCATGCAAGTGCACTGTCTACTGGCAAAAGGATACTCTCTCTAATCGCTCTGAGTTTGAATCTTGACGCTGAATTCTTTGAGAACATTGGTGCCTTCAGCTGCCCATCAGCATTTCTTCGATTATTGCACTACCCAGGTGAAGTAGACGACTCTGATGATGGCAATTATGGTGCATCAGCTCACTCTGATTATGGAATGATAACCCTCCTAGCAACAGACGGCACTCCTGGGCTACAGATATGCAGGGAAAAGAATAGGAATCCCCAGCTCTGGGAAGATGTTCATCACATTGATGGGGCCCTGATTGTTAACATTGGCGATTTGCTAGAAAGGTGGACGAATTGTATTTACAGGTCTACAGTGCACCGTGTTGTTGCAGTTGGTAAAGAGCGATATTCTGCGGCTTTTTTTCTTGACCCAAACCCTGATTTAGTGGTTCAGTGTTTGGAAAGCTGTTGCAGCGAGTCATGCCCACCGAGGTTCTCACCTATAAAGAGTGGCGACTATTTGAAAGAGCGATTGAGCGCTACATACAAA ->species1_35 -ATGGCGGCGGCCACCACGTCGCGGCGCGGCCCGGGCGCCATGGACGACGAGAACCTCACCTTCGAGACCTCCCCGGGGGTCGAGGTCATCAGCAGCTTCGACCAGATGGGGATCCGCGACGACCTCCTCCGCGGCATCTACGCCTACGGCTTCGAGAAGCCCTCCGCCATCCAGCAGCGCGCCGTCCTCCCCATCATCAGCGGCCGCGACGTCATCGCCCAGGCCCAGTCCGGGACCGGCAAGACCTCCATGATCTCGCTCTCCGTCTGCCAGATCGTAGACACCGCCGTCCGTGAGGTGCAGGCTTTAATACTGTCACCAACTAGAGAACTTGCTGCACAAACAGAAAGAGTTATGCTGGCTATCGGTGACTTCATCAATATCCAAGTGCATGCTTGTATTGGTGGCAAAAGTATTGGTGAGGATATTAGAAAGCTTGAGCACGGAGTGCATGTGGTGTCAGGAACACCTGGCAGAGTCTGTGATATGATCAAGAGAAGGACCTTGCGTACAAGAGCCATTAAGCTCCTAATTCTGGATGAAGCTGATGAGATGTTGGGCAGAGGCTTTAAGGATCAGATATATGATGTGTACAGATACCTCCCTCCAGAACTCCAGGTTTGCTTGATCTCCGCAACTCTGCCTCACGAGATCTTGGAAATGACCAGCAAGTTCATGACTGATCCAGTTCGGATCCTTGTGAAGCGTGATGAATTGACTCTAGAGGGCATCAAACAATTCTTTGTTGCTGTTGAGAAAGAAGAATGGAAGTTTGACACGCTTTGTGATCTTTATGATACACTGACAATCACCCAAGCTGTCATTTTCTGCAACACAAAGAGAAAGGTTGATTGGCTTACGGAAAGAATGCGCAGCAATAACTTCACAGTATCAGCTATGCATGGCGACATGCCTCAAAAGGAAAGGGATGCCATTATGGGTGAATTCAGGTCTGGTGCAACCCGTGTTCTAATCACGACAGATGTGTGGGCTCGAGGCCTCGATGTTCAGCAGGTCTCTCTTGTCATAAATTATGATCTCCCAAATAATCGTGAACTTTACATCCATCGCATTGGTCGCTCTGGACGTTTTGGTCGCAAGGGTGTGGCCATCAATTTTGTCAAAAAGGAAGACATCCGTATCCTGAGAGATATCGAGCAGTACTACAGCACGCAGATTGATGAAATGCCAATGAATGTTGCTGATCTAATT ->species1_36 -ATGGCGGCGGCCACCACGTCCCGGCGCGGCCCCGGCGCCATGGACGACGAGAACCTCACCTTCGAGACCTCCCCCGGGGTCGAGGTCATCAGCAGCTTCGACCAGATGGGGATCCGCGAGGACCTCCTCCGCGGCATCTACGCCTACGGCTTCGAGAAGCCCTCCGCCATCCAGCAGCGCGCCGTCCTCCCCATCATCAGCGGCCGCGACGTCATCGCCCAGGCCCAGTCCGGAACCGGCAAGACCTCCATGATCTCGCTCTCCGTCTGCCAGATCGTCGACACCGCCGTCCGAGAGGTTCAGGCCTTGATACTCTCACCAACTAGAGAACTTGCTGCACAAACAGAAAGAGTTATGCTGGCCATTGGTGATTACATCAATATCCAAGTGCATGCTTGTATTGGTGGCAAAAGTATTGGTGAGGATATTAGAAAGCTTGAGCATGGAGTGCATGTTGTGTCAGGAACACCTGGCAGAGTCTGTGATATGATCAAGAGAAGGACCTTGCGTACAAGAGCCATTAAGCTCCTAATTCTGGATGAAGCCGATGAGATGTTGGGCAGAGGCTTTAAGGATCAGATATATGATGTCTACAGATATCTACCCCCAGAGCTCCAGGTTTGCTTGATCTCCGCAACTCTGCCACATGAGATCTTGGAAATGACCAGCAAGTTCATGACTGACCCAGTCCGGATCCTTGTAAAGCGTGATGAATTGACCCTAGAGGGCATCAAACAATTCTTTGTTGCTGTTGAGAAAGAAGAATGGAAGTTTGATACTCTTTGTGATCTTTATGATACACTGACAATCACCCAAGCTGTCATTTTCTGCAACACGAAGAGAAAGGTTGATTGGCTTACAGAAAGAATGCGCAGCAATAACTTCACGGTATCAGCTATGCATGGTGACATGCCTCAAAAGGAAAGGGATGCCATTATGGGTGAATTCAGGTCTGGTGCAACCCGTGTTCTAATTACGACAGATGTGTGGGCTCGAGGCCTGGATGTTCAGCAGGTCTCTCTTGTCATAAACTATGATCTTCCAAATAATCGTGAACTTTACATCCATCGCATTGGTCGCTCTGGACGTTTTGGTCGCAAGGGTGTGGCCATCAATTTTGTCAAAAAGGAAGACATCCGTATCCTGAGAGATATTGAGCAGTACTACAGCACACAGATTGATGAAATGCCAATGAATGTTGCTGATCTAATT
--- a/test-data/species2.faa Tue Apr 11 13:36:06 2017 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,8 +0,0 @@ ->species2_5 -MESQSAVVPLIAELPEKRGGKTLVEEVWEESKKLWEVTGPAAFTGMVLYSMTIVSQAFAGHLGDRHLAAFSIANTVISGLNFGILLGMASALETLCGQAYGAKQYSMMGTYLQRSWLVLLAFAVLLAPTYIFSGQLLMVLGQPAELSREAGLLGMYLLPLHLMFAIQLPLNKFLQCQRKNWVIALSSVLGFPVHVVATWLLAQRFQLGVLGAAMSLNLSWALITGLQLAYAVGGGCPETWRGFSSSAFMGLKDFVSLSVASGVMTCLESWYYRLLIFLTAYAKNAELAVDALSICLSWAGWEMMIHFGFLAGTGVRVANELGANNGRAAKFATIVSTTTSFLICLLISSLALIFHDKLAILFTSSEAVIDAVDGISVLLALTILLNGIQPVLSGVAVGSGWQALVAYVNIGSYYIIGVPFGVLLAWGFHYGVLGIWVGMIGGTMVQTLILSFITLRCDWNEEALKASSRMRTWSSSK ->species2_6 -MEENRSDIPLISGSELPDRRGGGKISELAKEVWGESKKLWVVAGPAAFTRLTFYGMTVVSQAFAGHIGDLELAAFSIATTVISGLSFGFFVGMASAMETLCGQAYGAKQYHMMGIYLQRSWLILLSFAVLLTPTYIFSEQLLTALGQPAELSRQAGLVSLYMLPLHFVYAIVLPLNKFLQCQRKNWVAAVTTAAAFPVHVVATWLLVRCFRLGVFGAAMALTLSWALATVGLLSYALGGGCPETWRGFSASAFVDLKDFIKLSAASGVMLCLENWYYRILVFLTGYVKNAELAVDALSICISYAGWEMMIHLGFLAGTGVRVANELGAANGARARFATIVSMTTSFLISLFISLLILIFHDKLGMIFSSSQAVIDAVDNISFLLALTILLNGIQPVLSGVAVGSGWQALVAYVNIGSYYLIGVPFGFLLGWGLHYGVQGIWVGMIVGTMVQTLILAYITLRCDWNEEALKASTRMRRWSNSK ->species2_9 -MGTLGGHVAPGAFFFLIGLWHLFGHSRLFLLQRGSYVAPVWFPVPGVRHIELIMIIIGSVISVSMELVIVQPKHQPFDDDGTIPSVHLHNFEHASISLAWLVFAAATIHMDRVRAPMRDAVSQLAAAAAFAQQLLIFHFHSADHAGVQGRYHRLLEMVVAVTLAASLLLIPYQRSIALSLVRSASLVFQGVWFTVMGVMMWTPALVPKGCFMNDEDGLQVVRCRTDEALDRAKSLVNLQFNWYLTGTVAFVVVFYLQMAKQYQEQPQYAPLVKGGRGSDGRCTIGEVNDDEDDLEASKGGLGYIEIER ->species2_10 -MGTLVGHVAPGAGFLLIGLWQLFSHIRLFLLRPSSYSAPVWFPAPGVRHLELILIIIGAAMSILMELVIGPAKHQPFDDDGTIPSDHLHNFEHASISLALLVFAAVTIHLDRVKAPLRDAVSQLVAAAAFAQQLLIFHLHSADHMGVEGQFHWLLQTVIAVTLATTLLGIPYPRSIVVSLVRSASLVLQGVWFVVMGVMLWTPALIPKGCFLNLEEGHDVVRCRTDEALDRAKSLVNLQFSWYLTGTVVFVVLFYLQMAKLYPEEPQYLPLVKGGGGGGDDRDSRFSIGDDDHDDEDDVEAAKRGFGHVVSGTKPVEIER
--- a/test-data/species2.fna Tue Apr 11 13:36:06 2017 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,8 +0,0 @@ ->species2_5 -ATGGAGAGTCAGAGCGCCGTCGTCCCGCTCATTGCCGAGCTCCCGGAGAAGCGGGGAGGCAAAACCCTGGTGGAGGAGGTATGGGAGGAGTCCAAGAAGCTGTGGGAAGTCACCGGCCCGGCCGCCTTTACGGGGATGGTACTCTACAGCATGACCATCGTCAGCCAGGCCTTCGCCGGCCACCTTGGTGACCGCCATCTCGCCGCTTTCTCCATCGCCAACACCGTCATATCTGGCCTTAACTTTGGCATTTTGCTTGGCATGGCGAGTGCGCTGGAGACATTATGCGGCCAAGCCTACGGTGCAAAGCAGTACTCGATGATGGGCACCTATCTCCAGCGCTCATGGCTCGTCCTCCTCGCCTTCGCGGTGCTCCTTGCTCCGACGTACATCTTCAGCGGGCAGCTGCTCATGGTCCTGGGCCAGCCCGCCGAGCTGTCTCGCGAGGCGGGCTTGCTCGGCATGTACCTGCTCCCGCTGCACCTCATGTTTGCCATCCAGCTGCCGCTCAACAAGTTCTTGCAGTGCCAGCGCAAGAACTGGGTCATCGCGCTGTCCTCGGTGCTGGGTTTCCCGGTGCACGTCGTGGCGACCTGGCTGCTGGCGCAGCGCTTTCAGCTTGGCGTCCTGGGCGCAGCGATGTCACTCAACCTGTCCTGGGCGCTCATCACGGGCCTGCAGCTCGCGTACGCTGTTGGCGGTGGGTGCCCAGAGACGTGGAGAGGGTTCTCGTCGTCGGCATTCATGGGCTTGAAGGACTTCGTCAGCTTGTCCGTCGCGTCGGGAGTCATGACGTGCTTGGAGAGTTGGTACTACCGGTTATTGATTTTCCTAACGGCGTACGCGAAGAACGCAGAATTGGCTGTGGATGCACTGTCTATCTGCTTGAGTTGGGCTGGATGGGAGATGATGATTCATTTCGGGTTCTTAGCAGGCACTGGGGTGAGGGTTGCCAATGAGCTAGGCGCCAATAATGGACGAGCTGCAAAGTTTGCGACGATCGTGTCCACGACGACATCATTCCTGATCTGCCTCTTAATTAGTTCACTCGCACTCATTTTCCATGACAAACTCGCAATACTGTTCACGTCTAGTGAGGCTGTGATCGATGCAGTTGACGGTATTTCTGTTCTGCTAGCCCTCACCATCCTCCTCAATGGCATCCAACCTGTGCTATCCGGAGTTGCCGTTGGTTCAGGGTGGCAAGCGCTAGTTGCGTATGTGAACATTGGGAGCTACTACATTATCGGTGTTCCTTTCGGTGTTCTGCTAGCATGGGGTTTCCACTACGGGGTCCTTGGCATTTGGGTTGGAATGATCGGTGGCACGATGGTGCAAACTCTGATTCTTTCATTTATCACCTTACGATGCGACTGGAATGAAGAGGCACTGAAAGCTTCTAGCAGAATGCGGACATGGAGCAGCTCCAAG ->species2_6 -ATGGAGGAGAATCGGAGCGATATCCCGCTCATCTCCGGCTCCGAGCTGCCGGACAGGAGGGGAGGAGGCAAGATCTCCGAGCTTGCGAAGGAGGTATGGGGAGAGTCCAAGAAGCTGTGGGTGGTCGCCGGCCCGGCCGCGTTCACGAGGCTGACATTCTATGGCATGACCGTGGTCAGCCAGGCCTTTGCCGGGCACATCGGTGACCTCGAGCTCGCCGCCTTCTCCATAGCCACCACCGTCATTTCTGGTCTCAGCTTTGGCTTCTTTGTTGGCATGGCGAGTGCAATGGAGACGCTGTGCGGCCAAGCCTACGGTGCAAAGCAGTACCACATGATGGGCATCTACCTGCAGCGCTCGTGGCTCATCCTCCTCAGCTTCGCCGTGCTTCTTACTCCGACCTACATCTTCAGCGAGCAGCTGCTCACCGCGCTGGGCCAGCCCGCCGAGCTGTCGCGCCAGGCGGGCTTGGTCAGCCTGTACATGCTCCCGCTGCACTTCGTCTACGCCATCGTCCTGCCGCTCAACAAGTTCCTGCAGTGCCAGCGCAAGAACTGGGTCGCCGCGGTCACCACGGCCGCGGCGTTCCCCGTTCACGTCGTCGCCACCTGGCTGCTGGTGCGTTGCTTCCGGCTCGGGGTCTTTGGAGCAGCGATGGCGCTCACCCTGTCCTGGGCACTCGCCACGGTGGGTCTCCTCTCGTATGCCTTGGGCGGCGGGTGCCCGGAGACGTGGAGGGGATTCTCAGCTTCTGCCTTCGTGGACTTGAAGGACTTCATCAAGTTGTCCGCGGCGTCTGGTGTCATGCTCTGCTTGGAGAATTGGTACTACCGGATCTTGGTTTTCCTGACGGGCTATGTGAAGAACGCTGAACTGGCTGTCGATGCACTGTCCATCTGTATAAGTTATGCTGGATGGGAGATGATGATTCATTTGGGATTCTTAGCAGGCACTGGGGTGAGGGTGGCTAATGAGCTCGGTGCAGCCAACGGAGCACGAGCGAGATTTGCGACAATTGTGTCGATGACGACATCATTTCTGATCAGCCTATTCATTAGTTTGCTCATCCTGATTTTCCATGACAAACTCGGAATGATCTTCTCGTCGAGTCAGGCTGTGATTGATGCAGTAGACAACATTTCCTTTCTGCTGGCCCTCACCATCCTCCTCAACGGAATCCAACCTGTGCTCTCTGGAGTTGCTGTTGGCTCAGGGTGGCAGGCATTGGTTGCTTATGTCAACATTGGGAGCTATTACTTGATTGGTGTTCCTTTCGGTTTTCTGCTAGGATGGGGCTTGCATTATGGGGTTCAAGGAATTTGGGTCGGAATGATCGTTGGCACAATGGTGCAAACTCTAATACTGGCATATATCACTCTACGGTGTGATTGGAATGAAGAGGCATTGAAAGCTAGTACCCGAATGCGGAGATGGAGCAACTCCAAG ->species2_9 -ATGGGCACACTAGGCGGGCACGTCGCGCCGGGCGCCTTCTTCTTCCTCATCGGCCTGTGGCATCTGTTCGGCCACAGCCGCCTGTTCTTGCTACAGCGGGGCTCCTACGTGGCTCCGGTGTGGTTCCCGGTGCCGGGCGTCCGTCACATCGAGCTCATAATGATAATAATCGGCTCGGTGATCTCCGTCTCGATGGAGCTCGTCATCGTGCAGCCGAAGCACCAGCCGTTCGACGACGACGGCACCATCCCCAGCGTCCACCTGCACAACTTCGAGCACGCGTCCATCTCGCTGGCGTGGCTCGTCTTCGCCGCCGCCACCATCCACATGGACAGGGTCCGGGCGCCGATGCGGGACGCGGTGTCGCAGCTGGCGGCCGCGGCCGCGTTCGCGCAGCAGCTGCTCATCTTCCACTTCCACTCCGCGGACCACGCGGGCGTGCAGGGGCGGTACCACCGTCTGCTGGAGATGGTGGTCGCCGTCACGCTCGCCGCCTCGCTGCTCTTGATCCCCTACCAACGGAGCATCGCGCTGAGCCTGGTCCGCTCGGCCAGCCTCGTGTTCCAGGGCGTCTGGTTCACCGTCATGGGCGTCATGATGTGGACGCCGGCGCTCGTCCCCAAAGGCTGCTTCATGAACGACGAAGATGGCCTCCAAGTCGTCCGGTGCCGCACCGACGAGGCGCTCGACCGCGCCAAGTCGCTCGTCAACCTGCAGTTCAACTGGTACCTGACCGGCACCGTGGCGTTCGTCGTCGTGTTCTACCTCCAGATGGCCAAGCAGTACCAGGAGCAGCCGCAGTACGCTCCGCTGGTGAAGGGAGGGAGAGGCAGCGATGGCCGGTGCACCATCGGAGAGGTCAATGACGACGAGGATGACCTTGAGGCCTCCAAAGGAGGCTTAGGATATATCGAAATTGAGAGG ->species2_10 -ATGGGCACTCTCGTCGGGCACGTCGCGCCGGGCGCCGGCTTCCTCCTCATCGGCCTGTGGCAGCTATTCAGCCACATCCGCCTGTTCCTGCTGCGCCCGAGCTCGTACTCTGCTCCGGTCTGGTTCCCGGCGCCGGGCGTGCGCCACCTCGAGCTCATACTCATCATCATCGGCGCGGCGATGTCCATCCTGATGGAGCTCGTCATCGGCCCCGCGAAGCACCAGCCGTTCGACGACGACGGCACCATCCCGTCAGACCACCTCCACAACTTCGAGCACGCGTCCATCTCGCTGGCGCTGCTCGTCTTCGCCGCGGTCACCATCCACCTCGACAGGGTAAAGGCGCCCCTGCGTGACGCCGTGTCGCAGCTCGTCGCCGCCGCGGCGTTCGCGCAGCAGCTGCTCATCTTCCACCTCCACTCGGCGGACCACATGGGCGTGGAGGGGCAGTTCCACTGGCTGCTGCAGACGGTCATCGCCGTCACGCTCGCCACCACGCTGCTCGGGATCCCTTACCCGCGGAGCATCGTGGTGAGCCTTGTCCGGTCGGCCAGCCTCGTGCTCCAGGGCGTCTGGTTCGTCGTCATGGGCGTCATGCTGTGGACGCCGGCGCTCATACCCAAGGGCTGCTTCCTCAACCTCGAGGAAGGGCACGACGTCGTCCGGTGCCGCACCGACGAGGCGCTCGACCGCGCCAAGTCGCTCGTCAACCTGCAGTTCAGCTGGTACCTCACCGGCACGGTGGTGTTCGTCGTCCTGTTCTACCTCCAGATGGCGAAGCTCTACCCCGAGGAGCCGCAGTATTTGCCGCTGGTGAAGGGAGGAGGCGGCGGCGGCGATGACCGCGATAGCCGGTTCAGCATCGGAGACGATGATCACGACGATGAGGACGATGTCGAGGCTGCAAAACGTGGCTTCGGACACGTGGTTAGCGGCACAAAGCCTGTCGAAATCGAGAGG