MashTree_FASTA
creates a phylogenetic tree using Mash distances.
Mash distances are representations of how many kmers two sequences have in common. These distances are generated by transforming all kmers from a sequence into an integer value with hashing and Bloom filters. The hashed kmers are sorted and a “sketch” is created by only using the kmers that appear at the top of the sorted list. These sketches can be compared by counting the number of hashed kmers they have in common. Mashtree uses a neighbor-joining algorithm to cluster these “distances” into phylogenetic trees.
This workflow also features an optional module, summarize_data
, that creates a presence/absence matrix for the analyzed samples from a list of indicated columns (such as AMR genes, etc.) that can be used in Phandango.
MashTree_Fasta
is run on a set of assembly fastas and creates a phylogenetic tree and matrix. These outputs are passed to a task that will rearrange the matrix to match the order of the terminal ends in the phylogenetic tree.
The optional summarize_data
task performs the following only if all of the data_summary_*
and sample_names
optional variables are filled out:
"amrfinderplus_virulence_genes,amrfinderplus_stress_genes"
, etc. that can be found within the origin Terra data table.amrfinder_amr_genes
column for a sample contains these values: "aph(3')-IIIa,tet(O),blaOXA-193"
, the summarize_data
task will check each sample in the set to see if they also have those AMR genes detected.By default, this task appends a Phandango coloring tag to color all items from the same column the same; this can be turned off by setting the optional phandango_coloring
variable to false
.
✉️ [email protected] | X (formerly Twitter) | LinkedIn | 🌐 Website