seqkit
seqkit copied to clipboard
[Feature suggestion][amplicon] primer-pair searching limit
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.
It does reduce search time when using many primers, but most of the time, we need to search all reads.
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.