add ChewBBACA #5899
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,15 @@ | ||
| categories: | ||
| - Variant Analysis | ||
| description: BSR-Based Allele Calling Algorithm | ||
| long_description: chewBBACA is a comprehensive pipeline including a set of functions for the creation and validation of whole genome and core genome MultiLocus Sequence Typing (wg/cgMLST) schemas, providing an allele calling algorithm based on Blast Score Ratio that can be run in multiprocessor settings and a set of functions to visualize and validate allele variation in the loci. chewBBACA performs the schema creation and allele calls on complete or draft genomes. | ||
| homepage_url: https://github.com/B-UMMI/chewBBACA/tree/master | ||
| name: chewbbaca | ||
| owner: iuc | ||
| remote_repository_url: https://github.com/galaxyproject/tools-iuc/tree/master/tools/chewbbaca | ||
| auto_tool_repositories: | ||
| name_template: "{{ tool_id }}" | ||
| description_template: "Wrapper for {{ tool_name }}." | ||
| suite: | ||
| name: "suite_chewbbaca" | ||
| description: "A suite of Galaxy tools designed to work with the chewbbaca-tools collection." | ||
| type: repository_suite_definition |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,173 @@ | ||
| <tool id="AlleleCall" name="ChewBBACA AlleleCall" version="@CHEW_VERSION@+galaxy0" python_template_version="3.5" profile="21.05"> | ||
| <description>Determine the allelic profiles of a set of genomes</description> | ||
| <macros> | ||
| <import>macros.xml</import> | ||
| </macros> | ||
| <expand macro="requirements" /> | ||
| <command detect_errors="exit_code"><![CDATA[ | ||
| mkdir './input' && | ||
| mkdir './schema' && | ||
| #for $file in $input_file | ||
| cp $file './input/${file.element_identifier}' && | ||
| #end for | ||
| unzip $input_schema -d './schema' && | ||
nilchia marked this conversation as resolved.
Show resolved
Hide resolved
|
||
| chewBBACA.py AlleleCall | ||
| #if str($training_file) != 'None' | ||
| --ptf $training_file | ||
| #end if | ||
| --bsr $blast_score_ratio | ||
| --l $minimum_length | ||
| --t $translation_table | ||
| --st $size_threshold | ||
| $no_inferred | ||
| --pm $prodigal_mode | ||
| --mode $mode | ||
| --force-continue | ||
| -i './input' -g './schema/schema_seed/' -o './output' | ||
| ]]></command> | ||
| <inputs> | ||
| <param format="fasta" name="input_file" type="data" multiple="true" label="Genome assemblies in FASTA format"/> | ||
| <param format="zip" name="input_schema" type="data" multiple="true" label="Schema Files in zip format" help="The schema directory contains the loci FASTA files and a folder named 'short' that contains the FASTA files with the loci representative alleles."/> | ||
| <section name="advanced" title="Advanced options"> | ||
| <param argument="--training-file" type="data" format="binary" label="Prodigal training file" optional="true" /> | ||
|
There was a problem hiding this comment. Please add help text from the documentation something like "Default is to get training file from the schema"
nilchia marked this conversation as resolved.
Show resolved
Hide resolved
|
||
| <param argument="--blast-score-ratio" type="float" min="0.0" max="1.0" value="0.6" label="BLAST Score Ratio value" /> | ||
|
There was a problem hiding this comment. optional, please do not set any value |
||
| <param argument="--minimum-length" type="integer" min="0" value="201" label="Minimum sequence length value"/> | ||
| <param argument="--translation-table" type="integer" min="0" value="11" label="Genetic code used to predict genes and to translate coding sequences"/> | ||
|
There was a problem hiding this comment. This is also an optional one. Can you also please double check if it is an integer type? There was a problem hiding this comment. Yes it is an integer. |
||
| <param argument="--size-threshold" type="float" min="0" value="0.2" label="CDS size variation threshold"/> | ||
|
There was a problem hiding this comment. This is an optional param. Please do not set any default values to optional params rather use |
||
| <param argument="--no-inferred" type="boolean" truevalue="--no-inferred" falsevalue="" checked="false" label="Add the sequences of inferred alleles (INF) to the schema" help="Use this parameter if the schema is being accessed by multiple processes/users simultaneously." /> | ||
| <param argument="--prodigal-mode" type="select" label="Prodigal Mode" help="single for finished genomes, reasonable quality draft genomes and big viruses. meta for metagenomes, low quality draft genomes, small viruses, and small plasmids"> | ||
| <option value="single" selected="true"> | ||
| single | ||
| </option> | ||
| <option value="meta"> | ||
| meta | ||
| </option> | ||
| </param> | ||
| <param argument="--mode" type="select" label="Execution mode" > | ||
| <option value="1">Only exact matches at DNA level</option> | ||
| <option value="2">Exact matches at DNA and Protein level </option> | ||
| <option value="3">Exact matches and minimizer-based clustering to find similar alleles based on BSR+0.1 </option> | ||
| <option value="4" selected="true">Exact matches and minimizer-based clustering to find similar alleles based on BSR+0.1 </option> | ||
| </param> | ||
| </section> | ||
bgruening marked this conversation as resolved.
Show resolved
Hide resolved
|
||
| </inputs> | ||
| <outputs> | ||
| <collection name="allelecall_results" type="list" label="${tool.name} on ${on_string}: AlleleCall Results"> | ||
| <discover_datasets pattern="(?P<name>.+)\.tsv$" format="tabular" directory="./output"/> | ||
| </collection> | ||
| <collection name="allelcall_log" type="list" label="${tool.name} on ${on_string}: AlleleCall Logs"> | ||
| <discover_datasets pattern="(?P<name>.+)\.txt$" format="txt" directory="./output"/> | ||
| </collection> | ||
| </outputs> | ||
| <tests> | ||
| <test> | ||
| <param name="input_file" value="GCA_000007265.1_ASM726v1_genomic.fna"/> | ||
| <param name="input_schema" value="GCA_000007265.1_ASM726v1_schema_seed.zip"/> | ||
| <output_collection name="allelecall_results" type="list"> | ||
| <element name="cds_coordinates" file="cds_coordinates.tsv" compare="sim_size"/> | ||
| <element name="loci_summary_stats" file="loci_summary_stats.tsv" compare="sim_size"/> | ||
| <element name="paralogous_loci" ftype="tabular"> | ||
| <assert_contents> | ||
| <has_text_matching expression="Genome.*Loci.*CDS"/> | ||
| </assert_contents> | ||
| </element> | ||
| <element name="results_alleles" ftype="tabular"> | ||
| <assert_contents> | ||
| <has_text_matching expression="1.*1.*NIPHEM.*1.*1"/> | ||
| <has_text_matching expression="GCA_000007265.*1"/> | ||
| </assert_contents> | ||
| </element> | ||
| <element name="results_alleles" file="results_alleles.tsv" compare="sim_size"/> | ||
| <element name="results_statistics" file="results_statistics.tsv" compare="sim_size"/> | ||
| </output_collection> | ||
| <output_collection name="allelcall_log" type="list"> | ||
| <element name="logging_info" ftype="txt"> | ||
| <assert_contents> | ||
| <has_text_matching expression="Used a BSR of: 0.6"/> | ||
| </assert_contents> | ||
| </element> | ||
| </output_collection> | ||
| </test> | ||
| </tests> | ||
| <help><![CDATA[ | ||
|
|
||
| chewBBACA version: 3.3.3 | ||
| Authors: Rafael Mamede, Pedro Cerqueira, Mickael Silva, João Carriço, Mário Ramirez | ||
| Github: https://github.com/B-UMMI/chewBBACA | ||
| Documentation: https://chewbbaca.readthedocs.io/en/latest/index.html | ||
| Contacts: [email protected] | ||
|
|
||
| ========================== | ||
| chewBBACA - AlleleCall | ||
| ========================== | ||
| Performs allele calling to determine the allelic profiles of a set of samples in FASTA format. The | ||
| process identifies new alleles, assigns an integer identifier to those alleles and adds them to the | ||
| schema. | ||
|
|
||
| -i, --input-files [INPUT_FILES] Path to the directory with the genome FASTA files or | ||
| to a file that contains a list of full paths to the | ||
| FASTA files, one per line. (default: None) | ||
|
|
||
| -g, --schema-directory SCHEMA_DIRECTORY Path to the schema directory. The schema directory | ||
| contains the loci FASTA files and a folder named | ||
| "short" that contains the FASTA files with the loci | ||
| representative alleles. (default: None) | ||
|
|
||
| -o, --output-directory OUTPUT_DIRECTORY Output directory where the allele calling results | ||
| will be stored (will create a subdirectory named | ||
| "results_<TIMESTAMP>" if the path passed by the user | ||
| already exists). (default: None) | ||
|
|
||
| --ptf, --training-file PTF_PATH Path to the Prodigal training file. Default is to | ||
| get the training file from the schema's directory | ||
| (default: None) | ||
|
|
||
| --bsr, --blast-score-ratio BLAST_SCORE_RATIO BLAST Score Ratio value. Sequences with alignments | ||
| with a BSR value equal to or greater than this value | ||
| will be considered as sequences from the same gene. | ||
| (default: None) | ||
|
|
||
| --l, --minimum-length MINIMUM_LENGTH Minimum sequence length accepted for a coding | ||
| sequence to be included in the schema. (default: | ||
| None) | ||
|
|
||
| --t, --translation-table TRANSLATION_TABLE Genetic code used to predict genes and to translate | ||
| coding sequences. Must match the genetic code used | ||
| to create the training file. (default: None) | ||
|
|
||
| --st, --size-threshold SIZE_THRESHOLD CDS size variation threshold. At the default value | ||
| of 0.2, alleles with size that deviates +-20 percent | ||
| from the locus length mode will be classified as | ||
| ASM/ALM (default: None) | ||
|
|
||
| --pm, --prodigal-mode {single,meta} Prodigal running mode ("single" for finished | ||
| genomes, reasonable quality draft genomes and big | ||
| viruses. "meta" for metagenomes, low quality draft | ||
| genomes, small viruses, and small plasmids). | ||
| (default: single) | ||
|
|
||
| --no-inferred If provided, the process will not add the sequences | ||
| of inferred alleles (INF) to the schema. Allelic | ||
| profiles will still include the allele identifiers | ||
| attributed to the inferred alleles. Use this | ||
| parameter if the schema is being accessed by | ||
| multiple processes/users simultaneously. (default: | ||
| False) | ||
|
|
||
| --mode {1,2,3,4} Execution mode (1: only exact matches at DNA level; | ||
| 2: exact matches at DNA and Protein level; 3: exact | ||
| matches and minimizer-based clustering to find | ||
| similar alleles based on BSR+0.1; 4: runs the full | ||
| process to find exact matches and similar matches | ||
| based on BSR value, including the determination of | ||
| new representative alleles to add to the schema). | ||
| (default: 4) | ||
|
|
||
|
|
||
| It is strongly advised to perform allele calling with the default schema parameters to ensure more | ||
| consistent results. Module documentation available at | ||
| https://chewbbaca.readthedocs.io/en/latest/user/modules/AlleleCall.html | ||
|
|
||
| ]]></help> | ||
| <expand macro="citations" /> | ||
| </tool> | ||
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,93 @@ | ||
| <tool id="AlleleCallEvaluator" name="chewBBACA AlleleCallEvaluator" version="@CHEW_VERSION@+galaxy0" python_template_version="3.5" profile="21.05"> | ||
| <description>Build an interactive report for allele calling results evaluation</description> | ||
| <macros> | ||
| <import>macros.xml</import> | ||
| </macros> | ||
| <expand macro="requirements" /> | ||
| <command detect_errors="exit_code"><![CDATA[ | ||
| mkdir './input' && | ||
| mkdir './schema' && | ||
| #for $file in $input_file | ||
| cp $file './input/${file.element_identifier}' && | ||
| #end for | ||
| unzip $input_schema -d './schema' && | ||
| chewBBACA.py AlleleCallEvaluator | ||
| #if str($annotations) != 'None' | ||
| -a $annotations | ||
| #end if | ||
| $light | ||
| $no_pa | ||
| $no_dm | ||
| $no_tree | ||
| $cg_alignment | ||
| -i './input' -g './schema/schema_seed/' -o '${allelecall_report.files_path}' && | ||
| cp '${allelecall_report.files_path}'/allelecall_report.html $allelecall_report | ||
nilchia marked this conversation as resolved.
Show resolved
Hide resolved
|
||
| ]]></command> | ||
| <inputs> | ||
| <param name="input_file" type="data_collection" collection_type="list" label="AlleleCall Results" format="tabular"/> | ||
| <param name="input_schema" format="zip" type="data" multiple="true" label="Schema Files in zip format" help="The schema directory contains the loci FASTA files and a folder named 'short' that contains the FASTA files with the loci representative alleles."/> | ||
| <section name="advanced" title="Advanced options"> | ||
| <param argument="--annotations" type="data" format="tabular" label="Annotation file" help="File created by the UniprotFinder module" optional="true" /> | ||
| <param argument="--light" type="boolean" truevalue="--light" falsevalue="" checked="false" label="Do not compute the presence-absence matrix, the distance matrix and the Neighbor-Joining tree." /> | ||
| <param argument="--no-pa" type="boolean" truevalue="--no-pa " falsevalue="" checked="false" label="Do not compute the presence-absence matrix" /> | ||
| <param argument="--no-dm" type="boolean" truevalue="--no-dm " falsevalue="" checked="false" label="Do not compute the distance matrix" /> | ||
| <param argument="--no-tree" type="boolean" truevalue="--no-tree " falsevalue="" checked="false" label="Do not compute the Neighbor-Joining tree" /> | ||
| <param argument="--cg-alignment" type="boolean" truevalue="--cg-alignment " falsevalue="" checked="false" label="Compute the MSA of the core genome loci, even if `--no-tree` is provided" /> | ||
| </section> | ||
| </inputs> | ||
| <outputs> | ||
| <data format="html" name="allelecall_report" label="${tool.name} on ${on_string}: AlleleCall Report"/> | ||
| </outputs> | ||
| <tests> | ||
| <test> | ||
| <param name="input_file"> | ||
| <collection type="list"> | ||
| <element name="cds_coordinates.tsv" value="cds_coordinates.tsv" ftype="tabular"/> | ||
| <element name="loci_summary_stats.tsv" value="loci_summary_stats.tsv" ftype="tabular"/> | ||
| <element name="results_alleles.tsv" value="results_alleles.tsv" ftype="tabular"/> | ||
| <element name="results_statistics.tsv" value="results_statistics.tsv" ftype="tabular"/> | ||
| </collection> | ||
| </param> | ||
| <param name="input_schema" value="GCA_000007265.1_ASM726v1_schema_seed.zip"/> | ||
| <output name="allelecall_report" file="allelecall_report.html" compare="sim_size"/> | ||
| </test> | ||
|
There was a problem hiding this comment. please add a second test, if you have an annotation file |
||
| </tests> | ||
| <help><![CDATA[ | ||
| chewBBACA version: 3.3.3 | ||
| Authors: Rafael Mamede, Pedro Cerqueira, Mickael Silva, João Carriço, Mário Ramirez | ||
| Github: https://github.com/B-UMMI/chewBBACA | ||
| Documentation: https://chewbbaca.readthedocs.io/en/latest/index.html | ||
| Contacts: [email protected] | ||
|
|
||
| =================================== | ||
| chewBBACA - AlleleCallEvaluator | ||
| =================================== | ||
|
|
||
| -i, --input-files INPUT_FILES Path to the directory that contains the allele calling | ||
| results. (default: None) | ||
|
|
||
| -g, --schema-directory SCHEMA_DIRECTORY Path to the schema directory. (default: None) | ||
|
|
||
| -o, --output-directory OUTPUT_DIRECTORY Path to the output directory where the report HTML files | ||
| will be generated. (default: None) | ||
|
|
||
| -a, --annotations ANNOTATIONS Path to the TSV file created by the UniprotFinder module. | ||
| (default: None) | ||
|
|
||
| --light Do not compute the presence-absence matrix, the distance | ||
| matrix and the Neighbor-Joining tree. (default: False) | ||
|
|
||
| --no-pa Do not compute the presence-absence matrix. (default: | ||
| False) | ||
|
|
||
| --no-dm Do not compute the distance matrix. (default: False) | ||
|
|
||
| --no-tree Do not compute the Neighbor-Joining tree. (default: | ||
| False) | ||
|
|
||
| --cg-alignment Compute the MSA of the core genome loci, even if --no- | ||
| tree is provided. (default: False) | ||
|
|
||
| ]]></help> | ||
| <expand macro="citations" /> | ||
| </tool> | ||
Oops, something went wrong.
Oops, something went wrong.
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
please change for tools too