Skip to content
This repository has been archived by the owner on Aug 23, 2022. It is now read-only.

Using a fraction of read length for the maximum allowed mismatch #18

Open
mengzhou opened this issue Sep 29, 2015 · 1 comment
Open

Comments

@mengzhou
Copy link
Member

Currently the -m option is an integer indicating the number of mismatch allowed per read. But in practice it is better for this number to be automatically determined based on input read length. In rmap, this option can be a float number, e.g. 0.1, to indicate the fraction. Therefore by default the maximum mismatch per map is read_length * 0.1.

For walt it is necessary to calculate such number per read during mapping, because it is possible that the input reads are processed after adapter trimming, and have different lengths.

@acgtun
Copy link
Contributor

acgtun commented Sep 29, 2015

This is a good idea.

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Projects
None yet
Development

No branches or pull requests

2 participants