SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Demultiplexing MiSeq run reads. memento Illumina/Solexa 5 09-10-2014 05:24 AM
CLC genomics wokbench and illumina demultiplexing Etherella Bioinformatics 8 06-03-2014 07:48 AM
duplicate reads in Illumina short, single end reads of RNAseq data inbarpl Bioinformatics 4 05-22-2012 08:36 AM
Demultiplexing Casava 1.8 reads Metz Illumina/Solexa 6 01-26-2012 05:06 AM
Demultiplexing Illumina reads Bardj Bioinformatics 1 02-16-2011 10:17 AM

Reply
 
Thread Tools
Old 02-21-2014, 09:52 AM   #1
Bioinform
Member
 
Location: US

Join Date: May 2013
Posts: 17
Default Demultiplexing illumina reads

Hi there,

I have some basic questions. It would be really great if any one could help me out.

1. I wanted to know different ways of demultiplexing. I have worked with demultiplex.pl script targeting barcodes.
Is there any scripts or ways targeting sequence.
I want both the F and R reads in one file the one i am using gives me seperate file (demultiplex.pl).

2. After demultiplexing i get different number of reads in F and R files (i suppose it to be equal in number). What would be the threshold of difference in the number of reads in F and R files?

Thanks
Bioinform is offline   Reply With Quote
Old 02-21-2014, 10:03 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

I assume these are in-line barcodes because otherwise illumina's software takes care of demultiplexing. Not sure why you are getting unequal number of F/R reads.

Is there a specific reason you want F/R reads in one file?
GenoMax is offline   Reply With Quote
Old 02-21-2014, 11:08 AM   #3
JackieBadger
Senior Member
 
Location: Halifax, Nova Scotia

Join Date: Mar 2009
Posts: 381
Default

You shouldn't have different numbers of reads in R1 and R2 files. Maybe your demultiplex script is throwing some sequences out based on quality or Ns?

Are R1 and R2 derived from a targeted amplicon, and should they overlap?
If so you should stitch them together first. Then you can use jMHC to demultiplex (after conversion to fast).

A little more information on your system would help in getting a clearer answer.
JackieBadger is offline   Reply With Quote
Old 02-23-2014, 06:08 PM   #4
Bioinform
Member
 
Location: US

Join Date: May 2013
Posts: 17
Default

Hi Guys thanks for the quick reeply.

GenoMax

my demultiplex script is not giving me equal number of reads, they pull F and R in different file so they may not be equal.
I think there are demultiplex script that pulls both the reads in one file, pulls F and automatically pulls R, so that the F and R reads become paired.

JackieBadger
I was guessing so too.

I have miseq data. I pooled 8 libraries and send them for sequencing. Since the barcodes are not the illumina kit barcodes i told them to send the files and i would demultiplex them.

With CLC i get similar number of F and R reads (paired in one file).
But with the demultiplex.pl script i get different number of F and R reads separately.
I was looking for demultiplex script that pulls both the reads in one file, pulls F and automatically pulls R targeting the sequence, and F and R reads get paired.

Thanks
Bioinform is offline   Reply With Quote
Old 02-25-2014, 09:37 AM   #5
Bioinform
Member
 
Location: US

Join Date: May 2013
Posts: 17
Default

Hi,

Does anyone has perl or pythogn script that can pull out Reads (Forward) from R1 file and corresponding pair (Reverse) from R2 file. CLC workbench does this but want some script that works alike.

Many thanks
Bioinform 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 12:40 AM.


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