seqkit icon indicating copy to clipboard operation
seqkit copied to clipboard

What types of mismatches are supported in `seqkit locate`?

Open jeremymsimon opened this issue 2 years ago • 4 comments

Hi, could you please clarify what sorts of mismatches are supported in seqkit locate? Is it single base changes only or are insertions and/or deletions allowed as well? fish seems to support indels since the alignment parameters mention gap penalties, is that the preferred tool if I want to include indels?

Thanks!

jeremymsimon avatar Jan 24 '24 13:01 jeremymsimon

Yes, seqkit locate only supports base mismatches.

shenwei356 avatar Jan 24 '24 14:01 shenwei356

Got it, thanks. And is there documentation on how the MapQual is calculated in fish? I think RawScore is based on the length of the query and the alignment parameters (+4, -4, -2, -1), but I don't see a description of how that relates to the MapQual

jeremymsimon avatar Jan 24 '24 14:01 jeremymsimon

@botond-sipos might help. He implemented it. :stuck_out_tongue_winking_eye:

shenwei356 avatar Jan 24 '24 15:01 shenwei356

Thanks @shenwei356 @botond-sipos that would be very helpful!

jeremymsimon avatar Feb 01 '24 14:02 jeremymsimon