After doing this, I wanted to perform the quality filters with trimmomatic instead of using cutadapt, so exported the data and found out that many samples had 0 reads (145/170)
Maybe I am doing something wrong, but have no clue. I would really appreciate any help!
Note: after runnig cutadapt demux-paired, I get this error:
WARNING: Adapter 'ACACACAC' (regular 5') was specified multiple times! Please make sure that this is what you want.
note2: As I had different pools for all my sequences, I had to demultiplex each appart. Here I only attach the first one, which contains only 34 samples.
Thank you for your detailed post. I'm not 100% sure what's causing this issue, but hopefully we can find some clues and find your missing reads.
I noticed that all the samples that had reads, had agactatg in the seq_barcode_R column. This makes me think that something is messed up with barcodes. Maybe they are reversed, reversed complimented, or just copied from the wrong file?
Let me know if you find out more about the barcodes, or any other clues!
Other notes:
Because it's a warning (and not an error), you can continue with your analysis as long as that's what you expected!
I also think there is something wrong with the barcodes, but not sure what can it be, I tried to replace the barcodes in the metadata with the reverse and the complimented, but still doesn’t work
I also checked for the reverse barcodes and they are actually present on the data, not sure how to proceed
Here I placed the data from pool1, maybe we can find something with this. Thank you again!
The reverse barcodes were in the reverse.fastq file.
I found something curious, the reads with the only reverse barcodes that worked, had the barcode in the middle of the sequences, but the others (the ones that didn't work), are at the start.. here goes an example:
I wonder if the reads or the barcodes have been flipped...
For the barcode that did work (ACACATGT, near the end), I see multiple hits to the reverse-complement (ACATGTGT) in the seq_barcode_R. This makes sense, as it's near the end of your read and in the seq_barcode_R column.
For the barcode that didn't work (ACGACGAG, at the start), I see an exact match to that in the seq_barcode_R column. (I don't see any matches to the reverse complement, the reverse, or the complement.)
Something weird happed I fliped the barcodes when demultiplexing again as you suggested and I got not 5 as before, but 8 samples with reads, but still missing a lot here I attacht the .qzv
This makes me think that something is actually wrong with the barcodes... but they just sent me this table in excel to compare:
Yeah, something is still wrong up here. After flipping the barcodes, the 8 samples that demultiplexed all have acacacac in the seq_barcode_F column. (or perhaps in the seq_barcode_R column if that's how you flipped it)
This gives me less confidence that either way is correct
Was the sequencing center able to successfully demultiplex your samples? It's possible we are still missing something, but if they are using a custom protocol, they may have a way to solve this problem. If so, they could send you pre-split results with each sample in a pair of fastq files, and we could take it from there!
I fliped the barcodes when demultiplexing again as you suggested and I got not 5 as before, but 8 samples with reads, but still missing a lot 
here I attacht the .qzv
