COVID-19: variation analysis on ARTIC ONT data
This workflow for ONT-sequenced ARTIC data is modeled after the alignment/variant-calling steps of the ARTIC pipeline. It performs, essentially, the same steps as that pipeline’s minion command, i.e. read mapping with minimap2 and variant calling with medaka. Like the Illumina ARTIC workflow it uses ivar for primer trimming. Since ONT-sequenced reads have a much higher error rate than Illumina-sequenced reads and are therefor plagued more by false-positive variant calls, this workflow does make no attempt to handle amplicons affected by potential primer-binding site mutations.
Inputs
ID | Name | Description | Type |
---|---|---|---|
Maximum read length | Maximum read length | Discard reads that are longer than this length |
|
Minimum read length | Minimum read length | Discard reads that are shorter than this length |
|
NC_045512.2 FASTA sequence of SARS-CoV-2 | NC_045512.2 FASTA sequence of SARS-CoV-2 | Fasta sequence for Severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1, complete genome |
|
ONT-sequenced reads | ONT-sequenced reads | ONT reads from ARTIC assay with fastqsanger encoding |
|
Primer binding sites info in BED format | Primer binding sites info in BED format | BED file containing (ARTIC) primer positions. For ARTIC primer schemes the file can be obtained from https://usegalaxy.eu/u/wolfgang-maier/h/covid-19-resources. |
|
Steps
ID | Name | Description |
---|---|---|
5 | fastp | toolshed.g2.bx.psu.edu/repos/iuc/fastp/fastp/0.20.1+galaxy0 |
6 | Compute | toolshed.g2.bx.psu.edu/repos/devteam/column_maker/Add_a_column1/1.6 |
7 | Replace Text | toolshed.g2.bx.psu.edu/repos/bgruening/text_processing/tp_replace_in_column/1.1.3 |
8 | Map with minimap2 | toolshed.g2.bx.psu.edu/repos/iuc/minimap2/minimap2/2.17+galaxy2 |
9 | Datamash | toolshed.g2.bx.psu.edu/repos/iuc/datamash_ops/datamash_ops/1.1.0 |
10 | Samtools view | toolshed.g2.bx.psu.edu/repos/iuc/samtools_view/samtools_view/1.9+galaxy2 |
11 | Parse parameter value | param_value_from_file |
12 | Samtools stats | toolshed.g2.bx.psu.edu/repos/devteam/samtools_stats/samtools_stats/2.0.2+galaxy2 |
13 | BamLeftAlign | toolshed.g2.bx.psu.edu/repos/devteam/freebayes/bamleftalign/1.3.1 |
14 | ivar trim | toolshed.g2.bx.psu.edu/repos/iuc/ivar_trim/ivar_trim/1.3.1+galaxy2 |
15 | QualiMap BamQC | toolshed.g2.bx.psu.edu/repos/iuc/qualimap_bamqc/qualimap_bamqc/2.2.2d+galaxy3 |
16 | medaka consensus tool | toolshed.g2.bx.psu.edu/repos/iuc/medaka_consensus/medaka_consensus/1.0.3+galaxy2 |
17 | Filter failed | __FILTER_FAILED_DATASETS__ |
18 | medaka variant tool | toolshed.g2.bx.psu.edu/repos/iuc/medaka_variant/medaka_variant/1.3.2+galaxy1 |
19 | medaka variant tool | toolshed.g2.bx.psu.edu/repos/iuc/medaka_variant/medaka_variant/1.3.2+galaxy1 |
20 | Flatten Collection | __FLATTEN__ |
21 | bedtools Intersect intervals | toolshed.g2.bx.psu.edu/repos/iuc/bedtools/bedtools_intersectbed/2.30.0 |
22 | MultiQC | toolshed.g2.bx.psu.edu/repos/iuc/multiqc/multiqc/1.8+galaxy1 |
23 | bcftools annotate | toolshed.g2.bx.psu.edu/repos/iuc/bcftools_annotate/bcftools_annotate/1.10 |
24 | SnpEff eff covid19 version | toolshed.g2.bx.psu.edu/repos/iuc/snpeff_sars_cov_2/snpeff_sars_cov_2/4.5covid19 |
25 | Lofreq filter | toolshed.g2.bx.psu.edu/repos/iuc/lofreq_filter/lofreq_filter/2.1.5+galaxy0 |
26 | Replace | toolshed.g2.bx.psu.edu/repos/bgruening/text_processing/tp_find_and_replace/1.1.3 |
Outputs
ID | Name | Description | Type |
---|---|---|---|
length_filtered_reads | length_filtered_reads | n/a |
|
fastp_report | fastp_report | n/a |
|
minimap2_mapped_reads | minimap2_mapped_reads | n/a |
|
filtered_mapped_reads | filtered_mapped_reads | n/a |
|
mapped_reads_stats | mapped_reads_stats | n/a |
|
realigned_filtered_mapped_reads | realigned_filtered_mapped_reads | n/a |
|
primer_trimmed_realigned_filtered_mapped_reads | primer_trimmed_realigned_filtered_mapped_reads | n/a |
|
bamqc_raw_output | bamqc_raw_output | n/a |
|
bamqc_html_output | bamqc_html_output | n/a |
|
medaka_consensus_data | medaka_consensus_data | n/a |
|
medaka_variants_general | medaka_variants_general | n/a |
|
medaka_variants_pbs | medaka_variants_pbs | n/a |
|
bamqc_raw_output_flattened | bamqc_raw_output_flattened | n/a |
|
variants_pbs | variants_pbs | n/a |
|
preprocessing_and_mapping_reports | preprocessing_and_mapping_reports | n/a |
|
variants_combined | variants_combined | n/a |
|
annotated_variants | annotated_variants | n/a |
|
annotated_variants_stats | annotated_variants_stats | n/a |
|
lofreq_filtered | lofreq_filtered | n/a |
|
annotated_softfiltered_variants | annotated_softfiltered_variants | n/a |
|
Version History
v0.3.1 (latest) Created 21st Dec 2021 at 03:01 by WorkflowHub Bot
Updated to v0.3.1
Open
master
4948641
v0.1 (earliest) Created 12th Mar 2021 at 13:41 by WorkflowHub Bot
Added/updated 8 files
Frozen
master
3b99025
Creators
Not specifiedAdditional credit
Wolfgang Maier
Submitter
Views: 9265 Downloads: 1626 Runs: 8
Created: 12th Mar 2021 at 13:41
Last updated: 21st Dec 2021 at 03:01
None