seqkit icon indicating copy to clipboard operation
seqkit copied to clipboard

[Feature suggestion][amplicon] primer-pair searching limit

Open Subject-13 opened this issue 2 years ago • 2 comments

Hello!

Might have an idea for your tool:

Prerequisites

  • [x] make sure you're are using the latest version by seqkit version
  • [x] read the usage

Describe your issue

  • [x] describe the problem

It might be useful to collect only a certain number of hits for a primer pair from the given table and then remove the primer pair from the list so that it is no longer searched for in the reads.

  • [x] provide a reproducible example

If I already have 1.000 hits for a particular primer pair, it is no longer necessary to search for the pair. Search in the reads until the limit for each pair is reached or all reads have been searched.

Subject-13 avatar Nov 18 '22 11:11 Subject-13

It does reduce search time when using many primers, but most of the time, we need to search all reads.

shenwei356 avatar Nov 18 '22 12:11 shenwei356

We could use this for our workflow since the nanopore sequencer provides a lot of data, but a certain amount of it becomes redundant for us and just costs time.

I would try to do this adaptation by myself and could give you feedback if it's stable.

Subject-13 avatar Nov 18 '22 13:11 Subject-13