Question: Nucmer and delta-filter not symmetric for reference and query
gravatar for Fabio Marroni
2.7 years ago by
Fabio Marroni2.6k
Fabio Marroni2.6k wrote:

Dear all, I was surprised by the fact that running the following two sets of commands I obtained different results

Command 1:

nucmer --prefix=Ref_1_vs_Ref_2 Ref_1.fasta Ref_2.fasta
delta-filter -1 >
show-coords -rcl > Ref_1_vs_Ref_2_filtered.coords

Command 2:

nucmer --prefix=Ref_2_vs_Ref_1 Ref_2.fasta Ref_1.fasta
delta-filter -1 >
show-coords -rcl > Ref_2_vs_Ref_1_filtered.coords

The difference is already present in the first output file, with one file showing approximately 10% more hits than the second. At present, I am not too worried about this (I am going to use both sense of analysis), but I was really surprised because, theoretically nucmer should perform an all-vs-all alignment, and thus I was expecting it to be symmetrical with respect to query and reference choice.

However, I am afraid I am missing some important information about nucmer. Would you expect it to give different results when switching query and reference?


Also the delta-filter -1, being the intersection of -q and -r (respectively align each reference against the best hit, and align the query against the best hit), should be symmetrical (but as I already told, the problem already arose before the filter).

nucmer multifasta alignment • 1.1k views
ADD COMMENTlink written 2.7 years ago by Fabio Marroni2.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1668 users visited in the last hour