We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Hi! Thank you very much for your tool! Seqkit always makes me happy :)
But I found a small bug.
OS: ubuntu20.04, amd64 (personal computer) Command:
seqkit grep -s -f barcodes.txt input.fastq -o output.fastq -m 1 -R -30:-1
barcodes.txt file contains such sequences:
barcodes.txt
AGTAGGCT GTAGGCTC CTGTACGA TGTACGAC GCACCAAG CACCAAGC TACGTTTC ACGTTTCC
input.fastq contains such read:
input.fastq
@1 NNNNNNNNNNNNNNNNNNNNNNNNNNNN + ############################
This command successfully be done without any problems:
seqkit grep -s -f barcodes.txt input.fq -o output.fq -m 1 -R 1:30
But this command fails:
seqkit grep -s -f barcodes.txt input.fq -o out.fq -m 1 -R -30:-1 [INFO] 8 patterns loaded from file [ERRO] fail to build FMIndex for sequence: 1
The text was updated successfully, but these errors were encountered:
Fix negative regions longer than sequence length. #479
3aa3dc9
Actually, it's a bug in sequence region parsing.
For -R -30:-1, it returns an empty sequence if the sequence is shorter than 30 (28 for sequence 1).
-R -30:-1
Sorry, something went wrong.
No branches or pull requests
Hi! Thank you very much for your tool! Seqkit always makes me happy :)
But I found a small bug.
Prerequisites
OS: ubuntu20.04, amd64 (personal computer)
Command:
barcodes.txt
file contains such sequences:input.fastq
contains such read:Bug
This command successfully be done without any problems:
But this command fails:
The text was updated successfully, but these errors were encountered: