I’m hoping to get some help with a recurring error with a q2-demux emp-single script I am trying to run.
I have multiplexed 16S single-end reads from an Illumina HiSeq run. After I successfully run qiime tools import for EMPSingleEndSequences and get an output (single_end_sequences.qza), I run the following script:
qiime demux emp-single
The script runs successfully, but in the output demux.qza file, I end up missing ~200 samples from my metadata file.
I’m not sure what the problem is, and if it is some sort of indexing issue I may be facing.
One thing I wonder is if it has anything to do with my need to include the --p-no-golay-error-correction parameter. Without this parameter, my script results in the following error:
“Golay decoding requires 12nt barcodes. The barcode attempting to be decoded (TTTGCATTTTA) is of length 11nt.”
We use 12nt IDT primers, so I’m not sure why the script using default golay parameters results in this error. I’m wondering if this points to an indexing issue and the reason for my lost samples.
Any feedback would be appreciated!