- [Cite](#cite)
- [About](#about)
- [Installation](#installation)
- [Usage](#usage)
- [Interactive Mode](#interactive-mode)
- [Static Mode](#static-mode)
- [Standard arguments](#standard-arguments)
- [Interactive Mode Commands](#interactive-mode-commands)
- [Static Mode Commands](#static-mode-commands)
- [Sample run - Interactive Mode](#sample-run---interactive-mode)
- [Sample run - Port Forwarding](#sample-run---port-forwarding)
- [Sample run - Static Plots](#sample-run---static-plots)
- [Sample run - comparing two sequences](#sample-run---comparing-two-sequences)
- [Questions](#questions)
- [Known Issues](#known-issues)
## Cite
Alexander P Sweeten, Michael C Schatz, Adam M Phillippy, ModDotPlot—rapid and interactive visualization of tandem repeats, Bioinformatics, Volume 40, Issue 8, August 2024, btae493, https://doi.org/10.1093/bioinformatics/btae493
If you use ModDotPlot for your research, please cite our software!
---
## About
ModDotPlot is a dot plot visualization tool designed for large sequences and whole genomes. ModDotPlot outputs an identity heatmap similar to [StainedGlass](https://mrvollger.github.io/StainedGlass/) by rapidly approximating the Average Nucleotide Identity between pairwise combinations of genomic intervals. This significantly reduces the computational time required to produce these plots, enough to view multiple layers of resolution in real time!
---
## Installation
```
git clone https://github.com/marbl/ModDotPlot.git
cd ModDotPlot
```
Although optional, it's recommended to setup a virtual environment before using ModDotPlot:
```
python -m venv venv
source venv/bin/activate
```
Once activated, you can install the required dependencies:
```
python -m pip install .
```
Finally, confirm that the installation was installed correctly by running `moddotplot -h`:
```
__ __ _ _____ _ _____ _ _
| \/ | | | | __ \ | | | __ \| | | |
| \ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_
| |\/| |/ _ \ / _` | | | | |/ _ \| __| | ___/| |/ _ \| __|
| | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_
|_| |_|\___/ \__,_| |_____/ \___/ \__| |_| |_|\___/ \__|
usage: moddotplot [-h] {interactive,static} ...
ModDotPlot: Visualization of Complex Repeat Structures
positional arguments:
{interactive,static} Choose mode: interactive or static
interactive Interactive mode commands
static Static mode commands
options:
-h, --help show this help message and exit`
```
---
## Usage
ModDotPlot must be run either in `interactive` mode, or `static` mode:
### Interactive Mode
```
moddotplot interactive <ARGS>
```
This will launch a [Dash application](https://plotly.com/dash/) on your machine's localhost. Open any web browser and go to `http://127.0.0.1:<PORT_NUMBER>` to view the interactive plot (this should happen automatically, but depending on your environment you might need to copy and paste this URL into your web browser). Running `Ctrl+C` on the command line will exit the Dash application. The default port number used by Dash is `8050`, but this can be customized using the `--port` command (see [interactive mode commands](#interactive-mode-commands) for further info).
### Static Mode
```
moddotplot static <ARGS>
```
This skips running Dash and quickly creates plots under the specified output directory using [plotnine](https://plotnine.readthedocs.io/en/v0.12.4/). By default, running ModDotPlot in static mode this will produce the following files:
- A paired-end bed file, containing intervals alongside their corresponding identity estimates.
- A self-identity dotplot for each sequence.
- A histogram of identity values for each sequence.
See [static mode commands](#static-mode-commands) for further info.
---
### Standard arguments
The following arguments are the same in both interactive and static mode:
`-f / --fasta <file>`
Fasta files to input. Multifasta files are accepted. Interactive mode will only support a maximum of two sequences at a time.
`-k / --kmer <int>`
K-mer size to use. This should be large enough to distinguish unique k-mers with enough specificity, but not too large that sensitivity is removed. Default: 21.
`-o / --output-dir <string>`
Name of output directory for bed file & plots. Default is current working directory.
`-id / --identity <int>`
Minimum sequence identity cutoff threshold. Default is 86. While it is possible to go as low as 50% sequence identity, anything below 80% is not recommended.
`--delta <float>`
Each partition takes into account a fraction of its neighboring partitions k-mers. This is to avoid sub-optimal identity scores when partitons don't overlap identically. Default is 0.5, and the accepted range is between 0 and 1. Anything greater than 0.5 is not recommended.
`-m / --modimizer <int>`
Modimizer sketch size. Must be lower than window size `w`. A lower sketch size means less k-mers to compare (and faster runtime), at the expense of lower accuracy. Recommended to be kept > 1000.
`-r / --resolution <int>`
Dotplot resolution. This corresponds to the number of windows each input sequence is partitioned into. Default is 1000. Overrides the `--window` parameter.
`--compare <bool>`
If set when 2 or more sequences are input into ModDotPlot, this will show an a vs. b style plot, in addition to a self-identity plot. Note that interactive mode currently only supports a maximum of two sequences. If more than two sequences are input, only the first two will be shown.
`--compare-only <bool>`
If set when 2 or more sequences are input into ModDotPlot, this will show an a vs. b style plot, without showing self-identity plots.
`--ambiguous <bool>`
By default, k-mers that are homopolymers of ambiguous IUPAC codes (eg. NNNNNNNNNNN’s) are excluded from identity estimation. This results in gaps along the central diagonal for these regions. If desired, these can be kept by setting the `—-ambiguous` flag in both interactive and static mode.
`--no-plot <bool>`
Save matrix to file, but skip rendering of plots. In interactive mode, this must be used alongside the `--save` flag
---
### Interactive Mode Commands
`--port <int>`
Port to display ModDotPlot on. Default is 8050, this can be changed to any accepted port.
`-w / --window <int>`
Minimum window size. By default, interactive mode sets a minimum window size based on the sequence length `n/2000` (eg. a 3Mbp sequence will have a 1500bp window). The maximum window size will always be set to `n/1000` (3000bp under the same example). This means that 2 matrices will be created.
`-q / --quick <bool>`
This will automatically run interactive mode with a minimum window size equal to the maximum window size (`n/1000`). This will result in a quick launch, however the resolution of the plot will not improve upon zooming in.
`-s / --save <bool>`
Save the matrices produced in interactive mode. By default, a folder called `interactive_matrices` will be saved in `--output_dir`, containing each matrix in compressed NumPy format, as well as metadata for each matrix in a pickle. Modifying the files in `interactive_matrices` will cause errors when attempting to load them in the future.
`-l / --load <directory>`
Load previously saved matrices. Used instead of `-f/--fasta`
---
### Static Mode Commands
`-c / --config <.json file>`
Run `moddotplot static` with a config file, rather than (sample syntax). Recommended when creating a really customized plot. Used instead of `-f/--fasta`.
`-b / --bed <.bed file>`
Create a plot from a previously computed pairwise bed file. Skips Average Nucleotide Identity computation. Used instead of `-f/--fasta`.
`-w / --window <int>`
Window size. Unlike interactive mode, only one matrix will be created, so this represents the *only* window size. Default is set to `n/1000` (eg. 3000bp for a 3Mbp sequence).
`--no-bed <bool>`
Skip output of bed file.
`--no-hist <bool>`
Skip output of histogram legend.
`--width <int>`
Adjust width of self dot plots. Default is 9 inches.
`--dpi <bool>`
Image resolution in dots per inch (not to be confused with dotplot resolution). Default is 600.
`--palette <str>`
List of accepted palettes can be found [here](https://jiffyclub.github.io/palettable/colorbrewer/). The syntax is to have the name of the palette, followed by an underscore with the number of colors, eg. `OrRd_8`. Default is `Spectral_11`.
`--palette-orientation <bool>`
Flip sequential order of color palette. Set to `-` by default for divergent palettes.
`--color <list of hexcodes>`
List of custom colors in hexcode format can be entered sequentially, mapped from low to high identity.
`--breakpoints <list of ints>`
Add custom identity threshold breakpoints. Note that the number of breakpoints must be equal to the number of colors + 1.
`-t / --axes-ticks <list of ints>`
Custom tickmarks for x and y axis. Values outside of the axes-limits will not be shown.
`-a / --axes-limits <int>`
Change axis limits for x and y axis. Useful for comparing multiple plots, allowing them to stay in scale.
`--bin-freq <bool>`
By default, histograms are evenly spaced based on the number of colors and the identity threshold. Select this argument to bin based on the frequency of observed identity values.
---
### Sample run - Interactive Mode
```
$ moddotplot interactive -f sequences/Chr1_cen.fa
__ __ _ _____ _ _____ _ _
| \/ | | | | __ \ | | | __ \| | | |
| \ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_
| |\/| |/ _ \ / _` | | | | |/ _ \| __| | ___/| |/ _ \| __|
| | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_
|_| |_|\___/ \__,_| |_____/ \___/ \__| |_| |_|\___/ \__|
Running ModDotPlot in interactive mode
Retrieving k-mers from Chr1:14000000-18000000....
Progress: |████████████████████████████████████████| 100.0% Completed
Chr1:14000000-18000000 k-mers retrieved!
Building self-identity matrices for Chr1:14000000-18000000, using a minimum window size of 2000....
Layer 1 using window length 2000
Progress: |████████████████████████████████████████| 100.0% Completed
Layer 2 using window length 4000
Progress: |████████████████████████████████████████| 100.0% Completed
ModDotPlot interactive mode is successfully running on http://127.0.0.1:8050/
Dash is running on http://127.0.0.1:8050/
```
The plotly plot can be navigated using the zoom (magnifying glass) and pan (hand) icons. The plot can be reset by double-clicking or selecting the home button. The identity threshold can be modified by seelcting the slider. Colors can be readjusted according to the same gradient based on the new identity levels.
### Sample run - Port Forwarding
Running interactive mode on an HPC environment can be accomplished through the use of port forwarding. On your remote server, run ModDotPlot as normal:
```
moddotplot interactive -f INPUT_FASTA_FILE(S) --port HPC_PORT_NUMBER
```
Then on your local machine, set up a port forwarding tunnel:
```
ssh -N -f -L <LOCAL_PORT_NUMBER>:127.0.0.1:<HPC_PORT_NUMBER> HPC@LOGIN.CREDENTIALS
```
You should now be able to view interactive mode using `http://127.0.0.1:<LOCAL_PORT_NUMBER>`. Note that your own HPC environment may have specific instructions and/or restrictions for setting up port forwarding.
VSCode now has automatic port forwarding built into the terminal menu. See [VSCode documentation](https://code.visualstudio.com/docs/editor/port-forwarding) for further details
### Sample run - Static Plots
When running ModDotPlot to produce static plots, it is recommended to use a config file, especially when creating extremely customized plots. The config file is provided in JSON, and accepts the same syntax as the command line arguments shown above.
```
$ cat config/config.json
{
"identity": 90,
"sparsity": 10,
"palette": "OrRd_7",
"breakpoints": [
90,
91,
92,
93,
96,
98,
99,
100
],
"output_dir": "Chr1_cen_plots",
"fasta": [
"../sequences/Chr1_cen.fa"
]
}
```
```
$ moddotplot static -c config/config.json
__ __ _ _____ _ _____ _ _
| \/ | | | | __ \ | | | __ \| | | |
| \ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_
| |\/| |/ _ \ / _` | | | | |/ _ \| __| | ___/| |/ _ \| __|
| | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_
|_| |_|\___/ \__,_| |_____/ \___/ \__| |_| |_|\___/ \__|
Running ModDotPlot in static mode
Retrieving k-mers from Chr1:14M-18M....
Progress: |████████████████████████████████████████| 100.0% Completed
Chr1:14M-18M k-mers retrieved!
Computing self identity matrix for chr1:14M-18M...
Sequence length n: 4000000
Window size w: 4000
Modimizer sketch value: 1000
Plot Resolution r: 1000
Progress: |████████████████████████████████████████| 100.0% Completed
Saved bed file to Chr1_cen_plots/Chr1:14M-18M.bed
Plots created! Saving to Chr1_cen_plots/Chr1:14M-18M...
Chr1_cen_plots/Chr1:14M-18M_TRI.png, Chr1_cen_plots/Chr1:14M-18M_TRI.pdf, Chr1_cen_plots/Chr1:14M-18M_FULL.png, Chr1_cen_plots/Chr1:14M-18M_FULL.png, Chr1_cen_plots/Chr1:14M-18M_HIST.png and Chr1_cen_plots/Chr1:14M-18M_HIST.pdf, saved sucessfully.
```
---
### Sample run - comparing two sequences
ModDotPlot can produce an a vs. b style dotplot for each pairwise combination of input sequences. Use the `--compare` command line argument to include these plots. When running `--compare` in interactive mode, a dropdown menu will appear, allowing the user to switch between self-identity and pairwise plots. Note that a maximum of two sequences are allowed in interactive mode. If you want to skip the creation of self-identity plots, you can use `--compare-only`:
```
moddotplot interactive -f sequences/chr15_segment.fa sequences/chr21_segment.fa --compare-only
```
---
## Questions
For bug reports or general usage questions, please raise a GitHub issue, or email alex ~dot~ sweeten ~at~ nih ~dot~ gov
---
## Known Issues
- Plot width and xlim (limiting the x axis to a different amount) currently do not work. I plan to have those working in v0.9.0.
- Mac users might encounter the following unexpected command line output: `/bin/sh: lscpu: command not found`. This is a known issue with Plotnine, the Python plotting library used by ModDotPlot. This can be safely ignored.
- If you encounter an error with the following traceback: `rv = reductor(4) TypeError: cannot pickle 'generator' object`, ths means that you have a newer version of Plotnine that is incompatible with ModDotPlot. Please uninstall plotnine and reinstall version 0.12.4 `pip install plotnine==0.12.4`.
- In interactive mode, comparing sequences of two sizes will lead to errors in zooming for the larger sequence. I plan to fix this in v0.9.0.
Raw data
{
"_id": null,
"home_page": null,
"name": "ModDotPlot",
"maintainer": null,
"docs_url": null,
"requires_python": ">=3.7",
"maintainer_email": "Alex Sweeten <alex.sweeten@nih.gov>",
"keywords": "dotplot, sketching, modimizer, heatmap",
"author": "Michael Schatz, Adam Phillippy",
"author_email": "Alex Sweeten <alex.sweeten@nih.gov>",
"download_url": "https://files.pythonhosted.org/packages/dc/31/5daeea35cf11ff49732dd24b39d7538abb363031cb3ff843da8eb4f9c47d/moddotplot-0.9.0.tar.gz",
"platform": null,
"description": "\n\n- [Cite](#cite)\n- [About](#about)\n- [Installation](#installation)\n- [Usage](#usage)\n - [Interactive Mode](#interactive-mode)\n - [Static Mode](#static-mode)\n - [Standard arguments](#standard-arguments)\n - [Interactive Mode Commands](#interactive-mode-commands)\n - [Static Mode Commands](#static-mode-commands)\n - [Sample run - Interactive Mode](#sample-run---interactive-mode)\n - [Sample run - Port Forwarding](#sample-run---port-forwarding)\n - [Sample run - Static Plots](#sample-run---static-plots)\n - [Sample run - comparing two sequences](#sample-run---comparing-two-sequences)\n- [Questions](#questions)\n- [Known Issues](#known-issues)\n\n\n## Cite\n\nAlexander P Sweeten, Michael C Schatz, Adam M Phillippy, ModDotPlot\u2014rapid and interactive visualization of tandem repeats, Bioinformatics, Volume 40, Issue 8, August 2024, btae493, https://doi.org/10.1093/bioinformatics/btae493\n\nIf you use ModDotPlot for your research, please cite our software!\n\n---\n\n## About\n\nModDotPlot is a dot plot visualization tool designed for large sequences and whole genomes. ModDotPlot outputs an identity heatmap similar to [StainedGlass](https://mrvollger.github.io/StainedGlass/) by rapidly approximating the Average Nucleotide Identity between pairwise combinations of genomic intervals. This significantly reduces the computational time required to produce these plots, enough to view multiple layers of resolution in real time!\n\n\n\n--- \n\n## Installation\n\n```\ngit clone https://github.com/marbl/ModDotPlot.git\ncd ModDotPlot\n```\n\nAlthough optional, it's recommended to setup a virtual environment before using ModDotPlot:\n\n```\npython -m venv venv\nsource venv/bin/activate\n```\n\nOnce activated, you can install the required dependencies:\n\n```\npython -m pip install .\n```\n\nFinally, confirm that the installation was installed correctly by running `moddotplot -h`:\n``` \n __ __ _ _____ _ _____ _ _ \n | \\/ | | | | __ \\ | | | __ \\| | | | \n | \\ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_ \n | |\\/| |/ _ \\ / _` | | | | |/ _ \\| __| | ___/| |/ _ \\| __|\n | | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_ \n |_| |_|\\___/ \\__,_| |_____/ \\___/ \\__| |_| |_|\\___/ \\__|\n\nusage: moddotplot [-h] {interactive,static} ...\n\nModDotPlot: Visualization of Complex Repeat Structures\n\npositional arguments:\n {interactive,static} Choose mode: interactive or static\n interactive Interactive mode commands\n static Static mode commands\n\noptions:\n -h, --help show this help message and exit`\n ```\n--- \n\n## Usage\n\nModDotPlot must be run either in `interactive` mode, or `static` mode:\n\n### Interactive Mode\n\n```\nmoddotplot interactive <ARGS>\n```\n\nThis will launch a [Dash application](https://plotly.com/dash/) on your machine's localhost. Open any web browser and go to `http://127.0.0.1:<PORT_NUMBER>` to view the interactive plot (this should happen automatically, but depending on your environment you might need to copy and paste this URL into your web browser). Running `Ctrl+C` on the command line will exit the Dash application. The default port number used by Dash is `8050`, but this can be customized using the `--port` command (see [interactive mode commands](#interactive-mode-commands) for further info).\n\n### Static Mode\n\n```\nmoddotplot static <ARGS>\n```\n\nThis skips running Dash and quickly creates plots under the specified output directory using [plotnine](https://plotnine.readthedocs.io/en/v0.12.4/). By default, running ModDotPlot in static mode this will produce the following files:\n\n- A paired-end bed file, containing intervals alongside their corresponding identity estimates.\n- A self-identity dotplot for each sequence.\n- A histogram of identity values for each sequence.\n \nSee [static mode commands](#static-mode-commands) for further info.\n\n--- \n\n### Standard arguments\n\nThe following arguments are the same in both interactive and static mode:\n\n`-f / --fasta <file>`\n\nFasta files to input. Multifasta files are accepted. Interactive mode will only support a maximum of two sequences at a time.\n\n`-k / --kmer <int>`\n\nK-mer size to use. This should be large enough to distinguish unique k-mers with enough specificity, but not too large that sensitivity is removed. Default: 21.\n\n`-o / --output-dir <string>`\n\nName of output directory for bed file & plots. Default is current working directory.\n\n`-id / --identity <int>`\n\nMinimum sequence identity cutoff threshold. Default is 86. While it is possible to go as low as 50% sequence identity, anything below 80% is not recommended. \n\n`--delta <float>`\n\nEach partition takes into account a fraction of its neighboring partitions k-mers. This is to avoid sub-optimal identity scores when partitons don't overlap identically. Default is 0.5, and the accepted range is between 0 and 1. Anything greater than 0.5 is not recommended.\n\n`-m / --modimizer <int>`\n\nModimizer sketch size. Must be lower than window size `w`. A lower sketch size means less k-mers to compare (and faster runtime), at the expense of lower accuracy. Recommended to be kept > 1000.\n\n`-r / --resolution <int>`\n\nDotplot resolution. This corresponds to the number of windows each input sequence is partitioned into. Default is 1000. Overrides the `--window` parameter.\n\n`--compare <bool>`\n\nIf set when 2 or more sequences are input into ModDotPlot, this will show an a vs. b style plot, in addition to a self-identity plot. Note that interactive mode currently only supports a maximum of two sequences. If more than two sequences are input, only the first two will be shown.\n\n`--compare-only <bool>`\n\nIf set when 2 or more sequences are input into ModDotPlot, this will show an a vs. b style plot, without showing self-identity plots.\n\n`--ambiguous <bool>`\n\nBy default, k-mers that are homopolymers of ambiguous IUPAC codes (eg. NNNNNNNNNNN\u2019s) are excluded from identity estimation. This results in gaps along the central diagonal for these regions. If desired, these can be kept by setting the `\u2014-ambiguous` flag in both interactive and static mode. \n\n`--no-plot <bool>`\n\nSave matrix to file, but skip rendering of plots. In interactive mode, this must be used alongside the `--save` flag\n\n--- \n\n### Interactive Mode Commands\n\n`--port <int>`\n\nPort to display ModDotPlot on. Default is 8050, this can be changed to any accepted port. \n\n`-w / --window <int>`\n\nMinimum window size. By default, interactive mode sets a minimum window size based on the sequence length `n/2000` (eg. a 3Mbp sequence will have a 1500bp window). The maximum window size will always be set to `n/1000` (3000bp under the same example). This means that 2 matrices will be created.\n\n`-q / --quick <bool>`\n\nThis will automatically run interactive mode with a minimum window size equal to the maximum window size (`n/1000`). This will result in a quick launch, however the resolution of the plot will not improve upon zooming in.\n\n`-s / --save <bool>`\n\nSave the matrices produced in interactive mode. By default, a folder called `interactive_matrices` will be saved in `--output_dir`, containing each matrix in compressed NumPy format, as well as metadata for each matrix in a pickle. Modifying the files in `interactive_matrices` will cause errors when attempting to load them in the future.\n\n`-l / --load <directory>`\n\nLoad previously saved matrices. Used instead of `-f/--fasta`\n\n--- \n\n### Static Mode Commands\n\n`-c / --config <.json file>`\n\nRun `moddotplot static` with a config file, rather than (sample syntax). Recommended when creating a really customized plot. Used instead of `-f/--fasta`.\n\n`-b / --bed <.bed file>`\n\nCreate a plot from a previously computed pairwise bed file. Skips Average Nucleotide Identity computation. Used instead of `-f/--fasta`. \n\n`-w / --window <int>`\n\nWindow size. Unlike interactive mode, only one matrix will be created, so this represents the *only* window size. Default is set to `n/1000` (eg. 3000bp for a 3Mbp sequence). \n\n`--no-bed <bool>`\n\nSkip output of bed file.\n\n`--no-hist <bool>`\n\nSkip output of histogram legend.\n\n`--width <int>`\n\nAdjust width of self dot plots. Default is 9 inches.\n\n`--dpi <bool>`\n\nImage resolution in dots per inch (not to be confused with dotplot resolution). Default is 600.\n\n`--palette <str>`\n\nList of accepted palettes can be found [here](https://jiffyclub.github.io/palettable/colorbrewer/). The syntax is to have the name of the palette, followed by an underscore with the number of colors, eg. `OrRd_8`. Default is `Spectral_11`.\n\n`--palette-orientation <bool>`\n\nFlip sequential order of color palette. Set to `-` by default for divergent palettes. \n\n`--color <list of hexcodes>`\n\nList of custom colors in hexcode format can be entered sequentially, mapped from low to high identity. \n\n`--breakpoints <list of ints>`\n\nAdd custom identity threshold breakpoints. Note that the number of breakpoints must be equal to the number of colors + 1.\n\n`-t / --axes-ticks <list of ints>`\n\nCustom tickmarks for x and y axis. Values outside of the axes-limits will not be shown. \n\n`-a / --axes-limits <int>`\n\nChange axis limits for x and y axis. Useful for comparing multiple plots, allowing them to stay in scale. \n\n`--bin-freq <bool>`\n\nBy default, histograms are evenly spaced based on the number of colors and the identity threshold. Select this argument to bin based on the frequency of observed identity values.\n\n--- \n\n### Sample run - Interactive Mode\n\n```\n$ moddotplot interactive -f sequences/Chr1_cen.fa \n\n __ __ _ _____ _ _____ _ _ \n | \\/ | | | | __ \\ | | | __ \\| | | | \n | \\ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_ \n | |\\/| |/ _ \\ / _` | | | | |/ _ \\| __| | ___/| |/ _ \\| __|\n | | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_ \n |_| |_|\\___/ \\__,_| |_____/ \\___/ \\__| |_| |_|\\___/ \\__|\n\nRunning ModDotPlot in interactive mode\n\nRetrieving k-mers from Chr1:14000000-18000000.... \n\nProgress: |\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588| 100.0% Completed\n\nChr1:14000000-18000000 k-mers retrieved! \n\nBuilding self-identity matrices for Chr1:14000000-18000000, using a minimum window size of 2000.... \n\nLayer 1 using window length 2000\n\nProgress: |\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588| 100.0% Completed\n\n\nLayer 2 using window length 4000\n\nProgress: |\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588| 100.0% Completed\n\n\nModDotPlot interactive mode is successfully running on http://127.0.0.1:8050/ \n\nDash is running on http://127.0.0.1:8050/\n```\n\n\n\nThe plotly plot can be navigated using the zoom (magnifying glass) and pan (hand) icons. The plot can be reset by double-clicking or selecting the home button. The identity threshold can be modified by seelcting the slider. Colors can be readjusted according to the same gradient based on the new identity levels. \n\n### Sample run - Port Forwarding\n\nRunning interactive mode on an HPC environment can be accomplished through the use of port forwarding. On your remote server, run ModDotPlot as normal:\n\n```\nmoddotplot interactive -f INPUT_FASTA_FILE(S) --port HPC_PORT_NUMBER\n```\n\nThen on your local machine, set up a port forwarding tunnel:\n\n```\nssh -N -f -L <LOCAL_PORT_NUMBER>:127.0.0.1:<HPC_PORT_NUMBER> HPC@LOGIN.CREDENTIALS\n```\n\nYou should now be able to view interactive mode using `http://127.0.0.1:<LOCAL_PORT_NUMBER>`. Note that your own HPC environment may have specific instructions and/or restrictions for setting up port forwarding.\n\nVSCode now has automatic port forwarding built into the terminal menu. See [VSCode documentation](https://code.visualstudio.com/docs/editor/port-forwarding) for further details \n\n\n\n### Sample run - Static Plots\n\nWhen running ModDotPlot to produce static plots, it is recommended to use a config file, especially when creating extremely customized plots. The config file is provided in JSON, and accepts the same syntax as the command line arguments shown above. \n\n```\n$ cat config/config.json\n\n{\n \"identity\": 90,\n \"sparsity\": 10,\n \"palette\": \"OrRd_7\",\n \"breakpoints\": [\n 90,\n 91,\n 92,\n 93,\n 96,\n 98,\n 99,\n 100\n ],\n \"output_dir\": \"Chr1_cen_plots\",\n \"fasta\": [\n \"../sequences/Chr1_cen.fa\"\n ]\n}\n```\n\n```\n$ moddotplot static -c config/config.json \n __ __ _ _____ _ _____ _ _ \n | \\/ | | | | __ \\ | | | __ \\| | | | \n | \\ / | ___ __| | | | | | ___ | |_ | |__) | | ___ | |_ \n | |\\/| |/ _ \\ / _` | | | | |/ _ \\| __| | ___/| |/ _ \\| __|\n | | | | (_) | (_| | | |__| | (_) | |_ | | | | (_) | |_ \n |_| |_|\\___/ \\__,_| |_____/ \\___/ \\__| |_| |_|\\___/ \\__|\n\nRunning ModDotPlot in static mode\n\nRetrieving k-mers from Chr1:14M-18M.... \n\nProgress: |\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588| 100.0% Completed\n\nChr1:14M-18M k-mers retrieved! \n\nComputing self identity matrix for chr1:14M-18M... \n\n Sequence length n: 4000000\n\n Window size w: 4000\n\n Modimizer sketch value: 1000\n\n Plot Resolution r: 1000\n\nProgress: |\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588\u2588| 100.0% Completed\n\n\nSaved bed file to Chr1_cen_plots/Chr1:14M-18M.bed\n\nPlots created! Saving to Chr1_cen_plots/Chr1:14M-18M...\n\nChr1_cen_plots/Chr1:14M-18M_TRI.png, Chr1_cen_plots/Chr1:14M-18M_TRI.pdf, Chr1_cen_plots/Chr1:14M-18M_FULL.png, Chr1_cen_plots/Chr1:14M-18M_FULL.png, Chr1_cen_plots/Chr1:14M-18M_HIST.png and Chr1_cen_plots/Chr1:14M-18M_HIST.pdf, saved sucessfully. \n```\n\n\n--- \n\n\n### Sample run - comparing two sequences\n\nModDotPlot can produce an a vs. b style dotplot for each pairwise combination of input sequences. Use the `--compare` command line argument to include these plots. When running `--compare` in interactive mode, a dropdown menu will appear, allowing the user to switch between self-identity and pairwise plots. Note that a maximum of two sequences are allowed in interactive mode. If you want to skip the creation of self-identity plots, you can use `--compare-only`:\n\n```\nmoddotplot interactive -f sequences/chr15_segment.fa sequences/chr21_segment.fa --compare-only\n```\n\n--- \n\n## Questions\n\nFor bug reports or general usage questions, please raise a GitHub issue, or email alex ~dot~ sweeten ~at~ nih ~dot~ gov\n\n--- \n\n## Known Issues\n\n- Plot width and xlim (limiting the x axis to a different amount) currently do not work. I plan to have those working in v0.9.0.\n\n- Mac users might encounter the following unexpected command line output: `/bin/sh: lscpu: command not found`. This is a known issue with Plotnine, the Python plotting library used by ModDotPlot. This can be safely ignored.\n\n- If you encounter an error with the following traceback: `rv = reductor(4) TypeError: cannot pickle 'generator' object`, ths means that you have a newer version of Plotnine that is incompatible with ModDotPlot. Please uninstall plotnine and reinstall version 0.12.4 `pip install plotnine==0.12.4`. \n\n- In interactive mode, comparing sequences of two sizes will lead to errors in zooming for the larger sequence. I plan to fix this in v0.9.0.\n",
"bugtrack_url": null,
"license": null,
"summary": null,
"version": "0.9.0",
"project_urls": null,
"split_keywords": [
"dotplot",
" sketching",
" modimizer",
" heatmap"
],
"urls": [
{
"comment_text": null,
"digests": {
"blake2b_256": "482667f17e93575f10644318756dc2a93958bdf4752250f5753edf2996425930",
"md5": "aeccbb3d7a9c3dcc6ab067c94c08eb01",
"sha256": "1a411d7fb17efa47ce8ad487f6b22cc515577738a31dda09ae6a85570ba47699"
},
"downloads": -1,
"filename": "ModDotPlot-0.9.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "aeccbb3d7a9c3dcc6ab067c94c08eb01",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.7",
"size": 40082,
"upload_time": "2025-01-21T19:35:32",
"upload_time_iso_8601": "2025-01-21T19:35:32.037993Z",
"url": "https://files.pythonhosted.org/packages/48/26/67f17e93575f10644318756dc2a93958bdf4752250f5753edf2996425930/ModDotPlot-0.9.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": null,
"digests": {
"blake2b_256": "dc315daeea35cf11ff49732dd24b39d7538abb363031cb3ff843da8eb4f9c47d",
"md5": "797139e82e4e6a3fef42c3bb275fae6a",
"sha256": "8f41873df01d9eb429de862f2eced4f1860518c938e8658b301fcb767fd1a9c1"
},
"downloads": -1,
"filename": "moddotplot-0.9.0.tar.gz",
"has_sig": false,
"md5_digest": "797139e82e4e6a3fef42c3bb275fae6a",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.7",
"size": 38385,
"upload_time": "2025-01-21T19:35:34",
"upload_time_iso_8601": "2025-01-21T19:35:34.410236Z",
"url": "https://files.pythonhosted.org/packages/dc/31/5daeea35cf11ff49732dd24b39d7538abb363031cb3ff843da8eb4f9c47d/moddotplot-0.9.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2025-01-21 19:35:34",
"github": false,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"lcname": "moddotplot"
}
JSON