Add rule to create pre-miRNA annotation file

workflow/rules/prepare.smk

@@ -458,6 +458,35 @@ rule gfftobed:
         ) &> {log}"
 
 
+###############################################################################
+### Extract pre-miRNA
+###############################################################################
+
+
+rule filter_pre_mirs:
+    input:
+        bed=os.path.join(
+            config["output_dir"], "{organism}", "mirna_filtered.bed"
+        ),
+    output:
+        bed=os.path.join(
+            config["output_dir"], "{organism}", "pre_mirna_filtered.bed"
+        ),
+    params:
+        cluster_log=os.path.join(
+            config["cluster_log"], "{organism}", "filter_pre_mirs.log"
+        ),
+        precursor="miRNA_primary_transcript",
+    log:
+        os.path.join(
+            config["local_log"], "{organism}", "filter_pre_mirs.log"
+        ),
+    singularity:
+        "docker://zavolab/ubuntu:18.04"
+    shell:
+        "(grep {params.precursor} {input.bed} > {output.bed}) &> {log}"
+
+
 ###############################################################################
 ### Extract mature miRNA
 ###############################################################################
@@ -580,6 +609,7 @@ rule iso_anno:
         "(bedtools slop \
         -i {input.bed} \
         -g {input.chrsize} \
+        -s \
         -l {params.bp_5p} \
         -r {params.bp_3p} \
         > {output.bed} \