I received four files from the sequencing center with a recent 16S Miseq run:
R1, R2, I1 (all fastq.gz), and a Qiime-type mapping file (with barcodes in it)
I want to demultiplex these files so that I get one R1 and one R2 fastq file for each sample in the end.
However, I'm not sure how to proceed. Most workflows online seem to want to merge the reads, but I'd like to keep them separate.
The lab tech mentioned that the barcodes in the mapping file correspond to the data in the I1 index file somehow, but I'm not sure how exactly these match up.
If anyone knows of a tool that does this, or can suggest a general pseudocode algorithm that I could code up, I'd appreciate it!