Question: Paired ends demultiplexing Qiime Illumina
gravatar for carvin900
5.9 years ago by
carvin9000 wrote:


sorry in advance for my beginner level in bioinformatic. I need to do analysis of illumina reads paired end. The starting files i have are: fwd.fastq , reverse.fastq and mapping file.txt. 

In the sequences files, there are several samples with different barcodes (barcodes with different lengths) and different primers (3 different set of primers). I tried so far to:

1) I extracted barcodes with --bc1_len X --bc2_len X (where is x is the length of the barcode). I have barcodes of length 4,5,6,7 and 8 so I did 5 times this command with different lengths

2) I merged the reads1.fastq and reads2.fastq output files from 1) for each barcodes length.I had 1,78% in general of non-merging sequences.I also asked to get a fastjoin.join_barcodes.fastq file.

3) I tried to do with the merged sequences and the fastjoin.join_barcodes.fastq file. But then I get over 1 900 000 reads non associated to barcodes out of 2 500 000 total reads.

If someone could help me I would be really glad!


ADD COMMENTlink modified 5.9 years ago • written 5.9 years ago by carvin9000

Hi Geek_y,

thanks for the rapid answer. I am fairly new in bioinformatic and my research lab is mostly using Qiime, so haven't heard of CASAVA. Also, our local server is based on Qiime. Would it be possible to launch commands from Qiime to CASAVA? Thanks for your help!



ADD REPLYlink written 5.9 years ago by carvin9000

Thanks for your help Geek_y, i found out my mistake. There was a switch of two bases in our primer sequence registered and the primer sequence in the reads. Now all is doing fine. I only have done merging and already Thanks again!

ADD REPLYlink written 5.9 years ago by carvin9000
gravatar for geek_y
5.9 years ago by
geek_y11k wrote:

Why don't you use illumina cassava for demultiplexing the data ? Is it a GBS kind of data where you have internal barcodes other than standard barcodes used for multiplexing ?

ADD COMMENTlink written 5.9 years ago by geek_y11k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1914 users visited in the last hour