diff --git a/CHANGELOG.md b/CHANGELOG.md index fb8159fa..5203d463 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -3,7 +3,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/) and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html). -## v2.2.0dev - [31-Oct-2024] +## v2.2.0dev - [04-Nov-2024] ### `Added` @@ -16,6 +16,9 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0 7. Added parameter `hic_samtools_ext_args` and set its default value to `-F 3852` [#159](https://github.com/Plant-Food-Research-Open/assemblyqc/issues/159) 8. Added the HiC QC report to the final report so that users don't have to navigate to the results folder [#162](https://github.com/Plant-Food-Research-Open/assemblyqc/issues/162) 9. Added the fastp log to the final report [#163](https://github.com/Plant-Food-Research-Open/assemblyqc/issues/163) +10. Updated the tube map along with the tool list [#166](https://github.com/Plant-Food-Research-Open/assemblyqc/issues/166) +11. Added Orthofinder [#167](https://github.com/Plant-Food-Research-Open/assemblyqc/issues/167) +12. Changed order of tool options in the `nextflow.config` file ### `Fixed` diff --git a/CITATIONS.md b/CITATIONS.md index 22ae2fb0..bf095d43 100644 --- a/CITATIONS.md +++ b/CITATIONS.md @@ -18,15 +18,15 @@ > Gremme G, Steinbiss S, Kurtz S. 2013. "GenomeTools: A Comprehensive Software Library for Efficient Processing of Structured Genome Annotations," in IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 10, no. 3, pp. 645-656, May 2013, doi: -- SAMTOOLS, [MIT/Expat](https://github.com/samtools/samtools/blob/develop/LICENSE) +- samtools, [MIT/Expat](https://github.com/samtools/samtools/blob/develop/LICENSE) > Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H. 2021. Twelve years of SAMtools and BCFtools, GigaScience, Volume 10, Issue 2, February 2021, giab008, -- NCBI/FCS, [License](https://github.com/ncbi/fcs/blob/main/LICENSE.txt) +- NCBI FCS, [License](https://github.com/ncbi/fcs/blob/main/LICENSE.txt) > Astashyn A, Tvedte ES, Sweeney D, Sapojnikov V, Bouk N, Joukov V, Mozes E, Strope PK, Sylla PM, Wagner L, Bidwell SL, Clark K, Davis EW, Smith-White B, Hlavina W, Pruitt KD, Schneider VA, Murphy TD. 2023. Rapid and sensitive detection of genome contamination at scale with FCS-GX. bioRxiv 2023.06.02.543519; doi: -- KRONA, [License](https://github.com/marbl/Krona/blob/master/KronaTools/LICENSE.txt) +- Krona, [License](https://github.com/marbl/Krona/blob/master/KronaTools/LICENSE.txt) > Ondov BD, Bergman NH, Phillippy AM. 2011. Interactive metagenomic visualization in a Web browser. BMC Bioinformatics. 2011 Sep 30;12:385. doi: @@ -36,7 +36,7 @@ > > Forked from: -- GFASTATS, [MIT](https://github.com/vgl-hub/gfastats/blob/main/LICENSE) +- gfastats, [MIT](https://github.com/vgl-hub/gfastats/blob/main/LICENSE) > Giulio Formenti, Linelle Abueg, Angelo Brajuka, Nadolina Brajuka, Cristóbal Gallardo-Alba, Alice Giani, Olivier Fedrigo, Erich D Jarvis, Gfastats: conversion, evaluation and manipulation of genome sequences using assembly graphs, Bioinformatics, Volume 38, Issue 17, September 2022, Pages 4214–4216, @@ -44,15 +44,15 @@ > Manni M, Berkeley MR, Seppey M, Simão FA, Zdobnov EM. 2021. BUSCO Update: Novel and Streamlined Workflows along with Broader and Deeper Phylogenetic Coverage for Scoring of Eukaryotic, Prokaryotic, and Viral Genomes, Molecular Biology and Evolution, Volume 38, Issue 10, October 2021, Pages 4647–4654, -- GFFREAD, [MIT](https://github.com/gpertea/gffread/blob/master/LICENSE) +- GffRead, [MIT](https://github.com/gpertea/gffread/blob/master/LICENSE) > Pertea G, Pertea M. GFF Utilities: GffRead and GffCompare. F1000Res. 2020 Apr 28;9:ISCB Comm J-304. doi: . PMID: 32489650; PMCID: PMC7222033. -- TIDK, [MIT](https://github.com/tolkit/telomeric-identifier/blob/main/LICENSE) +- tidk, [MIT](https://github.com/tolkit/telomeric-identifier/blob/main/LICENSE) > -- SEQKIT, [MIT](https://github.com/shenwei356/seqkit/blob/master/LICENSE) +- SeqKit, [MIT](https://github.com/shenwei356/seqkit/blob/master/LICENSE) > Shen W, Le S, Li Y, Hu F. 2016. SeqKit: A Cross-Platform and Ultrafast Toolkit for FASTA/Q File Manipulation. PLoS ONE 11(10): e0163962. @@ -72,19 +72,19 @@ > Shujun O, Ning J 2018. LTR_retriever: A Highly Accurate and Sensitive Program for Identification of Long Terminal Repeat Retrotransposons, Plant Physiology, 176, 2 (2018). -- KRAKEN2, [MIT](https://github.com/DerrickWood/kraken2/blob/master/LICENSE) +- Kraken 2, [MIT](https://github.com/DerrickWood/kraken2/blob/master/LICENSE) > Wood DE, Salzberg SL, Wood DE, Lu J, Langmead B. 2019. Improved metagenomic analysis with Kraken 2. Genome Biol 20, 257 (2019). -- JUICEBOX.JS, [MIT](https://github.com/igvteam/juicebox.js/blob/master/LICENSE) +- juicebox.js, [MIT](https://github.com/igvteam/juicebox.js/blob/master/LICENSE) > Robinson JT, Turner D, Durand NC, Thorvaldsdóttir H, Mesirov JP, Aiden EL. 2018. Juicebox.js Provides a Cloud-Based Visualization System for Hi-C Data. Cell Syst. 2018 Feb 28;6(2):256-258.e1. doi: . Epub 2018 Feb 7. PMID: 29428417; PMCID: PMC6047755. -- FASTP, [MIT](https://github.com/OpenGene/fastp/blob/master/LICENSE) +- fastp, [MIT](https://github.com/OpenGene/fastp/blob/master/LICENSE) > Chen S, Zhou Y, Chen Y, Gu J. 2018. fastp: an ultra-fast all-in-one FASTQ preprocessor, Bioinformatics, Volume 34, Issue 17, 01 September 2018, Pages i884–i890, -- FASTQC, [GPL v3](https://github.com/s-andrews/FastQC/blob/master/LICENSE.txt) +- FastQC, [GPL v3](https://github.com/s-andrews/FastQC/blob/master/LICENSE.txt) > @@ -92,7 +92,7 @@ > Dudchenko O, Batra SS, Omer AD, Nyquist SK, Hoeger M, Durand NC, Shamim MS, Machol I, Lander, Aiden AP, Aiden EL 2017. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds.Science356, 92-95(2017). doi: . Available at: -- HIC_QC, [AGPL v3](https://github.com/phasegenomics/hic_qc/blob/master/LICENSE) +- hic_qc, [AGPL v3](https://github.com/phasegenomics/hic_qc/blob/master/LICENSE) > @@ -100,42 +100,46 @@ > -- BWA, [GPL v3](https://github.com/lh3/bwa/blob/master/COPYING) +- bwa-mem, [GPL v3](https://github.com/lh3/bwa/blob/master/COPYING) > Li H. 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. -- MATLOCK, [AGPL v3](https://github.com/phasegenomics/matlock/blob/master/LICENSE) +- Matlock, [AGPL v3](https://github.com/phasegenomics/matlock/blob/master/LICENSE) > ; -- SAMBLASTER, [MIT](https://github.com/GregoryFaust/samblaster/blob/master/LICENSE.txt) +- samblaster, [MIT](https://github.com/GregoryFaust/samblaster/blob/master/LICENSE.txt) > Faust GG, Hall IM. 2014. SAMBLASTER: fast duplicate marking and structural variant read extraction, Bioinformatics, Volume 30, Issue 17, September 2014, Pages 2503–2505, -- CIRCOS, [GPL v3](https://www.gnu.org/licenses/gpl-3.0.txt) +- Circos, [GPL v3](https://www.gnu.org/licenses/gpl-3.0.txt) > Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R. Horsman D, ... Marra MA. 2009. Circos: an information aesthetic for comparative genomics. Genome research, 19(9), 1639-1645. -- MUMMER, [Artistic 2.0](https://github.com/mummer4/mummer/blob/master/LICENSE.md) +- MUMmer, [Artistic 2.0](https://github.com/mummer4/mummer/blob/master/LICENSE.md) > Marçais G, Delcher AL, Phillippy AM, Coston R, Salzberg SL, Zimin A. 2018. MUMmer4: A fast and versatile genome alignment system. PLoS Comput Biol. 2018 Jan 26;14(1):e1005944. doi: . PMID: 29373581; PMCID: PMC5802927. -- PLOTSR, [MIT](https://github.com/schneebergerlab/plotsr/blob/master/LICENSE) +- Plotsr, [MIT](https://github.com/schneebergerlab/plotsr/blob/master/LICENSE) > Goel M, Schneeberger K. 2022. plotsr: visualizing structural similarities and rearrangements between multiple genomes. Bioinformatics. 2022 May 13;38(10):2922-2926. doi: . PMID: 35561173; PMCID: PMC9113368. -- SYRI, [MIT](https://github.com/schneebergerlab/syri/blob/master/LICENSE) +- Syri, [MIT](https://github.com/schneebergerlab/syri/blob/master/LICENSE) > Goel M, Sun H, Jiao WB, Schneeberger K. 2019. SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol. 2019 Dec 16;20(1):277. doi: . PMID: 31842948; PMCID: PMC6913012. -- MINIMAP2, [MIT](https://github.com/lh3/minimap2/blob/master/LICENSE.txt) +- Minimap2, [MIT](https://github.com/lh3/minimap2/blob/master/LICENSE.txt) > Li H. 2021. New strategies to improve minimap2 alignment accuracy, Bioinformatics, Volume 37, Issue 23, December 2021, Pages 4572–4574, doi: -- MERQURY, [United States Government Work](https://github.com/marbl/merqury?tab=License-1-ov-file#readme) +- Merqury, [United States Government Work](https://github.com/marbl/merqury?tab=License-1-ov-file#readme) > Rhie, A., Walenz, B.P., Koren, S. et al. 2020. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol 21, 245. doi: +- OrthoFinder, [GPL v3](https://github.com/davidemms/OrthoFinder/blob/master/License.md) + + > Emms, D.M., Kelly, S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol 20, 238 (2019). doi: 10.1186/s13059-019-1832-y + ## Software packaging/containerisation tools - [Anaconda](https://anaconda.com) diff --git a/README.md b/README.md index 2283fc0b..0b4c71e7 100644 --- a/README.md +++ b/README.md @@ -12,26 +12,43 @@ ## Introduction -**plant-food-research-open/assemblyqc** is a [Nextflow](https://www.nextflow.io/docs/latest/index.html) pipeline which evaluates assembly quality with multiple QC tools and presents the results in a unified html report. The tools are shown in the [Pipeline Flowchart](#pipeline-flowchart) and their references are listed in [CITATIONS.md](./CITATIONS.md). +**plant-food-research-open/assemblyqc** is a [Nextflow](https://www.nextflow.io/docs/latest/index.html) pipeline which evaluates assembly quality with multiple QC tools and presents the results in a unified html report. The tools are shown in the [Pipeline Flowchart](#pipeline-flowchart) and their references are listed in [CITATIONS.md](./CITATIONS.md). The pipeline includes skip flags to disable execution of many of it tools. ## Pipeline Flowchart

