I’m running QIIME2 2019.7 on a Conda environment. It is installed on my institution’s server so I am running it remotely via command line.
I have 280 multiplexed samples that I am trying to get demultiplexed. I think I have imported correctly using the following command. They were sequenced on an Illumina miseq, and I think they are EMP type data.
qiime tools import
–type EMPPairedEndSequences
–input-path fastq
–output-path multiplexed-seqs.qza
I am having difficulty with demultiplexing. I don’t know what to use for per-sample barcodes in m-barcodes-file and m-barcodes-column. I have a tsv with the forward and reverse barcodes in seperate columns, but I don’t know if these barcodes are the same as the ones I need for demultiplexing. If they are, I don’t know how what to do with them to get QIIME to properly demultiplex them.
I ran the following command telling QIIME to use the forward barcodes, but I got the following error message.
qiime demux emp-paired
–i-seqs multiplex-seqs.qza
–m-barcodes-file barcodes.txt
–m-barcodes-column forward
–o-per-sample-sequences demux.qza
–o-error-correction-details ErrorCorrectionDetails.qza
Plugin error from demux:
A duplicate barcode was detected. The barcode TAAGGCGA was observed for samples LSGBK6 and SGBK2.
Debug info has been saved to /tmp/qiime2-q2cli-err-_1y0f5wp.log
I also tried to import the data as “Multiplexed Paired End Barcode In Sequence,” and use cutadapt to demultiplex, but when I go this route I end up with only 24 samples in the end, which can’t be correct.
qiime tools import
–type ‘MultiplexedPairedEndBarcodeInSequence’
–input-path fastq
–output-path mux.qza
qiime cutadapt demux-paired
–i-seqs mux.qza
–m-forward-barcodes-file barcodes.txt
–m-forward-barcodes-column forward
–m-reverse-barcodes-file barcodes.txt
–m-reverse-barcodes-column reverse
–o-per-sample-sequences redemux.qza
–o-untrimmed-sequences re-untrimmed.qza
Thanks for any help!!