#load package
devtools::load_all('/gpfs42/robbyfs/scratch/lab_lpasquali/msubirana/marc/repos/DNAseqPipeline')
#parse arguments
args <- base::commandArgs(trailingOnly = TRUE)
#define variables
file <- args[1]
out_dir <- args[2]
threads <- parallel::detectCores()
#realign BAMs using DNAseqPipeline::realignBwa
fastqc(file=file,
out_dir=out_dir,
threads=threads)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.