This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Understanding MUMmer show-diff

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.

mummer

Hello Huggie,

Did you find the answer?

0 answers

No answers yet.

Log in to answer this question.