Introduction

This document describes the output produced by the pipeline.

Most of the plots are taken from the MultiQC report, which summarises results at the end of the pipeline.

Pipeline overview

The pipeline is built using Nextflow and predicts protein structures using the following methods:

See main README.md for a condensed overview of the steps in the pipeline, and the bioinformatics tools used at each step.

The directories listed below will be created in the output directory after the pipeline has finished. All paths are relative to the top-level results directory.

AlphaFold2

Output files
  • AlphaFold2/
    • <SEQUENCE NAME>/ that contains the computed MSAs, unrelaxed structures, relaxed structures, ranked structures, raw model outputs, prediction metadata, and section timings
    • <SEQUENCE NAME>.alphafold.pdb that is the structure with the highest pLDDT score (ranked first)
    • <SEQUENCE NAME>_plddt_mqc.tsv that presents the pLDDT scores per residue for each of the 5 predicted models
  • DBs/ that contains symbolic links to the downloaded database and parameter files

Below you can find an indicative example of the TSV file with the pLDDT scores per residue for each of the 5 predicted models produced by AlphaFold2, which is included in the MultiQC report:

Positionsrank_0rank_1rank_2rank_3rank_4
166.1760.6160.3264.2065.31
278.0174.2073.1177.3678.46
382.1678.1676.7080.2080.68
486.0382.7881.8882.1983.93
588.0884.3884.7385.5887.70
689.3786.0686.3186.8488.52
791.2788.2788.0987.0188.67
891.2889.4290.1787.4790.07
993.1090.0992.8690.7093.41
1093.2391.4293.0790.1392.91
1194.1292.4493.0089.9092.97
1295.1593.6394.2592.6694.38
1395.0993.7594.3692.5494.95
1494.0892.7293.4390.3193.63
1594.3493.7793.3191.7293.57
1695.5694.6294.4693.5595.20
1795.5494.7594.6593.6195.37
1893.9193.8993.3091.3392.95
1995.4895.7894.4893.9595.05
2095.9695.4695.1494.0195.83
2194.0694.0693.1391.6993.54
2292.9893.2891.1488.8091.25
2395.2895.1393.3991.4893.56
2493.4193.3892.3289.8592.40
2590.8891.4088.6085.6787.65
2689.3088.9084.5883.1184.52
2791.9690.9589.0486.4287.77
2891.2090.6888.7186.4387.62
2988.0187.5385.8383.1184.95
3081.2983.7277.7575.7674.84
3187.1486.9282.1082.3278.74
3292.3490.1389.0488.3186.49
3391.7088.9485.5285.9481.75
3490.1188.2384.3385.4780.01
3593.3591.4990.6089.4087.10
3694.1592.4790.1790.4886.77
3793.4092.0186.3887.8480.11
3892.7989.9789.3188.5585.15
3994.6691.2992.7490.6790.30
4095.9893.5894.3091.6990.73
4194.9492.5788.3188.4080.33
4292.8991.0384.0385.3174.66
4394.5493.4486.5084.9176.68
4496.9395.2392.4291.9886.11
4594.4092.2787.4089.0279.44
4691.7490.9481.3584.8874.93
4796.1994.4690.5189.8284.51
4894.8493.0491.0291.5787.72
4991.2389.3486.1087.6382.12
5091.6489.5884.9385.8879.38

ColabFold

Output files
  • colabfold/webserver/ or colabfold/local/ based on the selected mode that contains the computed MSAs, unrelaxed structures, relaxed structures, ranked structures, raw model outputs and scores, prediction metadata, logs and section timings
  • DBs/ that contains symbolic links to the downloaded database and parameter files

Below you can find some indicative examples of the output images produced by ColabFold, which are included in the MultiQC report:

Sequence coverage

Alt text

predicted Local Distance Difference Test (pLDDT)

Alt text

Predicted Aligned Error (PAE)

Alt text

MultiQC report

Output files
  • multiqc
    • multiqc_report.html: A standalone HTML file that can be viewed in your web browser.
    • multiqc_data/: Directory containing parsed statistics from the different tools used in the pipeline.
    • multiqc_plots/: Directory containing static images from the report in various formats.

MultiQC is a visualisation tool that generates a single HTML report summarising all samples in your project. Most of the pipeline QC results are visualised in the report and further statistics are available within the report data directory.

Results generated by MultiQC collate pipeline QC from AlphaFold2 or ColabFold.

The pipeline has special steps which also allow the software versions to be reported in the MultiQC output for future traceability. For more information about how to use MultiQC reports, see http://multiqc.info.

Pipeline information

Output files
  • pipeline_info/
    • Reports generated by Nextflow: execution_report.html, execution_timeline.html, execution_trace.txt and pipeline_dag.dot/pipeline_dag.svg.
    • Reports generated by the pipeline: pipeline_report.html, pipeline_report.txt and software_versions.yml. The pipeline_report* files will only be present if the --email / --email_on_fail parameter’s are used when running the pipeline.
    • Reformatted samplesheet files used as input to the pipeline: samplesheet.valid.csv.

Nextflow provides excellent functionality for generating various reports relevant to the running and execution of the pipeline. This will allow you to troubleshoot errors with the running of the pipeline, and also provide you with other information such as launch commands, run times and resource usage.