chromo plot

Use chromo plot when you already have a MUMmer coords or minimap2 PAF file and want a visual check without running mummerplot or re-aligning only for the plot. It writes PDF by default and can also write SVG or PNG. Forward-strand alignments are blue and reverse-strand alignments are red.

If you are learning how to read those patterns, start with How to Interpret Dot Plots. It explains axis orientation, forward and reverse segments, chimeric-looking contigs, inversions, duplicate or repeat-like signal, gaps, and common traps.

chromo plot never runs MUMmer or minimap2. It draws the alignment rows you give it. If those rows were generated from raw.fa, the plot still describes raw.fa, even when you also provide a chromo sort assignment report. The assignment report can reorder the query axis by kept sorted contigs, but it does not make a new alignment of <prefix>.ordered.fa.

After any FASTA-changing step, such as chromo sort outputting ordered.fa, chromo clean outputting clean.fa, chromo fix outputting fixed.fa, chromo cut, chromo manual export, or chromo scaffold, re-run MUMmer or minimap2 before using chromo plot to validate that changed FASTA.

Run chromo plot

Whole-genome plot from MUMmer coords:

chromo plot \
  --ref-fasta reference.fa \
  --assembly-fasta assembly.fa \
  --coords mummer/sample.coords \
  --output-prefix plots/sample

This writes plots/sample.pdf.

Whole-genome and per-reference plots from PAF, ordered by a chromo sort assignment report:

chromo plot \
  --ref-fasta reference.fa \
  --assembly-fasta assembly.fa \
  --paf paf/sample.paf \
  --assignments results/sample.contig_assignments.tsv \
  --output-prefix plots/sample \
  --per-ref \
  --sel-ref Gm6 Gm12 Gm15 \
  --formats pdf svg

When --assignments is provided, the query axis is ordered by the kept contigs in the assignment report. This is useful for reviewing a sorted ChromoSort order from the same alignment file that powered the sort. It is not equivalent to plotting the ordered.fa output; align ordered.fa separately if you want the plot to reflect the edited FASTA records.

Per-reference plots keep the x-axis on the full selected reference sequence, but compact each query row to the merged query coordinate spans covered by alignments to that selected reference. This keeps multi-reference contigs and large unaligned query gaps from adding unrelated vertical whitespace to each --per-ref panel. Axis tick units scale with the plotted length, using bp, kb, Mb, or Gb as appropriate for the current panel.

Use --sel-ref when you only want to replot one or a few reference sequences. For example, --sel-ref Gm6 Gm12 Gm15 limits the main plot to those reference sequences and, with --per-ref, writes only those per-reference panels.

chromo plot Outputs

Output Description
<prefix>.pdf Whole-genome PDF dot plot by default.
<prefix>.svg Whole-genome SVG dot plot when --formats svg is set.
<prefix>.png Whole-genome PNG dot plot when --formats png is set.
<prefix>.<ref>.<format> Per-reference plots when --per-ref is set; restricted to selected references when --sel-ref is also set.

Example chromo plot Output

Table 1. Example plot file set. A single run can write a whole-genome plot and, with --per-ref, one plot per reference sequence in each requested format.

Output path Meaning
plots/sample.pdf Whole-genome PDF dot plot.
plots/sample.svg Whole-genome SVG dot plot from the same alignment rows.
plots/sample.chr1.svg Per-reference SVG for chr1.
plots/sample.chr2.svg Per-reference SVG for chr2.
ChromoSort whole-genome dot plot with forward and reverse alignments.
Figure 1. Example whole-genome dot-plot output. PNG output from the synthetic coords fixture showing all reference sequences on the x-axis and assembly contigs on the y-axis; forward alignments are blue and reverse alignments are red.
ChromoSort per-reference dot plot for chr1.
Figure 2. Example per-reference dot-plot output. Per-reference output for chr1 generated with --per-ref, useful for inspecting one chromosome-level slice without re-running an aligner.

chromo plot Parameters

Parameter Default Meaning
--coords required unless --paf MUMmer show-coords alignment file.
--paf required unless --coords minimap2 PAF alignment file.
--formats pdf One or more output formats: pdf, svg, png.
--assignments none Optional chromo sort assignment report for ordering the query axis by kept sorted contigs.
--per-ref off Also write one plot per reference sequence with plotted alignments.
--sel-ref none Limit the main plot and --per-ref output to one or more reference IDs, such as --sel-ref Gm6 Gm12 Gm15.
--per-ref-query-order fasta Use FASTA order or first reference-hit order for per-reference query axes.
--min-segment-bp 0 Minimum query-aligned bp for a row to be drawn.
--min-segment-idy 0.0 Ignore individual alignment rows below this percent identity.
--min-mapq 0 Ignore PAF rows below this MAPQ. Ignored for coords.
--include-secondary-paf off Include PAF rows marked tp:A:S; skipped by default.
--max-segments 0 Maximum drawn alignment rows after filtering; 0 means no limit.