last-genome-alignments
last-genome-alignments copied to clipboard
Human self-alignments?
Does last
have the ability to efficiently compute self-alignments? lastz
has a mode specifically for this, but I was hoping to be able to replicate that more rapidly with last
.
Unfortunately last
does not (yet) have such a mode. Maybe one day, but probably not in the near future. Naive self-alignment with last
is not catastrophically inefficient, at least.
I'll also mention that it's unclear to me what genome self-alignment should ideally do. Find all pairs of paralogous sequences, no matter how ancient their common ancestor? Align all pairs of Alus - surely not?
It might help to detect segmental duplications. Lastz does have a self-alignment and dot plot function.