Question: How to demultiplex 16S PE Miseq reads?
0
gravatar for Lina F
3.2 years ago by
Lina F160
Boston, MA
Lina F160 wrote:

Hi all,

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!

Cheers,

~Lina

ADD COMMENTlink written 3.2 years ago by Lina F160

Seen this? http://qiime.org/tutorials/processing_illumina_data.html This looks like it explains how to work with Qiime mapping stuff and demultiplexing.

ADD REPLYlink written 3.2 years ago by WouterDeCoster41k

It would be simple for the core that ran your samples to demultiplex this data. At least ask to see if they would do it.

ADD REPLYlink written 3.2 years ago by genomax71k
Please log in to add an answer.

Help
Access

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