This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove sequences without alignment in mummerplot

Hi,

I aligned a scaffold assembly to a closely related species with Mummer. So, this is a "many vs one" approach.

I am happy with the tool. However, when I produce the dot plot on the filtered alignments (I used delta-filter using the "-1" option to get unique alignments) the plot contains scaffolds without any alignments (cf upper part of the plot below).

I would like to remove them. I can manually remove them from the "out.gp" file but is there a way to do it directly with mummerplot ? I did not find anything in the user manual. Another solution would be to write a script that modify the "out.gp" file but how to find the sequence ids that have no alignments ?

Mummerplot

mummer mummerplot
• 2,155 views • link
written by Corentin Corentin

Did you try --filter in mummerplot?

• 1 views • link
written by Sej Modha Sej Modha

I tried but the sequences still appear in the dot plot, like before as an empty line.

• 1 views • link
written by Corentin Corentin

0 answers

No answers yet.

Log in to answer this question.

More posts like this