This is a python snakemake file which creates feature counts table for miRNA sequences combining cutadapt function and bowtie2 reference alignment.
- python/py37-anaconda-2019.10 2) snakemake/5.7.1-py37 3) gbc-cutadapt/2.10 4) gbc-bowtie2/2.4.1 5) subread/2.0.1
-
Activate the python anaconda-snakemake environment. conda activate snakemake
-
Edit the config.json file according to your data and file paths.
-
Ensure that meta-table contains all the necessary fields.
** NOTE EXACT HEADERS HAVE TO BE ENFORCED WITH ONLY ONE SPACE BETWEEN THE ENTERIES or key errors will be thrown during processing**
-
Try the dry run of the snakefile created to check any potential errors. snakemake -n -s Snakefile
-
Launch all the jobs.
snakemake --latency-wait 120 -p -j 100 -s Snakefile
- The pipeline will produce a countfinal.txt file containing all the details of the miRNA sequences required along with an output directory with all the trimmed and mapped reads.