SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Demultiplexing using variable-length barcodes id0 Bioinformatics 1 08-10-2013 02:46 PM
Demultiplexing in-line barcodes saikumarkv Illumina/Solexa 6 03-28-2012 02:07 PM
Custom barcodes AND Paired end seq readlength GenomicIBK SOLiD 3 11-20-2011 06:02 PM
Splitting NuGen barcodes from paired-end sequences ester Bioinformatics 10 10-03-2011 09:36 AM
Library prep Solid4 paired-end BARCODES KNS General 0 08-04-2010 05:26 AM

Reply
 
Thread Tools
Old 04-23-2015, 10:43 AM   #1
ddelalca
Junior Member
 
Location: Philadelphia

Join Date: Apr 2015
Posts: 2
Default Demultiplexing combinatorial barcodes from paired end data

Hey, I've been trying to figure out how to demultiplex paired-end 250 miseq data that has three different barcodes.

The samples all come from 96 well plates. I have two 5' barcodes representing plate number as well as plate column. These two barcodes are linked via a common adapter sequence. I then have one 3' barcode that represents the plate column.

I am fairly sure I would be able to demultiplexing the data in 2 rounds in which I first demultiplex based upon each plate/column combination, and then demultiplex the result based upon the 3' barcode.

This seems fairly cumbersome as I have 2880 unique samples. Is there some sort of software in which I can provide simply the three unique barcodes, and then demultiplex in one step?

Also, I usually use the A5Pipeline (https://code.google.com/p/ngopt/wiki/A5PipelineREADME) to error correct sequence reads, remove Illumina adapter sequences, and assemble the paired-reads. Would I be able to assemble the reads using A5 and then demultiplex the data, or do I have to run the pipeline on each individual sample after first demultiplexing. Thanks!
ddelalca is offline   Reply With Quote
Old 04-27-2015, 06:01 AM   #2
ddelalca
Junior Member
 
Location: Philadelphia

Join Date: Apr 2015
Posts: 2
Default

Does anyone have experience with combinatorial barcodes?
ddelalca is offline   Reply With Quote
Old 04-27-2015, 06:11 AM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,091
Default

How long are these barcodes? It may be possible to write some custom code to do this in one step but without seeing what you reads look like (or a good schematic) no one is going to be able to help. It may be best to find someone locally who can work with you interactively.

If the A5Pipeline changes the context of the barcodes it will make things complicated. You may have to run the pipeline post demultiplexing.
GenoMax is offline   Reply With Quote
Old 04-27-2015, 08:07 AM   #4
JackieBadger
Senior Member
 
Location: Halifax, Nova Scotia

Join Date: Mar 2009
Posts: 380
Default

This will do it https://code.google.com/p/jmhc/
JackieBadger is offline   Reply With Quote
Old 04-27-2015, 08:30 AM   #5
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,091
Default

Quote:
Originally Posted by JackieBadger View Post
Can you import fastq files? Description states that it operates on fasta files (internally?). I don't think this will do 3 barcodes in a single pass.
GenoMax is offline   Reply With Quote
Old 05-01-2015, 01:42 AM   #6
relipmoc
Member
 
Location: Los Angeles, CA

Join Date: Jul 2011
Posts: 58
Default use skewer

Hi ddelalca,
You may use skewer for this purpose. An example command line is as follows:
Code:
skewer -x V3-forward.txt -y V4-reverse.txt --mode ap -t 4 ampliconReads-pair1.fastq ampliconReads-pair2.fastq assigned/
If you only use some of the combinations, you may specify the matrix by -M as well.
Good luck!
Attached Files
File Type: txt V3-forward.txt (312 Bytes, 31 views)
File Type: txt V4-reverse.txt (483 Bytes, 13 views)
File Type: txt matrix.txt (302 Bytes, 27 views)
relipmoc is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:43 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO