Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

How to force NUCmer to use exact match instead of unique match?
0
0
Entering edit mode
5.1 years ago
nchuang ▴ 260

I am trying to snip off the flanks of my query sequence when aligning to a reference sequence using Nucmer. I noticed that the alignment:

ACTAGGGAGGA
-----GGGGGA

gets trimmed at:

ACTAGGGA|GGA
-----GGG|GGA

when I want it to trim like:

ACTAG|GGAGGA
-----|GGGGGA

I have tried --maxmatch and --nosimplifiy but it doesn't seem to change the coordinates that it outputs in the .coords file. It has this similar issue with the 3' flank as well.

nucmer mummer alignment • 1.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 106 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6