Skip to content
Snippets Groups Projects

Resolve "get unmapped reads"

Merged Kenzo-Hugo Hillion requested to merge 32-unmapped-reads into master
6 files
+ 124
0
Compare changes
  • Side-by-side
  • Inline
Files
6
 
configfile: "config.yaml"
 
 
# ==== Snakefile path ====
 
__bedtools_bamtofastq_pairs_rules = config.get("snakefiles", {}).get("bedtools_bamtofastq_pairs")
 
 
__main_output_dir = config.get('output_dir', 'output')
 
 
# ==== Main config ====
 
SAMPLES = config.get('samples')
 
__input_dir = config.get('input_dir', 'data')
 
 
# ==== Bedtools bamtofastq ====
 
__bedtools_bamtofastq_pairs_output_dir = f"{__main_output_dir}/bedtools_bamtofastq_pairs"
 
__bedtools_bamtofastq_pairs_input = "{dir}/{{sample}}.bam".format(dir=__input_dir, sample="{sample}")
 
__bedtools_bamtofastq_pairs_output_forward_name = "{dir}/{{sample}}_1".format(
 
dir=__bedtools_bamtofastq_pairs_output_dir, sample="{sample}"
 
)
 
__bedtools_bamtofastq_pairs_output_reverse_name = "{dir}/{{sample}}_2".format(
 
dir=__bedtools_bamtofastq_pairs_output_dir, sample="{sample}"
 
)
 
include: __bedtools_bamtofastq_pairs_rules
 
 
rule all:
 
input:
 
forward = expand(f"{__bedtools_bamtofastq_pairs_output_forward_name}.fq.gz", sample=SAMPLES),
 
reverse = expand(f"{__bedtools_bamtofastq_pairs_output_reverse_name}.fq.gz", sample=SAMPLES)
Loading