--- a/fasta-stats.xml	Thu Nov 22 04:16:35 2018 -0500
+++ b/fasta-stats.xml	Wed Apr 21 09:10:46 2021 +0000
@@ -1,4 +1,4 @@
-<tool id="fasta-stats" name="Fasta Statistics" version="1.0.1">
+<tool id="fasta-stats" name="Fasta Statistics" version="1.0.2">
     <description>Display summary statistics for a fasta file.</description>
     <requirements>
         <requirement type="package" version="5.26">perl</requirement>
@@ -6,11 +6,15 @@
     <command detect_errors="exit_code"><![CDATA[
         perl '${__tool_directory__}/fasta-stats.pl'
         '$dataset'
+        #if $genome_size:
+            $genome_size
+        #end if
         > '$stats'
         ]]>
     </command>
     <inputs>
         <param name="dataset" type="data" format="fasta" label="fasta or multifasta file" help="fasta dataset to get statistics for."/>
+        <param name="genome_size" type="float" optional="True" label="Genome size estimate (optional)" help="Estimate of the genome size in megabases (MB). If specified, NG50 will be calculated."/>
     </inputs>
     <outputs>
         <data name="stats" format="tabular" label="${tool.name} on ${on_string}: Fasta summary stats"/>
@@ -20,6 +24,11 @@
             <param name="dataset" value="test.fasta"/>
             <output name="stats" file="test_out.txt"/>
         </test>
+        <test>
+            <param name="dataset" value="test.fasta"/>
+            <param name="genome_size" value="5.0"/>
+            <output name="stats" file="ng50_out.txt"/>
+        </test>
     </tests>
     <help>
 **Fasta Stats**
@@ -36,6 +45,8 @@
 
     GC content in %
 
+    If an optional genome size estimate is specified, then the NG50 length will also be calculated. 
+
 ------
 
 Inputs: