Consider also that if have a read that has 1 match with 0 errors plus, say, 3 additional matches with 1 error than:
"-m 1 --best --strata" will report the match with 0 errors (because the match is unique in the best error stratum)
While:
"-m 1" or "-m 1 --best" will suppress all the alignments and report the read as unmapped.
So "-m 1" without "--best --strata" gives the strongest guarantee that a match is unique, although in my opinion it is too conservative.
Please correct me if I'm wrong...
Dario
|