I am trying to compare different databases and the amount of reads that map from my sample, being a collection of short reads from the environment. After looking in depth in the manual, they say the BWA MEM is "accurate", but how accurate is it actually? How close does a sequence have to be to be considered a match to the reference? Additionally, can anyone explain the difference between the -T option, and the -A option? And what the 'match score' is? Is this related to % identity?
Thanks for all the help in advance!