It's hard to say for sure what is going on without looking at your barcodes file and associated sequences - but I suspect there may be some barcode mismatching going on if you are losing almost all of your reads when demultiplexing. I would run a quick manual inspection of a few barcodes within your barcode file to make sure they match up in your sequence file, that you're identifying the correct barcodes-column (i.e. that barcode is the correct column), and that the barcode sequence reads should in fact be reverse complemented prior to demultiplexing.
If all of that seems reasonable - if you're comfortable with sharing your barcodes file and sequence data (private message is fine if you'd prefer not sharing publicly on the forum), I can take a closer look and see if anything stands out on my end!
I attempted to manually inspect a few barcodes and got a number of hits on each, though I am not certain that means everything is fit to go. When I run the demux, it is not that I get no hits, but that I get only a couple hundred reads per sample rather than tens of thousands. I have attempted to run in in both rev-com directions and --p-rev-comp-mapping-barcodes
returns the limited read could while --p-no-rev-comp-mapping-barcodes does not return anything. Thank you for any help with this.
The files are too large to upload here, but the sequences and qiime objects I have created are in this google drive folder!
Thanks for providing all of your data - I took a look at your demux details, and did see that high read loss. I tried re-running demux emp-single with both --p-rev-comp-mapping-barcodes and --p-rev-comp-barcodes and the sequence retention was much higher (around 75%).
I double checked that golay error correction wasn't leaving out any sequences as well by re-running emp-single with just --p-no-golay-error-correction but only retained around 60% of the sequences.
It seems like the quality of your barcodes is why we're only seeing around a 75% sequence retention when demultiplexing - but with that being said, you should still be able to work with this vs. losing almost all of your reads.
So, just to recap - you'll see much higher sequence retention when running the following command:
Another thing that you may want to do is reach out to your sequencing center to confirm what the correct orientation of your barcode sequences is, and whether or not they are Golay barcodes. It seems as though they are Golay and should be reverse complimented (based on the results from my previous message) - but it would be best to confirm with them so that you know for sure!