Skip to content

Futile gap opening #2

Description

@jnvaughn

Hi,

When using progressiveMauve (default parameters), we often observe an issue I've call "futile gap opening". Here is an example:

ACTAGCTGATG--------CTGACGTAATCGTGATGATCGATGCTGA 
ACTAGCTGATGCTGACGTA--------ATCGTGATGATCGATGCTGA
ACTAGCTGATGCTGACGTA--------ATCGTGATGATCGATGCTGA

As you can see, these gaps are clearly not warranted.

The issue is quite common. The histogram below is from whole genome alignments. I won't describe the x-axis in detail but the sharp peak at 0.5 reflect these futile gap openings as a fraction of all gap events in the total alignment.

Screen Shot 2021-03-16 at 5 55 44 PM

Does anyone know the nature of this issue and how to avoid?

Thanks.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions