Homec4science
Diffusion eSCT pipeline interoperability (public) Browse

R2915/example_pipelines/bioinformatics_pipeline2/slurm-submit-script-WITHOUT-containers.sh
a0d2698cc758public

Last Change
Raw File

slurm-submit-script-WITHOUT-containers.sh
View Options

#!/bin/bash
#SBATCH --job-name=pipeline2-no-containers
#SBATCH --partition=shi
#SBATCH --time=02:00:00
#SBATCH --cpus-per-task=4
#SBATCH --mem=10G
#SBATCH --output=no-containers.log
module load Bowtie2/2.2.9-goolf-1.7.20
module load SAMtools/1.3.1-goolf-1.7.20
module load freebayes/1.1.0-goolf-1.7.20_39e5e4b
# define the number of cpus to use
# you can use env vars provided by your scheduler or hardcode the number of cpus to use
export NUM_CPUS=${SLURM_CPUS_PER_TASK}
##########################################################
# YOU SHOULD NOT NEED TO CHANGE ANYTHING BELOW THIS LINE
##########################################################
# create output directories if they don't exist
[ -d output/bowtie2 ] || mkdir -p output/bowtie2
[ -d output/samtools ] || mkdir -p output/samtools
[ -d output/freebayes ] || mkdir -p output/freebayes
RefGenome='./scicore-pipeline2-input-data/RefGenome/igenomes/Saccharomyces_cerevisiae/Ensembl/R64-1-1/Sequence/'
echo 'aligning reads with bowtie2...'
bowtie2 -p ${NUM_CPUS} -x $RefGenome/Bowtie2Index/genome -q -1 scicore-pipeline2-input-data/fastq/SRR5511076_1.ds.fastq -2 scicore-pipeline2-input-data/fastq/SRR5511076_2.ds.fastq -S output/bowtie2/yeast_reseq_ds-bt2aln.sam
echo 'SAM to sorted BAM...'
samtools view -bS output/bowtie2/yeast_reseq_ds-bt2aln.sam | samtools sort - --threads ${NUM_CPUS} -o output/samtools/yeast_reseq_ds-bt2aln.s.bam
#rm output/bowtie2/yeast_reseq_ds-bt2aln.sam
echo 'calling variants with FreeBayes...'
freebayes -f $RefGenome/WholeGenomeFasta/genome.fa output/samtools/yeast_reseq_ds-bt2aln.s.bam > output/freebayes/yeast_reseq_ds-bt2aln.s.vcf
c4science · Help