Dotplot basics#

A technique (Gibbs and McIntyre) for comparing sequences. All cogent3 sequence collections classes (SequenceCollection, Alignment and ArrayAlignment) have a dotplot method.

The method returns a drawable, as demonstrated below between unaligned sequences.

from cogent3 import load_unaligned_seqs

seqs = load_unaligned_seqs("data/SCA1-cds.fasta", moltype="dna")
draw = seqs.dotplot()

If sequence names are not provided, two randomly chosen sequences are selected (see below). The plot title reflects the parameter values for defining a match. window is the size of the sequence segments being compared. threshold is the number of exact matches within window required for the two sequence segments to be considered a match. gap is the size of a gap between adjacent matches before merging.

Modifying the matching parameters#

If we set window and threshold to be equal, this is equivalent to an exact match approach.

draw = seqs.dotplot(name1="Human", name2="Mouse", window=8, threshold=8)

Displaying dotplot for the reverse complement#

draw = seqs.dotplot(name1="Human", name2="Mouse", rc=True)


clicking on an entry in the legend turns it off

Setting plot attributes#

I’ll modify the title and figure width.

draw = seqs.dotplot(name1="Human", name2="Mouse", rc=True, title="SCA1", width=400)