Question: rmdup - remove PCR duplicates from ChIP-seq data
0
gravatar for david.fear
4.6 years ago by
david.fear • 10
European Union
david.fear • 10 wrote:

I have aligned my chip-seq data with bow tie and wanted to remove PCR duplicates. I ran rmdup in SAM tools on the aligned file but the number or reads I got back (i.e.. 8787256 / 9006759 = 0.9756 in library) was more than the number of single aligned reads I started with? I assume these numbers are number of unique reads vs total reads tested, is this correct?

 

Many thanks

Dave

samtools chip-seq • 4.9k views
ADD COMMENT • link • modified 4.6 years ago by Jennifer Hillman Jackson ♦ 25k • written 4.6 years ago by david.fear • 10
1
gravatar for Jennifer Hillman Jackson
4.6 years ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hi,

The tool you are looking for is "NGS: Picard (beta) -> Mark Duplicate reads".

Thanks! Jen, Galaxy team

ADD COMMENT • link written 4.6 years ago by Jennifer Hillman Jackson ♦ 25k
0
gravatar for david.fear
4.6 years ago by
david.fear • 10
European Union
david.fear • 10 wrote:

Hi,

Thanks for your reply, how is this different from rmdup in SAM tools?

Dave

 

ADD COMMENT • link written 4.6 years ago by david.fear • 10
0
gravatar for Jennifer Hillman Jackson
4.6 years ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

There is much discussion about the use these two tools. This is a good place to start: http://seqanswers.com/forums/showthread.php?t=5424

Best, Jen, Galaxy team

ADD COMMENT • link written 4.6 years ago by Jennifer Hillman Jackson ♦ 25k
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 16.09
Traffic: 171 users visited in the last hour