This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can a duplicate read be interchromosomal?

Hi all,

I've seen many Q/As on here and elsewhere that state that the main advantage of Picard Mark Duplicates over samtools rmdup is that the former removes interchromosomal duplicates while the latter does not.

My question is, if duplicate reads are determined by them having the same 5' start coordinate, how can reads mapping to different chromosomes (interchromosomal) ever be considered duplicates? They would have a different coordinate.

Cheers!

interchromosomal duplicate ngs samtools picard
• 1,650 views • link
written by maxrwjones • maxrwjones

1 answer

If the coordinates point to identical regions (duplicate chromosomal regions are common), yes, it can be considered duplicated

• 0 views • link
written by JC JC

Thanks for your answer :)

This makes sense and I understand there can be duplicated regions of sequence. It was my understanding though that these tools did not detect sequence identity - they simply flag up reads if the 5' mapped coordinate is the same as that of another read. Even duplicated regions on different chromosomes would have different coordinates... they will be millions of bases apart in the concatenated genomic reference.

But maybe the tools are sorting by both coordinates and sequence?

• 0 views • link
written by maxrwjones • maxrwjones

Log in to answer this question.

More posts like this