--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml	Wed Mar 11 19:19:24 2026 +0000
@@ -0,0 +1,25 @@
+<macros>
+    <token name="@TOOL_VERSION@">3.2.0</token>
+    <token name="@VERSION_SUFFIX@">0</token>
+    <token name="@PROFILE@">25.1</token>
+
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@TOOL_VERSION@">preseq</requirement>
+        </requirements>
+    </xml>
+
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1038/nmeth.2375</citation>
+        </citations>
+    </xml>
+    
+    <xml name="macro_input_bam">
+        <param name="input_bam" type="data" format="bam" label="Input BAM file" help="IMPORTANT: This BAM file must be sorted by coordinate!"/>
+    </xml>
+
+    <xml name="macro_verbose">
+        <param argument="-v" name="verbose" type="boolean" truevalue="-v" falsevalue="" checked="false" label="Verbose Output" help="Print more information to the tool standard error."/>
+    </xml>
+</macros>
\ No newline at end of file

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/preseq_lc_extrap.xml	Wed Mar 11 19:19:24 2026 +0000
@@ -0,0 +1,58 @@
+<tool id="preseq_lc_extrap" name="Preseq lc_extrap" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
+    <description>Yield Extrapolation</description>
+    
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+
+    <command detect_errors="exit_code"><![CDATA[
+        ## 1. File Preparation: Symlink the input BAM to a standard name
+		ln -s '$input_bam' input.bam &&
+        
+		## 2. lc_extrap Analysis Execute
+        preseq lc_extrap 
+            -B input.bam
+            $verbose
+            -e $extrap_limit
+            -s $step_size
+            -o '$output_tsv' 
+    ]]></command>
+
+    <inputs>
+        <expand macro="macro_input_bam" />
+        <param argument="-e" name="extrap_limit" type="integer" value="10000000" label="Extrapolation Limit" help="Total number of reads to extrapolate to."/>
+        <param argument="-s" name="step_size" type="integer" value="100000" label="Step Size" help="Step size for extrapolation."/>
+        <expand macro="macro_verbose" />
+    </inputs>
+
+    <outputs>
+        <data name="output_tsv" format="tsv"/>
+    </outputs>
+
+    <tests>
+        <test expect_num_outputs="1">
+            <param name="input_bam" value="test.bam" ftype="bam"/>
+            <param name="extrap_limit" value="100000"/>
+            <param name="step_size" value="10000"/>
+            <output name="output_tsv" file="test_lc_extrap.tsv" ftype="tsv" lines_diff="4"/>
+        </test>
+    </tests>
+
+    <help><![CDATA[
+**Preseq: lc_extrap**
+
+The `lc_extrap` tool from the Preseq package predicts the expected yield of distinct reads for future sequencing efforts. By analyzing the current complexity of a genomic library, it calculates an extrapolation curve. This helps answer the question: "If I sequence this library deeper, how many new, unique reads will I discover?"
+
+**Input**
+
+- A **coordinate-sorted** BAM file containing aligned reads.
+
+**Output**
+
+- A tabular (TSV) file containing the extrapolated yield curve. This file is suitable for plotting and can be integrated into MultiQC reports.
+    ]]></help>
+    
+    <expand macro="citations"/>
+
+</tool>

Binary file test-data/test.bam has changed

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_c_curve.tsv	Wed Mar 11 19:19:24 2026 +0000
@@ -0,0 +1,3 @@
+total_reads	distinct_reads
+0	0
+1000	794

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_lc_extrap.tsv	Wed Mar 11 19:19:24 2026 +0000
@@ -0,0 +1,11 @@
+TOTAL_READS	EXPECTED_DISTINCT	LOWER_0.95CI	UPPER_0.95CI
+0	0	0	0
+10000.0	4419.7	3867.5	4892.9
+20000.0	6306.1	5185.3	7363.1
+30000.0	7389.4	5798.2	8944.7
+40000.0	8097.4	6178.2	10055.6
+50000.0	8587.5	6429.4	10874.3
+60000.0	8955.8	6609.1	11503.0
+70000.0	9233.2	6743.9	12001.0
+80000.0	9450.5	6848.8	12405.2
+90000.0	9627.0	6932.8	12739.8