ASPICov: Automated and Standardized Pipeline for Identification of SARS-Cov2 nucleotidic Variants

public public 1yr ago Version: Version 1 0 bookmarks

ASPICov was created to provide biologists with a fast, reliable, and comprehensive analysis of NGS SARS-Cov2 samples. This versatile application tool enables the processing of samples from either the capture or amplicon method, as well as Illumina or Ion Torrent technology. This Nextflow pipeline adheres to nf-core principles and employs Singularity containers to ensure FAIR data analysis.

Code Snippets

17
18
19
"""
bamtools.sh ${bam_file}  ${sample_id}  bamtools.cmd  >& bamtools.log 2>&1
"""
19
20
21
"""
bcftools.sh ${vcf_file}  ${sample_id}  bcftools.cmd  >& bcftools.log 2>&1
"""
18
19
20
"""
bwa.sh  ${params.ref_fasta} bwa.cmd ${trimmed_files} ${sample_id} >& bwa.log 2>&1
"""
NextFlow From line 18 of modules/bwa.nf
43
44
45
"""
bwa.sh  ${params.ref_fasta} bwa2.cmd ${primerclipped_file} ${sample_id} >& bwa2.log 2>&1
"""
NextFlow From line 43 of modules/bwa.nf
16
17
18
"""
bgzip.sh ${vcf_annot_filter_pass_file}  ${sample_id}  bgzip.cmd  >& bgzip.log 2>&1
"""
37
38
39
"""
tabix.sh ${vcf_zip_file}  ${sample_id}  tabix.cmd  >& tabix.log 2>&1
"""
57
58
59
"""
consensus.sh ${vcf_zip_file}  ${sample_id}  ${params.ref_fasta}  consensus.cmd  >& consensus.log 2>&1
"""
81
82
83
"""
extraction.sh  ${consensus_file_nn}  ${sample_id}  extraction.cmd  >& extraction.log 2>&1
"""
15
16
17
"""
CoveragePlot.py ${wig_file} >& cvplot.log 2>&1
"""
19
20
21
"""
filter.sh  ${vcf_annot_filter_file} ${sample_id} filter.cmd  >& filter.log 2>&1
"""
16
17
18
"""
freebayes.sh ${params.ref_fasta}  ${sample_id} ${bam_ch}  freebayes.cmd  >& freebayes.log 2>&1
"""
 8
 9
10
"""
get_singularity_images.sh ${baseDir} singularity_images_ok >& get_singularity_images.log 2>&1
"""
14
15
16
"""
get_test_data.sh ${baseDir} test_data_ok  >& get_test_data.log 2>&1
"""
17
18
19
"""
igvtools.sh ${sample_id}  ${bam_ch} ${params.ref_genome}  igvtools.cmd  >& igvtools.log 2>&1
"""
17
18
19
"""
index.sh ${ref_fasta}  index.cmd  >& index.log 2>&1
"""
20
21
22
""" 
megahit.sh ${params.technology} ${sample_id} megahit.cmd ${trimmed_files.join(' ')} >& megahit.log 2>&1
"""
25
26
27
"""
megahit.sh ${params.technology} ${sample_id} megahit.cmd ${trimmed_files} >& megahit.log 2>&1
"""
18
19
20
21
"""
normalization.sh  ${params.ref_fasta}  ${sample_id}  ${vcf_tmp_file} normalization.cmd  >& normalization.log 2>&1

"""
19
20
21
"""
pangolin.sh "${consensus_file}"  pangolin.cmd  >& pangolin.log 2>&1
"""
21
22
23
"""
trimming-primers.py ${params.bedpe_file} ${samfile} ${params.ref_fasta} ${sample_id} >& primer_trimming.log 2>&1
"""
16
17
18
"""
prinseq.sh ${fastq_file}  ${sample_id}  prinseq.cmd  >& prinseq.log 2>&1
"""
15
16
17
"""
plots-variants.py -p . -o plot-variant.pdf >& pvariant.log 2>&1
"""
17
18
19
""" 
fastqc.sh ${params.technology}  ${sample_id} fastqc.cmd ${trimmed_files.join(' ')} >& fastqc.log 2>&1
"""
22
23
24
"""
fastqc.sh  ${params.technology} ${sample_id} fastqc.cmd ${trimmed_files}  >& fastqc.log 2>&1
"""
42
43
44
"""
multiqc.sh multiqc.cmd >& multiqc.log 2>&1
"""
22
23
24
"""
quast.sh ${denovo_file} ${sample_id} quast.cmd >& quast.log 2>&1
"""
29
30
31
"""
samtools.sh ${sample_id} ${sam_file} ${params.ref_fasta} ${params.ref_bed} samtools.cmd >& samtools.log 2>&1
"""
62
63
64
"""
samtools.sh ${sample_id} ${sam2_file} ${params.ref_fasta} ${params.ref_bed} samtools2.cmd >& samtools2.log 2>&1
"""
19
20
21
"""
snpeff.sh ${vcf_file} ${sample_id} snpeff.cmd  >& snpeff.log 2>&1
"""
17
18
19
"""
snpsift.sh ${vcf_annot_filter_pass_file}  ${sample_id}  snpsift.cmd  >& snpsift.log 2>&1
"""
40
41
42
"""
bcftools_extr.sh ${vcf_annot_filter_pass_file}  ${sample_id}  ${xls_fields_file}  bcftools_extr.cmd  >& bcftools_extr.log 2>&1
"""
16
17
18
"""
svh.sh . svh.cmd  >& svh.log 2>&1
"""
NextFlow From line 16 of modules/svh.nf
16
17
18
"""
synthesis1.sh . synthesis1.cmd >& synthesis1.log 2>&1
"""
37
38
39
"""
synthesis2.sh . synthesis2.cmd >& synthesis2.log 2>&1
"""
63
64
65
"""
synthesis3.sh . synthesis3.cmd >& synthesis3.log 2>&1
"""
18
19
20
"""
trimmomatic.sh ${task.cpus} ${reads[0]} ${reads[1]} ${sample_id} ${params.trimmomatic_adapter_file} trimmomatic.cmd >& trimmomatic.log 2>&1
"""
17
18
19
"""
vcfallelicprimitives.sh ${vcf_norm_file} ${sample_id} vcfallelicprimitives.cmd  >& vcfallelicprimitives.log 2>&1
"""
ShowHide 37 more snippets with no or duplicated tags.

Login to post a comment if you would like to share your experience with this workflow.

Do you know this workflow well? If so, you can request seller status , and start supporting this workflow.

Free

Created: 1yr ago
Updated: 1yr ago
Maitainers: public
URL: https://gitlab.com/vtilloy/aspicov
Name: aspicov
Version: Version 1
Badge:
workflow icon

Insert copied code into your website to add a link to this workflow.

Downloaded: 0
Copyright: Public Domain
License: Apache License 2.0
  • Future updates

Related Workflows

cellranger-snakemake-gke
snakemake workflow to run cellranger on a given bucket using gke.
A Snakemake workflow for running cellranger on a given bucket using Google Kubernetes Engine. The usage of this workflow ...