SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Splitting NuGen barcodes from paired-end sequences ester Bioinformatics 10 10-03-2011 09:36 AM
paired-end reads mapped to genome.. gene with only one direction of paired-end reads? danwiththeplan Bioinformatics 2 09-22-2011 02:06 AM
Barcode splitter - for paired end reads & specifiy prefix in output file spoonman Bioinformatics 2 08-10-2011 03:54 PM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-04-2011 11:51 PM
Splitting 454 paired reads in a FASTQ file sjackman Bioinformatics 5 09-10-2010 11:09 AM

Reply
 
Thread Tools
Old 10-25-2011, 12:05 PM   #1
pbatzel
Junior Member
 
Location: OR, United States

Join Date: Aug 2011
Posts: 5
Default Illumina1.8 Paired-End Barcode Splitting?

Hi everyone,

Hopefully this is a pretty simple question with a pretty simple answer. I've been using a pipeline for pre-processing Illumina1.8, fresh off the sequencer, data. The pipeline puts sequences with the correct first 5bp (we have a list of 5bp barcodes denoting different samples), including however many mismatches we want to allow, into separate files. We run it first on the /1 half of the data and then on the /2 half and end up with AACCC_1.fq and AACCC_2.fq files for each different barcode.

My question comes in when a given sequence from the /1 side matches a barcode in my list exactly, but the /2 side is not the same as the /1 and doesn't match any barcode. Am I able to safely assume that because the /1 and /2 have the same sequenceID from Illumina and are supposedly read from different ends of the same molecule, that I can put the /2 into the same barcode file as the /1?

My biology background isn't very strong, so sorry ahead of time is this is a simple question. Right now if only one side matches a given barcode, the other side is just thrown out. That would be a terrible waste of PE data if my assumption above is safe. Let me know what you think.
pbatzel is offline   Reply With Quote
Old 10-25-2011, 01:46 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

I guess you are running your own custom barcodes somehow? Like they are embedded in the sequence at both ends?

But in theory, if the reads came from the same flowcell, and lane, and tile, and xy-coordinates, (which is all in the original name given by Illumina to the read)they should be the two ends of the same DNA fragment.
swbarnes2 is offline   Reply With Quote
Old 10-25-2011, 02:08 PM   #3
pbatzel
Junior Member
 
Location: OR, United States

Join Date: Aug 2011
Posts: 5
Default

Hi swbarnes thanks for your reply.

Yes, not knowing the biology behind it, the barcodes are somewhat magical to me (although I suppose a lot of this is). My very rough understanding of it is that we use a specific primer for each different sample and therefore can set our own 3' and 5' "barcodes" or adapters, or whatever else someone might call them. I also should have mentioned this is RNA-Seq stuff, so everything is cDNA.

You're spot on with the illumina "ID" I mentioned being something like lane, title, xy and I am indeed hoping that I can make the assumption that they are from the same DNA fragment if they have the same ID. If they are from the same fragment and one side has a barcode that is correct, I should be able to fix the barcode on the second side and bin it into the correct file.

After I posted this question, a colleague mentioned something about "chimeric molecules". I'm assuming from the PCR step, that would botch up this assumption, but only for a very small % of the reads unless something has gone horribly wrong. Again, I'm very new to all of this and have very little biological background so any info is very helpful.
pbatzel is offline   Reply With Quote
Reply

Tags
barcode, illumina, paired-end, preprocessing

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 11:50 AM.


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