This pipeline starts with different individual assemblies (fasta files) and their respective reads (fastq.gz files). It will map all reads against all the assemblies using BWA.
Then per sample, Metabat2, Concoct and Metabinner are ran to create genomic bins. These bins are then consolidated using DasTool. After that, the bins of all samples are dereplicated with dRep to form MAGs. CheckM and GtdbTk are then used to estimate the quality and taxonomy of the MAGs
- Place your preprocessed/trim reads (e.g.
sample_R2.fastq.gzandsample_R2.fastq.gzfiles) in areadsfolder - Place the individual assemblies (e.g.
sample.fasta) into anassemblyfolder - Modify the
config.yamlfile to change the different paths and eventually the different options - Modify the
sample_list.txtfile to include your samples
snakemake -rp --cores 28 --use-conda --rerun-incomplete
This part was mainly taken from @susheelbhanu nomis_pipeline
- Modify the
slurm.yamlfile by checkingpartition,qosandaccountthat heavily depends on your system - Modify the
sbatch.shfile by checking#SBATCH -p,#SBATCH --qos=and#SBATCH -Aoptions that heavily depends on your system
sbatch config/sbatch.sh