Skip to content

Mylotools overview

Mylotools is a set of utilities for working with myloasm assemblies. Source code for mylotools can be found here.

The main use of mylotools is to visualize assembled contigs for quality control. There are also functions for processing myloasm's contigs.

A QC plot from `mylotools report` or `mylotools plot`

Mylotools currently has the following useful functionality:

usage: mylotools [-h] [--version] {report,sanitize-headers,plot,strain-viz,extract-contigs} ...

positional arguments:
  {report,sanitize-headers,plot,strain-viz,extract-contigs}
    report              Generate comprehensive report with plots for all long contigs
    sanitize-headers    Sanitize FASTA headers by replacing underscores with spaces (creates backup)
    plot                Generate a plot of various statistics for one of myloasm's output contigs (`report` does `plot` for all contigs)
    strain-viz          Visualize overlaps between and within two or more similar contigs
    extract-contigs     Extract contigs > X bp into a folder as its own fasta

Installing mylotools

Conda install

mamba install -c bioconda mylotools

Clone + pip

git clone https://github.com/bluenote-1577/mylotools
cd mylotools

## assuming you have python3 installed
pip install . 

Using mylotools

It's easiest to use mylotools by cd-ing into the assembly output.

myloasm reads.fq .... -o myloasm_results
cd myloasm_results

Generating an HTML report

mylotools report -o mylo-report
ls mylo-report/contig_summary_report.html

This extracts all contigs of length > 300 kbp OR designated circular (by default) and then runs mylotools plot to generate QC plots. contig_summary_report.html is an interactive plot -- clicking a contig redirects to its QC plot.

A comprehensive HTML report. Clicking on the contig leads to a QC plot.