For studying Qiime2, I have forward.fastq.gz, reverse.fastq.gz, barcodes.fastq.gz and mapping_file.tsv downloaded from Qiita (ID 11635). I think it matches the type EMP-paired and it looks like that there are no barocdes and primers in forward.fastq.gz and reverse.fastq.gz.
Here are my operation steps and error reporting
Step1: I have imported my files as EMP-paired using the command:
qiime tools import
--type EMPPairedEndSequences
--input-path /result/emp-paired-end-sequences
--output-path /result
I succeed.
Step2: I used the following code to demultiplex and got the "demux03.qza" successfully.
qiime demux emp-paired
--i-seqs emp-paired-end-sequences.qza
--m-barcodes-file mapping_file.txt
--m-barcodes-column barcode
--p-rev-comp-mapping-barcodes
--o-per-sample-sequences demux03.qza
--o-error-correction-details demux-details03.qza
--verbose
qiime demux summarize
--i-data demux03.qza
--o-visualization demux03.qzv
However, when viewing "demux-full.qzv", the sequence counts dropped dramatically after the first two samples and looked so strange.
I have try the code "--p-no-golay-error-correction" and "--p-rev-comp-barcodes" to replace and I failed (I put my attempts and errors below).
Attempt1:
I think Qiita provides counts per sample, so you could check that to see if the sample counts are truly wrong, or if some samples simply have more reads inside of them.
You could also extract and summarize barcode counts inside the barcodes.fastq.gz file to see how they compare, but that's a lot more manual processing.
Thanks for your reply,
I have checked the counts provided by Qiita, and they are truly more than that in our outcomes.
Also, I compared the number of lines contained in barcodes.fastq.gz, forward.fastq.gz and reverse.fastq.gz. The three raw files have the same number of lines.
I am not sure if I get what you mean. Could you tell me what should I do next?
Anxious anticipates your reply.
Qiita reports that both rev_comp_barcode and rev_comp_mapping_barcodes was specified. What happens if you rerun using --p-rev-comp-barcodes and --p-rev-comp-mapping-barcodes?
Thank you very much for your time and explanation.
The problem has been well solved when I rerun using --p-rev-comp-barcodes and --p-rev-comp-mapping-barcodes.
However, I still have a confusion. Could you please tell me how can I know whether there have rev_comp_barcode and rev_comp_mapping_barcodes as shown in your reply.
Sincere greetings.
