The Identify ARTIC V3 SARS-CoV-2 Low Frequency and Shared Variants (Illumina) workflow includes all necessary steps for processing paired-end reads from SARS-CoV-2 samples, such as sample QC, trimming of adapters (note that it might be necessary to adjust trim adapters depending on propocol) and primers, variant calling relative to reference MN908947.3 and extraction of a consensus sequence. Default the workflow is configured to trim ARTIC V3 primers and adapters of the reads, however if these steps are unnecessary they can be removed or adjusted.
The workflow can be found in the Toolbox at:
Template Workflows | Biomedical Workflows () | SARS-CoV-2 Workflows () | Identify ARTIC V3 SARS-CoV-2 Low Frequency and Shared Variants (Illumina) ()
If you are connected to a CLC Server via your Workbench, you will be asked where you would like to run the analysis. We recommend that you run the analysis on a CLC Server when possible.
In the next step, select the data to be analyzed. These can be sequence lists containing paired-end reads selected from the Navigation Area, or by using the "Select files for import" option, where files containing paired-end read data can be selected from disk. These will be imported as part of the workflow run. When importing Illumina paired-end data, the "Paired reads" option needs to be enabled.
The workflow contains an Iterate element, allowing each sample to be analyzed individually, before the results are combined for comparison. The "Batch" check box, at the bottom of the dialog, should normally remain unchecked when launching this workflow.
In the next step the relevant Reference Data Set is selected. The workflow uses SARS-CoV-2 reference MN908947.3 by default. Note that alternative reference data sets can be created, as described in http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Custom_Sets.html. Depending on protocol trim adapters adjustments might be necessary. Either extend the current data element or create a new element containing the correct trim adapter sequences and add it to a Custom Sets as described in the link above.
In the next step, you specify how the batch units are defined, that is, which data files come from each individual sample and thus should be analyzed together. Batch units can be defined through the organization of the input data or by using metadata. Further information can be found at http://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/current/index.php?manual=Running_workflows_in_batch_mode.html. When the metadata option is chosen, selecting an Excel file that describes the data will often be the most convenient method, and it is the only option available when input data will be imported as part of the workflow run. When using already-imported data as input, existing metadata tables, where associations from the input data are already in place, can also be selected.
In the next step, a preview of the batch units is shown. If this looks as expected, you can proceed to configure the analysis settings.
In the next step the Remove False Positives (high frequency) quality filter can be adjusted. High frequency variants (>=50%) with an average base quality lower than the specified value will be discarded (figure 4.2).
In the next step, the Remove False Positives (low frequency) quality filter applying filters for both quality and frequency. Lower frequency variants (Default >=10%) with an average base quality lower than the specified value will be discarded by default.
Finally, choose where to save the results.
The outputs generated are described in SARS-CoV-2 workflow output.