add rasusa wrapper

bio/rasusa/environment.yaml

@@ -0,0 +1,5 @@
+channels:
+  - bioconda
+  - conda-forge
+dependencies:
+  - rasusa ==0.3.0

bio/rasusa/meta.yaml

@@ -0,0 +1,4 @@
+name: rasusa
+description: Randomly subsample sequencing reads to a specified coverage using `rasusa <https://github.com/mbhall88/rasusa>`_.
+authors:
+  - Michael Hall

bio/rasusa/test/Snakefile

@@ -0,0 +1,15 @@
+rule subsample:
+    input:
+        r1="{sample}.r1.fq",
+        r2="{sample}.r2.fq",
+    output:
+        r1="{sample}.subsampled.r1.fq",
+        r2="{sample}.subsampled.r2.fq",
+    params:
+        options="--seed 15",
+        genome_size="3mb", # required
+        coverage=20, # required
+    log:
+        "logs/subsample/{sample}.log",
+    wrapper:
+        "master/bio/rasusa"

bio/rasusa/test/a.r1.fq

@@ -0,0 +1,4 @@
+@read1
+ACGT
++
+$$$$

bio/rasusa/test/a.r2.fq

@@ -0,0 +1,4 @@
+@read1
+ACGT
++
+!!!!

bio/rasusa/wrapper.py

@@ -0,0 +1,17 @@
+__author__ = "Michael Hall"
+__copyright__ = "Copyright 2020, Michael Hall"
+__email__ = "[email protected]"
+__license__ = "MIT"
+
+
+from snakemake.shell import shell
+
+
+options = snakemake.params.get("options", "")
+
+
+shell(
+    "rasusa {options} -i {snakemake.input} -o {snakemake.output} "
+    "-c {snakemake.params.coverage} -g {snakemake.params.genome_size} "
+    "2> {snakemake.log}"
+)