Skip to content
New issue

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

Mismatches allowed set to 2 does not handle single barcode mismatches #2

Open
anniqueclaringbould opened this issue Apr 18, 2024 · 0 comments

Comments

@anniqueclaringbould
Copy link

anniqueclaringbould commented Apr 18, 2024

Hi developers!

Thanks for making defastq, it's super fast and easy to use. However, as I was comparing to a previous result with a different tool, I realised there is a bug when the option --allowed-mismatches set to 2.

I'm adding a minimal working example here of a set of fastq files that have 1 read in them and an index file with just one barcode in it. The reads contains the barcode with 1 mismatch. When I run defastq with --allowed-mismatches 1, it works as expected and the read is demultiplexed. However, when I run with --allowed-mismatches 2, the read ens up in undecoded.R2.fastq.gz.

Code:

/g/scb/zaugg/claringb/scz_CRISPR_screen/tools/defastq \
 -1 test_reads.R1.fastq.gz \
 -2 test_reads.R2.fastq.gz \
 -b read2 \
 -s 27 \
 -l 11 \
 -i test_barcodes.csv \
 -o mismatch2 \
 -a 2 \
 -n 8

/g/scb/zaugg/claringb/scz_CRISPR_screen/tools/defastq \
 -1 test_reads.R1.fastq.gz \
 -2 test_reads.R2.fastq.gz \
 -b read2 \
 -s 27 \
 -l 11 \
 -i test_barcodes.csv \
 -o mismatch1 \
 -a 1 \
 -n 8

Could you look into what's going on here?
Thanks, Annique

test_barcodes.csv
test_reads.R1.fastq.gz
test_reads.R2.fastq.gz

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant