From VCFs to SNP distance matrix 1.0
Visit source
Download RO-Crate
Run on GalaxyGiven a set of VCF files and the reference genome used to do the mapping and SNP calling, create a multifasta file containing the genomes of all samples and calculate the matrix of pairwise SNP distances
Associated Tutorial
This workflows is part of the tutorial Identifying tuberculosis transmission links: from SNPs to transmission clusters, available in the GTN
Thanks to...
Tutorial Author(s): Galo A. Goig, Daniela Brites, Christoph Stritt
Tutorial Contributor(s): Wolfgang Maier, Saskia Hiltemann, Helena Rasche, Galo A. Goig, Björn Grüning, Peter van Heusden, Christoph Stritt, Lucille Delisle
Inputs
| ID | Name | Description | Type |
|---|---|---|---|
| Collection of VCFs to analyze | #main/Collection of VCFs to analyze | n/a |
|
| Reference genome of the MTBC ancestor | #main/Reference genome of the MTBC ancestor | n/a |
|
Steps
| ID | Name | Description |
|---|---|---|
| 2 | Filter TB variants | We will ensure at this step that variants to build the MSA are fixed variants and that we low-confidence filter repetitive regions of the MTB genome toolshed.g2.bx.psu.edu/repos/iuc/tb_variant_filter/tb_variant_filter/0.1.3+galaxy0 |
| 3 | Generate the complete genome of each of the samples | The complete genome of each of the samples is generated by inserting the SNPs defined in the respective VCF in the reference genome that was used for mapping and SNP calling toolshed.g2.bx.psu.edu/repos/iuc/bcftools_consensus/bcftools_consensus/1.9+galaxy2 |
| 4 | Concatenate genomes to build a MSA | All genomes are concatenated in a single multifasta file. Because all o them have the same length, this may be seen as a multiple sequence alignment. toolshed.g2.bx.psu.edu/repos/bgruening/text_processing/tp_cat/0.1.1 |
| 5 | Keep only variable positions | Discard invariant positions from the MSA to simplify the file so only contains positions with at least one SNP in at least one strain. toolshed.g2.bx.psu.edu/repos/iuc/snp_sites/snp_sites/2.5.1+galaxy0 |
| 6 | Calculate SNP distances | From the MSA. Calculate pairwise SNP distances between samples. toolshed.g2.bx.psu.edu/repos/iuc/snp_dists/snp_dists/0.6.3+galaxy0 |
Outputs
| ID | Name | Description | Type |
|---|---|---|---|
| {input_file} | #main/{input_file} | n/a |
|
| _anonymous_output_1 | #main/_anonymous_output_1 | n/a |
|
| _anonymous_output_2 | #main/_anonymous_output_2 | n/a |
|
| _anonymous_output_3 | #main/_anonymous_output_3 | n/a |
|
| _anonymous_output_4 | #main/_anonymous_output_4 | n/a |
|
Version History
master
Creators and Submitter
Views: 46 Downloads: 6 Runs: 0
Created: 2nd Jun 2025 at 10:59
AttributionsNone
Related items
Teams: Galaxy Training Network
Organizations: Galaxy
I'm a bot managed by @hexylena to upload workflows from the Galaxy Training Network to the WorkflowHub. If you have any issues with my behaviour please let her know by filing an issue.
Galaxy is an open, web-based platform for accessible, reproducible, and transparent computational biological research.
- Accessible: Users can easily run tools without writing code or using the CLI; all via a user-friendly web interface.
- Reproducible: Galaxy captures all the metadata from an analysis, making it completely reproducible.
- Transparent: Users share and publish analyses via interactive pages that can enhance analyses with user annotations.
- Scalable: Galaxy can run ...
Teams: Galaxy Training Network
Web page: https://galaxyproject.org/
The Galaxy Training Network (GTN) is a collection of hands-on tutorials that are designed to be interactive and are built around Galaxy.
These tutorials can be used for learning and teaching how to use Galaxy for general data analysis, as well as a wide array of hands-on tutorials covering specific domains such as assembly, RNA-Seq analysis, deep learning, climate analysis, and more!