-- [FASTA VALIDATOR](https://github.com/linsalrob/fasta_validator) + [SEQKIT RMDUP](https://github.com/shenwei356/seqkit): FASTA validation -- [GENOMETOOLS GT GFF3VALIDATOR](https://genometools.org/tools/gt_gff3validator.html): GFF3 validation -- [ASSEMBLATHON STATS](https://github.com/PlantandFoodResearch/assemblathon2-analysis/blob/a93cba25d847434f7eadc04e63b58c567c46a56d/assemblathon_stats.pl), [GFASTATS](https://github.com/vgl-hub/gfastats): Assembly statistics -- [GENOMETOOLS GT STAT](https://genometools.org/tools/gt_stat.html): Annotation statistics -- [NCBI FCS ADAPTOR](https://github.com/ncbi/fcs): Adaptor contamination pass/fail -- [NCBI FCS GX](https://github.com/ncbi/fcs): Foreign organism contamination pass/fail -- [BUSCO](https://gitlab.com/ezlab/busco): Gene-space completeness estimation -- [TIDK](https://github.com/tolkit/telomeric-identifier): Telomere repeat identification -- [LAI](https://github.com/oushujun/LTR_retriever/blob/master/LAI): Continuity of repetitive sequences -- [KRAKEN2](https://github.com/DerrickWood/kraken2): Taxonomy classification -- [HIC CONTACT MAP](https://github.com/igvteam/juicebox.js): Alignment and visualisation of HiC data -- [MUMMER](https://github.com/mummer4/mummer) → [CIRCOS](http://circos.ca/documentation/) + [DOTPLOT](https://plotly.com) & [MINIMAP2](https://github.com/lh3/minimap2) → [PLOTSR](https://github.com/schneebergerlab/plotsr): Synteny analysis -- [MERQURY](https://github.com/marbl/merqury): K-mer completeness, consensus quality and phasing assessment -- [ORTHOFINDER](https://github.com/davidemms/OrthoFinder): Phylogenetic orthology inference for comparative genomics +- `Assembly` + - [fasta_validator](https://github.com/linsalrob/fasta_validator) + [SeqKit rmdup](https://github.com/shenwei356/seqkit): FASTA validation + - [assemblathon_stats](https://github.com/PlantandFoodResearch/assemblathon2-analysis/blob/a93cba25d847434f7eadc04e63b58c567c46a56d/assemblathon_stats.pl), [gfastats](https://github.com/vgl-hub/gfastats): Assembly statistics + - [NCBI FCS-adaptor](https://github.com/ncbi/fcs): Adaptor contamination pass/fail + - [NCBI FCS-GX](https://github.com/ncbi/fcs): Foreign organism contamination pass/fail + - [tidk](https://github.com/tolkit/telomeric-identifier): Telomere repeat identification + - [BUSCO](https://gitlab.com/ezlab/busco): Gene-space completeness estimation + - [LAI](https://github.com/oushujun/LTR_retriever/blob/master/LAI): Continuity of repetitive sequences + - [Kraken 2](https://github.com/DerrickWood/kraken2), [Krona](https://github.com/marbl/Krona): Taxonomy classification + - `Alignment and visualisation of HiC data` + - [sra-tools](https://github.com/ncbi/sra-tools): HiC data download from SRA or use of local FASTQ files + - [fastp](https://github.com/OpenGene/fastp), [FastQC](https://github.com/s-andrews/FastQC): QC and trimming + - [SeqKit sort](https://github.com/shenwei356/seqkit): Alphanumeric sorting of FASTA by sequence ID + - [bwa-mem](https://github.com/lh3/bwa): HiC read alignment + - [samblaster](https://github.com/GregoryFaust/samblaster): Duplicate marking + - [hic_qc](https://github.com/phasegenomics/hic_qc): HiC read and alignment statistics + - [Matlock](https://github.com/phasegenomics/matlock): BAM to juicer conversion + - [3d-dna/visualize](https://github.com/aidenlab/3d-dna/tree/master/visualize): `.hic` file creation + - [juicebox.js](https://github.com/igvteam/juicebox.js): HiC contact map visualisation + - `K-mer completeness, consensus quality and phasing assessment` + - [sra-tools](https://github.com/ncbi/sra-tools): Assembly, maternal and paternal data download from SRA or use of local FASTQ files + - [Merqury hapmers](https://github.com/marbl/merqury/blob/master/trio/hapmers.sh): Hapmer generation if parental data is available + - [Merqury](https://github.com/marbl/merqury): Completeness, consensus quality and phasing assessment + - `Synteny analysis` + - [MUMmer](https://github.com/mummer4/mummer) → [Circos](http://circos.ca/documentation/) + [dotplot](https://plotly.com): One-to-all and all-to-all synteny analysis at the contig level + - [Minimap2](https://github.com/lh3/minimap2) → [Syri](https://github.com/schneebergerlab/syri)/[Plotsr](https://github.com/schneebergerlab/plotsr): One-one to chromosome synteny analysis at the chromosome level +- `Annotation` + - [GenomeTools gt gff3validator](https://genometools.org/tools/gt_gff3validator.html) + [FASTA/GFF correspondence](subworkflows/gallvp/gff3_gt_gff3_gff3validator_stat/main.nf): GFF3 validation + - [GenomeTools gt stat](https://genometools.org/tools/gt_stat.html): Annotation statistics + - [GffRead](https://github.com/gpertea/gffread), [BUSCO](https://gitlab.com/ezlab/busco): Gene-space completeness estimation in annotation proteins + - [OrthoFinder](https://github.com/davidemms/OrthoFinder): Phylogenetic orthology inference for comparative genomics ## Usage @@ -50,9 +67,9 @@ Now, you can run the pipeline using: ```bash nextflow run plant-food-research-open/assemblyqc \ -revision \ - -profile \ - --input assemblysheet.csv \ - --outdir + -profile \ + --input assemblysheet.csv \ + --outdir ``` > [!WARNING] diff --git a/docs/images/assemblyqc.png b/docs/images/assemblyqc.png index 54ef8b9d..3fd97dd1 100644 Binary files a/docs/images/assemblyqc.png and b/docs/images/assemblyqc.png differ diff --git a/docs/images/orthofinder.png b/docs/images/orthofinder.png new file mode 100644 index 00000000..4a53e3c0 Binary files /dev/null and b/docs/images/orthofinder.png differ diff --git a/docs/output.md b/docs/output.md index dd5dd998..68adc287 100644 --- a/docs/output.md +++ b/docs/output.md @@ -8,26 +8,27 @@ The directories listed below will be created in the results directory after the ## Pipeline overview -The pipeline is built using [Nextflow](https://www.nextflow.io/) and processes data using the following steps: +The pipeline is built using [Nextflow](https://www.nextflow.io/) and processes data to produce following outputs: -- [FASTA and GFF3 validation](#fasta-and-gff3-validation) +- [Format validation](#format-validation) - [Assemblathon stats](#assemblathon-stats) -- [Gfastats](#gfastats) -- [GenomeTools gt stat](#genometools-gt-stat) -- [NCBI FCS adaptor](#ncbi-fcs-adaptor) -- [NCBI FCS GX](#ncbi-fcs-gx) +- [gfastats](#gfastats) +- [NCBI FCS-adaptor](#ncbi-fcs-adaptor) +- [NCBI FCS-GX](#ncbi-fcs-gx) +- [tidk](#tidk) - [BUSCO](#busco) -- [TIDK](#tidk) - [LAI](#lai) -- [Kraken2](#kraken2) +- [Kraken 2](#kraken-2) - [HiC contact map](#hic-contact-map) -- [Synteny](#synteny) - [Merqury](#merqury) +- [Synteny](#synteny) +- [GenomeTools gt stat](#genometools-gt-stat) +- [OrthoFinder](#orthofinder) - [Pipeline information](#pipeline-information) -### FASTA and GFF3 validation +### Format validation The pipeline prints a warning in the pipeline log if FASTA or GFF3 validation fails. The error log from the validator is reported in the `report.html`. The remaining QC tools are skipped for the assembly with invalid fasta file. @@ -46,7 +47,7 @@ The pipeline prints a warning in the pipeline log if FASTA or GFF3 validation fa > [!WARNING] > Contig-related stats are based on the assumption that `assemblathon_stats_n_limit` is specified correctly. If you are not certain of the value of `assemblathon_stats_n_limit`, please ignore the contig-related stats. -### Gfastats +### gfastats
Output files @@ -56,35 +57,21 @@ The pipeline prints a warning in the pipeline log if FASTA or GFF3 validation fa
-Gfastats is a fast and exhaustive tool for summary statistics. +gfastats is a fast and exhaustive tool for summary statistics. -### GenomeTools gt stat - -
-Output files - -- `genometools_gt_stat/` - - `*.gt.stat.yml`: Assembly annotation stats in yaml format. - -
- -GenomeTools `gt stat` tool calculates a basic set of statistics about features contained in GFF3 files. - -
AssemblyQC - GenomeTools gt stat gene length distribution
AssemblyQC - GenomeTools gt stat gene length distribution
- -### NCBI FCS adaptor +### NCBI FCS-adaptor
Output files - `ncbi_fcs_adaptor/` - - `*_fcs_adaptor_report.tsv`: NCBI FCS adaptor report in CSV format. + - `*_fcs_adaptor_report.tsv`: NCBI FCS-adaptor report in CSV format.
-[FCS-adaptor detects](https://github.com/ncbi/fcs/wiki/FCS-adaptor#rules-for-action-assignment) adaptor and vector contamination in genome sequences. +[FCS-adaptor](https://github.com/ncbi/fcs/wiki/FCS-adaptor#rules-for-action-assignment) detects adaptor and vector contamination in genome sequences. -### NCBI FCS GX +### NCBI FCS-GX
Output files @@ -98,42 +85,47 @@ GenomeTools `gt stat` tool calculates a basic set of statistics about features c
-[FCS-GX detects](https://github.com/ncbi/fcs/wiki/FCS-GX#outputs) contamination from foreign organisms in genome sequences. +[FCS-GX](https://github.com/ncbi/fcs/wiki/FCS-GX#outputs) detects contamination from foreign organisms in genome sequences. -### BUSCO +### tidk
Output files -- `busco/` - - `busco_figure.png`: Summary figure created from all the BUSCO summaries. - - `tag` - - `short_summary.specific.*_odb10.tag_*.txt`: BUSCO summary for the assembly represented by `tag`. +- `tidk/` + - `*.apriori.tsv`: Frequencies for successive windows in forward and reverse directions for the pre-specified telomere-repeat sequence. + - `*.apriori.svg`: Plot of `*.apriori.tsv` + - `*.tidk.explore.tsv`: List of the most frequent repeat sequences. + - `*.top.sequence.txt`: The top sequence from `*.tidk.explore.tsv`. + - `*.aposteriori.tsv`: Frequencies for successive windows in forward and reverse directions for the top sequence from `*.top.sequence.txt`. + - `*.aposteriori.svg`: Plot of `*.aposteriori.tsv`.
-[BUSCO estimates](https://busco.ezlab.org/busco_userguide.html) the completeness and redundancy of processed genomic data based on universal single-copy orthologs. +tidk toolkit is designed to [identify and visualize](https://github.com/tolkit/telomeric-identifier) telomeric repeats for the Darwin Tree of Life genomes. -
AssemblyQC - BUSCO summary plot
AssemblyQC - BUSCO summary plot
+
AssemblyQC - tidk plot
AssemblyQC - tidk plot
-### TIDK +### BUSCO
Output files -- `tidk/` - - `*.apriori.tsv`: Frequencies for successive windows in forward and reverse directions for the pre-specified telomere-repeat sequence. - - `*.apriori.svg`: Plot of `*.apriori.tsv` - - `*.tidk.explore.tsv`: List of the most frequent repeat sequences. - - `*.top.sequence.txt`: The top sequence from `*.tidk.explore.tsv`. - - `*.aposteriori.tsv`: Frequencies for successive windows in forward and reverse directions for the top sequence from `*.top.sequence.txt`. - - `*.aposteriori.svg`: Plot of `*.aposteriori.tsv`. +- `busco/` + - `fasta` + - `busco_figure.png`: Summary figure created from all the BUSCO summaries. + - `tag` + - `short_summary.specific.*_odb10.tag_*.txt`: BUSCO summary for the assembly represented by `tag`. + - `gff` + - `busco_figure.png`: Summary figure created from all the BUSCO summaries. + - `tag` + - `short_summary.specific.*_odb10.tag_*.txt`: BUSCO summary for the annotation of the assembly represented by `tag`.
-TIDK toolkit is designed to [identify and visualize](https://github.com/tolkit/telomeric-identifier) telomeric repeats for the Darwin Tree of Life genomes. +[BUSCO](https://busco.ezlab.org/busco_userguide.html) estimates the completeness and redundancy of processed genomic data based on universal single-copy orthologs. -
AssemblyQC - TIDK plot
AssemblyQC - TIDK plot
+
AssemblyQC - BUSCO summary plot
AssemblyQC - BUSCO summary plot
### LAI @@ -154,22 +146,22 @@ LTR Assembly Index (LAI) is a reference-free genome metric that [evaluates assem > [!WARNING] > Soft masked regions are unmasked when calculating LAI. However, hard masked regions are left as is. The pipeline will fail to calculate LAI if all the LTRs are already hard masked. -### Kraken2 +### Kraken 2
Output files - `kraken2/` - - `*.kraken2.report`: [Kraken2 report](https://github.com/DerrickWood/kraken2/wiki/Manual#output-formats). - - `*.kraken2.cut`: [Kraken2 output](https://github.com/DerrickWood/kraken2/wiki/Manual#output-formats). + - `*.kraken2.report`: [Kraken 2 report](https://github.com/DerrickWood/kraken2/wiki/Manual#output-formats). + - `*.kraken2.cut`: [Kraken 2 output](https://github.com/DerrickWood/kraken2/wiki/Manual#output-formats). - `*.kraken2.krona.cut`: [Select columns](../modules/local/kraken2_krona_plot.nf) from `*.kraken2.cut` used for generation of a Krona taxonomy plot. - `*.kraken2.krona.html`: Interactive Krona taxonomy plot.
-Kraken2 [assigns taxonomic labels](https://ccb.jhu.edu/software/kraken2/) to sequencing reads for metagenomics projects. Further reading regarding performance of Kraken2: +Kraken 2 [assigns taxonomic labels](https://ccb.jhu.edu/software/kraken2/) to sequencing reads for metagenomics projects. Further reading regarding performance of Kraken 2: -
AssemblyQC - Interactive Krona plot from Kraken2 taxonomy
AssemblyQC - Interactive Krona plot from Kraken2 taxonomy
+
AssemblyQC - Interactive Krona plot from Kraken 2 taxonomy
AssemblyQC - Interactive Krona plot from Kraken 2 taxonomy
### HiC contact map @@ -178,17 +170,17 @@ Kraken2 [assigns taxonomic labels](https://ccb.jhu.edu/software/kraken2/) to seq - `hic/` - `fastqc_raw/` - - `*_1_fastqc.html/*_2_fastqc.html`: FASTQC html report for the raw reads - - `*_1_fastqc.zip/*_2_fastqc.zip`: FASTQC stats for the raw reads + - `*_1_fastqc.html/*_2_fastqc.html`: FastQC html report for the raw reads + - `*_1_fastqc.zip/*_2_fastqc.zip`: FastQC stats for the raw reads - `fastp/` - - `*.fastp.html`: FASTP HTML report - - `*.fastp.json`: FASTP statistics in JSON format - - `*.fastp.log`: FASTP log - - `*_1.fastp.fastq.gz/*_2.fastp.fastq.gz`: Reads passed by FASTP - - `*_1.fail.fastq.gz/*_2.fail.fastq.gz`: Reads failed by FASTP + - `*.fastp.html`: fastp HTML report + - `*.fastp.json`: fastp statistics in JSON format + - `*.fastp.log`: fastp log + - `*_1.fastp.fastq.gz/*_2.fastp.fastq.gz`: Reads passed by fastp + - `*_1.fail.fastq.gz/*_2.fail.fastq.gz`: Reads failed by fastp - `fastqc_trim/` - - `*_1_fastqc.html/*_2_fastqc.html`: FASTQC html report for the reads passed by FASTP. - - `*_1_fastqc.zip/*_2_fastqc.zip`: FASTQC stats for the reads passed by FASTP. + - `*_1_fastqc.html/*_2_fastqc.html`: FastQC html report for the reads passed by FASTP. + - `*_1_fastqc.zip/*_2_fastqc.zip`: FastQC stats for the reads passed by FASTP. - `hicqc` - `*.on.*_qc_report.pdf`: HiC QC report for reads mapped to an assembly. - `assembly/` @@ -206,6 +198,28 @@ Hi-C contact mapping experiments measure the frequency of physical contact betwe AssemblyQC - HiC results +### Merqury + +
+Output files + +- `merqury/` + - `tag1-and-tag2`: Results folder for haplotype `tag1` and `tag2`. + - `*.completeness.stats`: Assembly completeness statistics + - `*.qv`: Assembly consensus quality QV statistics + - `*.fl.png`: Spectra plots + - `*.hapmers.blob.png`: Hap-mer blob plot +
+ +[Merqury](https://github.com/marbl/merqury) is used for the k-mer analysis. + +
+AssemblyQC - Spectra-cn plot +AssemblyQC - Plotsr synteny plot +
+AssemblyQC - Merqury plots +
+ ### Synteny
@@ -225,37 +239,42 @@ Hi-C contact mapping experiments measure the frequency of physical contact betwe - `plotsr.png`: Plotsr synteny plot
-[Circos](https://circos.ca) and linear synteny plots are created from genome-wide alignments performed with [MUMMER](https://github.com/mummer4/mummer?tab=readme-ov-file) and [`dnadiff.pl`](https://github.com/mummer4/mummer/blob/master/scripts/dnadiff.pl). +[Circos](https://circos.ca) and dotplots are created from genome-wide alignments performed with [MUMmer](https://github.com/mummer4/mummer?tab=readme-ov-file). Whereas, [Plotsr](https://github.com/schneebergerlab/plotsr) plots are created from genome-wide alignments performed with [Minimap2](https://github.com/lh3/minimap2).
AssemblyQC - Circos synteny plot AssemblyQC - Plotsr synteny plot -AssemblyQC - Dotplot synteny plot +AssemblyQC - dotplot synteny plot
AssemblyQC - Synteny plots
-### Merqury +### GenomeTools gt stat
Output files -- `merqury/` - - `tag1-and-tag2`: Results folder for haplotype `tag1` and `tag2`. - - `*.completeness.stats`: Assembly completeness statistics - - `*.qv`: Assembly consensus quality QV statistics - - `*.fl.png`: Spectra plots - - `*.hapmers.blob.png`: Hap-mer blob plot -
+- `genometools_gt_stat/` + - `*.gt.stat.yml`: Assembly annotation stats in yaml format. -[MERQURY](https://github.com/marbl/merqury) is used for the k-mer analysis. + -
-AssemblyQC - Spectra-cn plot -AssemblyQC - Plotsr synteny plot -
-AssemblyQC - Merqury plots -
+GenomeTools `gt stat` tool calculates a basic set of statistics about features contained in GFF3 files. + +
AssemblyQC - GenomeTools gt stat gene length distribution
AssemblyQC - GenomeTools gt stat gene length distribution
+ +### OrthoFinder + +
+Output files + +- `orthofinder/assemblyqc`: OrthoFinder output folder. + +
+ +If more than one assemblies are included along with their annotations, OrthoFinder is executed on the annotation proteins to perform a phylogenetic orthology inference for comparative genomics. + +
AssemblyQC - OrthoFinder species tree
AssemblyQC - OrthoFinder species tree
### Pipeline information diff --git a/docs/parameters.md b/docs/parameters.md index 7ef7843a..1c8fd3f6 100644 --- a/docs/parameters.md +++ b/docs/parameters.md @@ -34,6 +34,15 @@ A Nextflow pipeline which evaluates assembly quality with multiple QC tools and | `ncbi_fcs_gx_db_path` | Path to NCBI FCS GX database. See: https://github.com/ncbi/fcs/wiki/FCS-GX | `string` | | | | | `contamination_stops_pipeline` | Skip remaining QC steps for an assembly which has adaptor or GX contamination | `boolean` | True | | | +## tidk options + +| Parameter | Description | Type | Default | Required | Hidden | +| --------------------- | ---------------------------------------------------------------------------------------------------------- | --------- | ------- | -------- | ------ | +| `tidk_skip` | Skip telomere identification | `boolean` | True | | | +| `tidk_repeat_seq` | Telomere repeat sequence. Typical values for plant: TTTAGGG, fungus, vertebrates: TTAGGG and Insect: TTAGG | `string` | | | | +| `tidk_filter_by_size` | Filter assembly sequences smaller than the specified length | `boolean` | | | | +| `tidk_filter_size_bp` | Filter size in base-pairs | `integer` | 1000000 | | | + ## BUSCO options | Parameter | Description | Type | Default | Required | Hidden | @@ -43,22 +52,13 @@ A Nextflow pipeline which evaluates assembly quality with multiple QC tools and | `busco_lineage_datasets` | BUSCO lineages. It should be provided as a space-separated list of lineages: 'fungi_odb10 microsporidia_odb10' | `string` | | | | | `busco_download_path` | Download path for BUSCO | `string` | | | | -## TIDK options - -| Parameter | Description | Type | Default | Required | Hidden | -| --------------------- | ---------------------------------------------------------------------------------------------------------- | --------- | ------- | -------- | ------ | -| `tidk_skip` | Skip telomere identification | `boolean` | True | | | -| `tidk_repeat_seq` | Telomere repeat sequence. Typical values for plant: TTTAGGG, fungus, vertebrates: TTAGGG and Insect: TTAGG | `string` | | | | -| `tidk_filter_by_size` | Filter assembly sequences smaller than the specified length | `boolean` | | | | -| `tidk_filter_size_bp` | Filter size in base-pairs | `integer` | 1000000 | | | - ## LAI options | Parameter | Description | Type | Default | Required | Hidden | | ---------- | ------------------- | --------- | ------- | -------- | ------ | | `lai_skip` | Skip LAI estimation | `boolean` | True | | | -## Kraken2 options +## Kraken 2 options | Parameter | Description | Type | Default | Required | Hidden | | ----------------- | --------------------- | --------- | ------- | -------- | ------ | @@ -75,6 +75,13 @@ A Nextflow pipeline which evaluates assembly quality with multiple QC tools and | `hic_fastp_ext_args` | Additional parameters for fastp trimming | `string` | --qualified_quality_phred 20 --length_required 50 | | | | `hic_samtools_ext_args` | Additional parameters for samtools view command run after samblaster | `string` | -F 3852 | | | +## Merqury options + +| Parameter | Description | Type | Default | Required | Hidden | +| --------------------- | -------------------------------- | --------- | ------- | -------- | ------ | +| `merqury_skip` | Skip merqury analysis | `boolean` | True | | | +| `merqury_kmer_length` | kmer length for merqury analysis | `integer` | 21 | | | + ## Synteny options | Parameter | Description | Type | Default | Required | Hidden | @@ -93,14 +100,7 @@ A Nextflow pipeline which evaluates assembly quality with multiple QC tools and | `synteny_plotsr_seq_label` | Sequence label prefix for plotsr synteny | `string` | Chr | | | | `synteny_plotsr_assembly_order` | The order in which the assemblies should be compared, provided as space separated string of assembly tags. If absent, assemblies are ordered by their tags alphabetically. | `string` | | | | -## Merqury options - -| Parameter | Description | Type | Default | Required | Hidden | -| --------------------- | -------------------------------- | --------- | ------- | -------- | ------ | -| `merqury_skip` | Skip merqury analysis | `boolean` | True | | | -| `merqury_kmer_length` | kmer length for merqury analysis | `integer` | 21 | | | - -## Orthofinder options +## OrthoFinder options | Parameter | Description | Type | Default | Required | Hidden | | ------------------ | ---------------- | --------- | ------- | -------- | ------ | diff --git a/docs/usage.md b/docs/usage.md index 3e8b314b..02b198b3 100644 --- a/docs/usage.md +++ b/docs/usage.md @@ -2,17 +2,17 @@ - [Assemblysheet input](#assemblysheet-input) - [External databases](#external-databases) - - [NCBI FCS GX database](#ncbi-fcs-gx-database) - - [Kraken2](#kraken2) + - [NCBI FCS-GX database](#ncbi-fcs-gx-database) - [BUSCO](#busco) + - [Kraken 2](#kraken-2) - [Other parameters](#other-parameters) - [Assemblathon stats](#assemblathon-stats) - - [NCBI FCS GX](#ncbi-fcs-gx) + - [NCBI FCS-GX](#ncbi-fcs-gx) + - [tidk](#tidk) - [BUSCO](#busco-1) - - [TIDK](#tidk) - [HiC](#hic) - - [Synteny analysis](#synteny-analysis) - [Merqury K-mer analysis](#merqury-k-mer-analysis) + - [Synteny analysis](#synteny-analysis) - [Minimum System Requirements](#minimum-system-requirements) - [Running the pipeline](#running-the-pipeline) - [Updating the pipeline](#updating-the-pipeline) @@ -43,9 +43,9 @@ See the [Merqury](#merqury-k-mer-analysis) section For description of assemblysh ## External databases -### NCBI FCS GX database +### NCBI FCS-GX database -If NCBI FCS GX foreign organism contamination check is executed by setting `ncbi_fcs_gx_skip` to `false`, the path to the GX database must be provided with option `ncbi_fcs_gx_db_path`. The user must ensure that the database is correctly downloaded and placed in a location accessible to the pipeline. Setup instructions are available at . The database path must contain following files: +If NCBI FCS-GX foreign organism contamination check is executed by setting `ncbi_fcs_gx_skip` to `false`, the path to the GX database must be provided with option `ncbi_fcs_gx_db_path`. The user must ensure that the database is correctly downloaded and placed in a location accessible to the pipeline. Setup instructions are available at . The database path must contain following files: ```bash all.assemblies.tsv @@ -59,14 +59,14 @@ all.seq_info.tsv.gz all.taxa.tsv ``` -### Kraken2 - -Path to Kraken2 database is provided by the `kraken2_db_path` parameter. This can be a URL to a public `.tar.gz` file such as `https://genome-idx.s3.amazonaws.com/kraken/k2_pluspfp_20240112.tar.gz`. The pipeline can download and extract the database. This is not the recommended practice owing to the size of the database. Rather, the database should be downloaded, extracted and stored in a read-only location. The path to that location can be provided by the `kraken2_db_path` parameter such as `/workspace/ComparativeDataSources/kraken2db/k2_pluspfp_20230314`. - ### BUSCO BUSCO lineage databases are downloaded and updated by the BUSCO tool itself. A persistent location for the database can be provided by specifying `busco_download_path` parameter. +### Kraken 2 + +Path to Kraken 2 database is provided by the `kraken2_db_path` parameter. This can be a URL to a public `.tar.gz` file such as `https://genome-idx.s3.amazonaws.com/kraken/k2_pluspfp_20240112.tar.gz`. The pipeline can download and extract the database. This is not the recommended practice owing to the size of the database. Rather, the database should be downloaded, extracted and stored in a read-only location. The path to that location can be provided by the `kraken2_db_path` parameter such as `/workspace/ComparativeDataSources/kraken2db/k2_pluspfp_20230314`. + ## Other parameters This section provides additional information for parameters. It does not list all the pipeline parameters. For an exhaustive list, see [parameters.md](./parameters.md). @@ -75,35 +75,22 @@ This section provides additional information for parameters. It does not list al `assemblathon_stats_n_limit` is the number of 'N's for the unknown gap size. This number is used to split the scaffolds into contigs to compute contig-related stats. NCBI's recommendation for unknown gap size is 100 . -### NCBI FCS GX +### NCBI FCS-GX - `ncbi_fcs_gx_tax_id` is the taxonomy ID for all the assemblies listed in the assemblysheet. A taxonomy ID can be obtained by searching a _Genus species_ at . -### BUSCO +### tidk -- `busco_lineage_datasets`: A space-separated list of BUSCO lineages. Any number of lineages can be specified such as "fungi_odb10 hypocreales_odb10". Each assembly is assessed against each of the listed lineage. To select a lineage, refer to . +- `tidk_repeat_seq`: The telomere search sequence. To select an appropriate sequence, see . Commonly used sequences are TTTAGGG (Plant), TTAGGG (Fungus, Vertebrates) and TTAGG (Insect). Further reading: -### TIDK +### BUSCO -- `tidk_repeat_seq`: The telomere search sequence. To select an appropriate sequence, see . Commonly used sequences are TTTAGGG (Plant), TTAGGG (Fungus, Vertebrates) and TTAGG (Insect). Further reading: +- `busco_lineage_datasets`: A space-separated list of BUSCO lineages. Any number of lineages can be specified such as "fungi_odb10 hypocreales_odb10". Each assembly is assessed against each of the listed lineage. To select a lineage, refer to . ### HiC - `hic`: Path to reads provided as a SRA ID or as a path to paired reads such as 'hic_reads{1,2}.fastq.gz'. These reads are applied to each assembly listed by `input`. -### Synteny analysis - -- `synteny_xref_assemblies`: Similar to `--input`, this parameter also provides a CSV sheet listing external reference assemblies which are included in the synteny analysis but are not analysed by other QC tools. See the [example xrefsheet](../assets/xrefsheet.csv) included with the pipeline. Its fields are: - - - `tag:` A unique tag which represents the reference assembly in the final report - - `fasta:` FASTA file - - `synteny_labels:` A two column tsv file listing fasta sequence ids (first column) and their labels for the synteny plots (second column) - -- `synteny_plotsr_assembly_order`: The order in which Minimap2 alignments are performed and, then, plotted by Plotsr. For assembly A, B and C; if the order is specified as 'B C A', then, two alignments are performed. First, C is aligned against B as reference. Second, A is aligned against C as reference. The order of these assemblies on the Plotsr figure is also 'B C A' so that B appears on top, C in the middle and A at the bottom. If this parameter is `null`, the assemblies are ordered alphabetically. All assemblies from `input` and `synteny_xref_assemblies` are included by default. If an assembly is missing from this list, that assembly is excluded from the analysis. - -> [!WARNING] -> PLOTSR performs a sequence-wise (preferably chromosome-wise) synteny analysis. The order of the sequences for each assembly is inferred from its `synteny_labels` file and the order of sequences in the FASTA file is ignored. As all the assemblies are included in a single plot and the number of sequences from each assembly should be same, sequences after the common minimum number are excluded. Afterwards, the sequences are marked sequentially as `Chr1`, `Chr2`, `Chr3`,... If a label other than `Chr` is desirable, it can be configured with the `synteny_plotsr_seq_label` parameter. - ### Merqury K-mer analysis Additional assemblysheet columns: @@ -124,6 +111,19 @@ See following assemblysheet examples for MERQURY analysis. The data for these examples comes from: [umd.edu](https://obj.umiacs.umd.edu/marbl_publications/triobinning/index.html) +### Synteny analysis + +- `synteny_xref_assemblies`: Similar to `--input`, this parameter also provides a CSV sheet listing external reference assemblies which are included in the synteny analysis but are not analysed by other QC tools. See the [example xrefsheet](../assets/xrefsheet.csv) included with the pipeline. Its fields are: + + - `tag:` A unique tag which represents the reference assembly in the final report + - `fasta:` FASTA file + - `synteny_labels:` A two column tsv file listing fasta sequence ids (first column) and their labels for the synteny plots (second column) + +- `synteny_plotsr_assembly_order`: The order in which Minimap2 alignments are performed and, then, plotted by Plotsr. For assembly A, B and C; if the order is specified as 'B C A', then, two alignments are performed. First, C is aligned against B as reference. Second, A is aligned against C as reference. The order of these assemblies on the Plotsr figure is also 'B C A' so that B appears on top, C in the middle and A at the bottom. If this parameter is `null`, the assemblies are ordered alphabetically. All assemblies from `input` and `synteny_xref_assemblies` are included by default. If an assembly is missing from this list, that assembly is excluded from the analysis. + +> [!WARNING] +> PLOTSR performs a sequence-wise (preferably chromosome-wise) synteny analysis. The order of the sequences for each assembly is inferred from its `synteny_labels` file and the order of sequences in the FASTA file is ignored. As all the assemblies are included in a single plot and the number of sequences from each assembly should be same, sequences after the common minimum number are excluded. Afterwards, the sequences are marked sequentially as `Chr1`, `Chr2`, `Chr3`,... If a label other than `Chr` is desirable, it can be configured with the `synteny_plotsr_seq_label` parameter. + ## Minimum System Requirements All the modules have been tested to work on a single machine with 10 CPUs + 32 GBs of memory, except NCBI FCS GX and Kraken2. Their minimum requirements are: diff --git a/nextflow.config b/nextflow.config index b313f4ac..0e7c24bb 100644 --- a/nextflow.config +++ b/nextflow.config @@ -29,22 +29,22 @@ params { contamination_stops_pipeline = true + // tidk options + tidk_skip = true + tidk_repeat_seq = null + tidk_filter_by_size = false + tidk_filter_size_bp = 1000000 + // BUSCO options busco_skip = true busco_mode = null busco_lineage_datasets = null busco_download_path = null - // TIDK options - tidk_skip = true - tidk_repeat_seq = null - tidk_filter_by_size = false - tidk_filter_size_bp = 1000000 - // LAI options lai_skip = true - // kraken2 options + // kraken 2 options kraken2_skip = true kraken2_db_path = null @@ -55,6 +55,10 @@ params { hic_fastp_ext_args = '--qualified_quality_phred 20 --length_required 50' hic_samtools_ext_args = '-F 3852' + // Merqury options + merqury_skip = true + merqury_kmer_length = 21 + // Synteny options synteny_skip = true synteny_mummer_skip = true @@ -70,11 +74,7 @@ params { synteny_plotsr_seq_label = 'Chr' synteny_plotsr_assembly_order = null - // Merqury options - merqury_skip = true - merqury_kmer_length = 21 - - // Orthofinder options + // OrthoFinder options orthofinder_skip = true // Output options diff --git a/nextflow_schema.json b/nextflow_schema.json index 164e8c9b..a457f73f 100644 --- a/nextflow_schema.json +++ b/nextflow_schema.json @@ -113,6 +113,37 @@ } } }, + "tidk_options": { + "title": "tidk options", + "type": "object", + "description": "", + "default": "", + "properties": { + "tidk_skip": { + "type": "boolean", + "description": "Skip telomere identification", + "default": true, + "fa_icon": "fas fa-forward" + }, + "tidk_repeat_seq": { + "type": "string", + "description": "Telomere repeat sequence. Typical values for plant: TTTAGGG, fungus, vertebrates: TTAGGG and Insect: TTAGG", + "pattern": "^[ACGT]+$", + "fa_icon": "fas fa-dna" + }, + "tidk_filter_by_size": { + "type": "boolean", + "description": "Filter assembly sequences smaller than the specified length", + "fa_icon": "fas fa-question-circle" + }, + "tidk_filter_size_bp": { + "type": "integer", + "default": 1000000, + "description": "Filter size in base-pairs", + "fa_icon": "fas fa-ruler-horizontal" + } + } + }, "busco_options": { "title": "BUSCO options", "type": "object", @@ -145,37 +176,6 @@ } } }, - "tidk_options": { - "title": "TIDK options", - "type": "object", - "description": "", - "default": "", - "properties": { - "tidk_skip": { - "type": "boolean", - "description": "Skip telomere identification", - "default": true, - "fa_icon": "fas fa-forward" - }, - "tidk_repeat_seq": { - "type": "string", - "description": "Telomere repeat sequence. Typical values for plant: TTTAGGG, fungus, vertebrates: TTAGGG and Insect: TTAGG", - "pattern": "^[ACGT]+$", - "fa_icon": "fas fa-dna" - }, - "tidk_filter_by_size": { - "type": "boolean", - "description": "Filter assembly sequences smaller than the specified length", - "fa_icon": "fas fa-question-circle" - }, - "tidk_filter_size_bp": { - "type": "integer", - "default": 1000000, - "description": "Filter size in base-pairs", - "fa_icon": "fas fa-ruler-horizontal" - } - } - }, "lai_options": { "title": "LAI options", "type": "object", @@ -191,7 +191,7 @@ } }, "kraken2_options": { - "title": "Kraken2 options", + "title": "Kraken 2 options", "type": "object", "description": "", "default": "", @@ -246,6 +246,27 @@ } } }, + "merqury_options": { + "title": "Merqury options", + "type": "object", + "description": "", + "default": "", + "properties": { + "merqury_skip": { + "type": "boolean", + "default": true, + "description": "Skip merqury analysis", + "fa_icon": "fas fa-forward" + }, + "merqury_kmer_length": { + "type": "integer", + "default": 21, + "description": "kmer length for merqury analysis", + "minimum": 3, + "fa_icon": "fas fa-ruler-horizontal" + } + } + }, "synteny_options": { "title": "Synteny options", "type": "object", @@ -333,29 +354,8 @@ } } }, - "merqury_options": { - "title": "Merqury options", - "type": "object", - "description": "", - "default": "", - "properties": { - "merqury_skip": { - "type": "boolean", - "default": true, - "description": "Skip merqury analysis", - "fa_icon": "fas fa-forward" - }, - "merqury_kmer_length": { - "type": "integer", - "default": 21, - "description": "kmer length for merqury analysis", - "minimum": 3, - "fa_icon": "fas fa-ruler-horizontal" - } - } - }, "orthofinder_options": { - "title": "Orthofinder options", + "title": "OrthoFinder options", "type": "object", "description": "", "default": "", @@ -466,10 +466,10 @@ "$ref": "#/$defs/ncbi_fcs_options" }, { - "$ref": "#/$defs/busco_options" + "$ref": "#/$defs/tidk_options" }, { - "$ref": "#/$defs/tidk_options" + "$ref": "#/$defs/busco_options" }, { "$ref": "#/$defs/lai_options" @@ -481,10 +481,10 @@ "$ref": "#/$defs/hic_options" }, { - "$ref": "#/$defs/synteny_options" + "$ref": "#/$defs/merqury_options" }, { - "$ref": "#/$defs/merqury_options" + "$ref": "#/$defs/synteny_options" }, { "$ref": "#/$defs/orthofinder_options"