I have many questions regarding the output of MUMMer's show-diff script. I've read description on show-diff -h as well as in the README of the top-level directory. However I do not understand the output. I'll start with GAP reported by show-diff.
According to README:
[GAP] A gap between two mutually consistent ordered and oriented alignments. gap-length-R is the length of the alignment gap in the reference, gap-length-Q is the length of the alignment gap in the query, and gap-diff is the difference between the two gap lengths. If gap-diff is positive, sequence has been inserted in the reference. If gap-diff is negative, sequence has been deleted from the reference. If both gap-length-R and gap-length-Q are negative, the indel is tandem duplication copy difference.
I made the following illustration trying to understand it.
I do not understand what it means by "If both gap-length-R and gap-length-Q are negative, the indel is tandem duplication copy difference, and how that is different from the reported DUP. Any help is appreciated, thanks.
Hello Huggie,
Did you find the answer?