SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
duplicate reads removal vasvale Bioinformatics 19 01-08-2015 12:59 AM
duplicate read removal gab0 Bioinformatics 14 10-02-2014 03:20 PM
some basic questions about duplicate removal ? a_mt Bioinformatics 1 06-27-2014 09:14 AM
Trinity-Duplicate removal reema Bioinformatics 2 02-27-2014 01:49 AM
threshold for duplicate removal? mard Bioinformatics 2 03-21-2010 03:45 PM

Reply
 
Thread Tools
Old 05-07-2015, 05:04 PM   #1
DNATECH
Member
 
Location: Davis, CA

Join Date: Mar 2015
Posts: 29
Default HiSeq3000/4000 duplicate removal

Hi All,

Our Illumina FAS mentioned that in underloaded lanes on the patterned flow-cells even PCR-free libraries can generate significant numbers of duplicates (clusters swapping over the "walls" into the next empty nanowell if given enough time?).
http://dnatech.genomecenter.ucdavis....data-download/
Would you know of any software that can analyze/remove duplicates based on the coordinate information in the read headers (without any alignments)?

Thanks in advance,
Lutz

Last edited by DNATECH; 05-07-2015 at 06:01 PM.
DNATECH is offline   Reply With Quote
Old 05-08-2015, 01:11 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

I doubt that software for this exists in the wild since you are asking for a very specific operation based on the location of the cluster for a condition that may occur with a relatively new technology (unless those sites that have had HiSeq X Ten have implemented something for this). Wonder if illumina can include mitigation for that in bcl2fastq itself.
GenoMax is offline   Reply With Quote
Old 05-08-2015, 08:47 AM   #3
DNATECH
Member
 
Location: Davis, CA

Join Date: Mar 2015
Posts: 29
Default

Hi GenoMax,

Thanks! Currently there are no local duplicate metrics in the SAV viewer and no options in bcl2fastq for flagging or filtering for these cases.
People working with X Ten data get their data already aligned and might be able to use picard etc.? However the HiSeq 3000/4000 will be used for all kinds of applications were alignment based analyses are not an option.
DNATECH is offline   Reply With Quote
Old 05-08-2015, 09:12 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

Quote:
Originally Posted by DNATECH View Post
Thanks! Currently there are no local duplicate metrics in the SAV viewer and no options in bcl2fastq for flagging or filtering for these cases.
People working with X Ten data get their data already aligned and might be able to use picard etc.? However the HiSeq 3000/4000 will be used for all kinds of applications were alignment based analyses are not an option.
There are tools to address duplicates (e.g. dedupe.sh from BBMap) but they are not cluster location/adjacency aware. Picard should be in this category too unless it has been modified.

I was wondering if this would be something illumina should implement in bcl2fastq software, specially if duplicates will be a problem with less than perfect libraries.
GenoMax is offline   Reply With Quote
Old 05-08-2015, 10:33 AM   #5
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 501
Default

They should merge duplicate reads in contiguous wells for extra-high quality!
__________________
Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com
SNPsaurus is offline   Reply With Quote
Old 05-11-2015, 07:13 AM   #6
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,315
Default

These would look like "optical duplicates", no?
I remember a recent thread where these were found in some MiSeq runs at high frequency.
--
Phillip
pmiguel is offline   Reply With Quote
Old 05-11-2015, 07:21 AM   #7
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,315
Default

Quote:
Originally Posted by pmiguel View Post
These would look like "optical duplicates", no?
I remember a recent thread where these were found in some MiSeq runs at high frequency.
--
Phillip
Ah, the thread is here:

http://seqanswers.com/forums/showthread.php?t=50115

Looks like "Picard MarkDuplicates" was used.

--
Phillip
pmiguel 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 07:18 PM.


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