SHOW-DIFF(1)

๐Ÿ’… SHOW-DIFF(1) โ€” user commands, but make it cute

okay bestie, so show-diff is the little tool that spills the tea on structural differences โ€” it's part of the mummer package, and yes it's iconic for that ๐Ÿ’–


โœจ NAME โœจ

show-diff โ€” show diff information (part of mummer package) ๐Ÿฆ‹


๐ŸŒธ SYNOPSIS ๐ŸŒธ

show-diff [options]<deltafile>

๐Ÿ’œ DESCRIPTION ๐Ÿ’œ

  • -f โ€” output diff information as AMOS features โœจ
  • -h โ€” display help information ๐Ÿ™‹โ€โ™€๏ธ
  • -H โ€” do not show header ๐Ÿšซ
  • -q โ€” show diff information for queries ๐Ÿ”
  • -r โ€” show diff information for references (this is the default, no cap) ๐Ÿ‘‘

It outputs a list of structural differences for each sequence in the reference and the query, sorted by position. So for a reference sequence R and its matching query sequence Q, the differences get sorted into these lil categories:

  • GAP โ€” a gap between two mutually consistent alignments ๐Ÿ•ณ๏ธ
  • DUP โ€” an inserted duplication ๐Ÿ‘ฏ
  • BRK โ€” other inserted sequence ๐Ÿงฉ
  • JMP โ€” a rearrangement ๐Ÿฆ˜
  • INV โ€” a rearrangement with inversion ๐Ÿ™ƒ
  • SEQ โ€” a rearrangement with another sequence ๐ŸŒ€

The first five columns of the output are: seq ID, feature type, feature start, feature end, and feature length. Extra columns get added depending on the feature type, because of course they do ๐Ÿ’…

heads up! Negative feature lengths mean there are overlapping adjacent alignment blocks. Not a glitch, just math ๐Ÿ˜ญ

Here's what the output looks like for each type:

IDR GAP gap-start gap-end gap-length-R gap-length-Q gap-diff

IDR DUP dup-start dup-end dup-length

IDR BRK gap-start gap-end gap-length

IDR JMP gap-start gap-end gap-length

IDR INV gap-start gap-end gap-length

IDR SEQ gap-start gap-end gap-length prev-sequence next-sequence

Positions always reference the sequence with the given ID, fr fr.

The sum of the fifth column (ignoring negative values) is the total amount of inserted sequence ๐Ÿงฎโœจ Summing the fifth column after removing DUP features gives you the total unique inserted sequence ๐Ÿซถ

important bestie!! Note that unaligned sequence is not counted, and could represent additional "unique" sequences ๐Ÿ‘€ See the documentation for tips on how to interpret these alignment break features, iykyk ๐Ÿ“–

and that's the whole vibe, bestie ๐ŸŽ€ this page is from mummer 3.23, dated April 2015 โ€” go off, you got this ๐Ÿ’–โœจ๐Ÿ‘‘