I've tried to get every possible hits for the mapped read using bwa samse.
The problem with the multiple hits output, is that the output does not have the mismatch information (MD:Z field in the SAM format). It just gives us chromosome positions and edit scores for all hits.
Is there any way to get the mismatch information for all multiple hits from bwa samse -n?
Thanks for the answer in advance.
The problem with the multiple hits output, is that the output does not have the mismatch information (MD:Z field in the SAM format). It just gives us chromosome positions and edit scores for all hits.
Is there any way to get the mismatch information for all multiple hits from bwa samse -n?
Thanks for the answer in advance.
Comment