Create K-mer Tree
The Create K-mer Tree tool may be helpful for identification of the closest common reference across samples. The tool uses reads, single sequences or sequence list as input and creates a distance-based phylogenetic tree. There are two ways to initiate creation of a k-mer tree: either from the Result Metadata Table (see the section on Running analysis directly from the Result Metadata Table), or from the Toolbox.
To run the Create K-mer Tree from the toolbox:
Toolbox | Microbial Genomics Module (
)| Typing and Epidemiology (beta) (
) | Trees (
) | Create K-mer Tree (
)
Input files can be specified step-by-step like shown in figure 14.11 or by selecting data recursively by right-clicking on the folder name and selecting Add folder contents (recursively). If using the recursive option, remember to double check that files relevant for the downstream analysis are selected. Click Next.
Figure 14.11: Selection of individual reads and single sequences or sequence list to be included in the K-mer tree analysis.
Now specify the following parameters (figure 14.12).
- K-mer parameters
- K-mer length is the fixed number (k) of DNA bases to search across.
- Only index k-mers with prefix allows specification of the initial bases of the k-mer sequence to limit the search space. Reduction of prefix size increases the RAM requirements, and therefore decrease the search speed.
- Method may be specified by either of the two statistical methods: Jaccard Distance or Feature Frequency Profile via Jensen-Shannon divergences (FFP). You can read more about the Jaccard Distance and FFP at https://en.wikipedia.org/wiki/Jaccard_index and https://en.wikipedia.org/wiki/Alignment-free_sequence_analysis, respectively.
- Strand may be specified as either only the Plus strand or Both strands.
- Result metadata. Specify location of the Result metadata table file.
- Tree view. Select a standard tree setting (i.e., None, K-mer Tree Default or SNP Tree Default) or your own custom tree setting. Read more on creating your customized tree settings: http://www.clcsupport.com/clcgenomicsworkbench/current/index.php?manual=Tree_Settings.html.
Figure 14.12: Various parameters may be set before generation of a K-mer tree.
Subsections