cufflinks and non-uniquely mapped reads

clariet · 05-07-2010, 07:40 AM

Cufflinks appear to take not only uniquely mapped reads. I only filtered out low quality mapped reads and ran cufflinks without any problem. But how cufflinks take care of those reads mapped to multiple regions? Does it only pick the region with highest mapping score? Anyone has ideas?

I am wondering how the results will differ if I use uniquely mapped reads only.

Any input will be highly appreciated.
Clariet

lpachter · 05-08-2010, 12:13 PM

The correct way to deal with reads that map to multiple places is explained in the paper by Li et al:
http://bioinformatics.oxfordjournals...tract/26/4/493

Mathematically, the ambiguity that arises from multiple hits in alignment, is no different than the ambiguity that arises when reads can belong to multiple isoforms at one location- in fact the statistical deconvolution is identical.

We plan to extend Cufflinks to correctly solve this problem in the near future.
Currently, multiple hits are assigned equal probabilities- details are in the supplement to the Cufflinks paper:
http://www.nature.com/nbt/journal/va...bt.1621-S1.pdf

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Tophat uniquely mapped reads	mrfox	Bioinformatics	2	05-23-2013 06:58 AM
How to get uniquely mapped reads from Tophat	subeet	Bioinformatics	10	11-28-2012 06:56 AM
not uniquely mapped reads	unidodo	RNA Sequencing	2	04-22-2011 03:07 PM
Uniquely mapped reads with bowtie	mapper	Bioinformatics	2	11-22-2010 11:44 PM
BWA Uniquely Mapped Reads	NF_seq	Bioinformatics	0	09-06-2010 04:32 AM

05-07-2010, 07:40 AM	#1
clariet Member Location: NJ Join Date: Mar 2010 Posts: 18	cufflinks and non-uniquely mapped reads Cufflinks appear to take not only uniquely mapped reads. I only filtered out low quality mapped reads and ran cufflinks without any problem. But how cufflinks take care of those reads mapped to multiple regions? Does it only pick the region with highest mapping score? Anyone has ideas? I am wondering how the results will differ if I use uniquely mapped reads only. Any input will be highly appreciated. Clariet

05-08-2010, 12:13 PM	#2
lpachter Member Location: Berkeley, cA Join Date: Feb 2010 Posts: 40	The correct way to deal with reads that map to multiple places is explained in the paper by Li et al: http://bioinformatics.oxfordjournals...tract/26/4/493 Mathematically, the ambiguity that arises from multiple hits in alignment, is no different than the ambiguity that arises when reads can belong to multiple isoforms at one location- in fact the statistical deconvolution is identical. We plan to extend Cufflinks to correctly solve this problem in the near future. Currently, multiple hits are assigned equal probabilities- details are in the supplement to the Cufflinks paper: http://www.nature.com/nbt/journal/va...bt.1621-S1.pdf