Question: How Does Cufflinks Compute The Fpkm For Reads That Do Not Map In The Proper Pair?
gravatar for le456900
5.4 years ago by
le45690040 wrote:

I want to use cufflinks in our study. I want to know if read1 of paired reads mapped to chr1, and read2 mapped to chr2, how cufflinks calculate the FPKM value.

cufflinks • 2.4k views
ADD COMMENTlink modified 5.3 years ago by vj390 • written 5.4 years ago by le45690040
gravatar for Phil S.
5.4 years ago by
Phil S.660
Stuttgart, Germany
Phil S.660 wrote:


in their manual they state:

'...the RNA-Seq fragment counts can be used as a measure of relative abundance of transcripts, and Cufflinks measures transcript abundances in Fragments Per Kilobase of exon per Million fragments mapped (FPKM), which is analagous to single-read "RPKM", originally proposed in:....'

so, I guess they just calucalte for each of the fragments (which is here chr1 and chr2) the FPKM, thats it. Since they don't incorporate the paired end info here....

imho, they don't have to as long as those reads map to a certain position in a unique fashion...

ADD COMMENTlink written 5.4 years ago by Phil S.660
gravatar for vj
5.3 years ago by
vj390 wrote:

I have the same question and I went through their website. What I believe from their explanation is that they estimate fragment lengths for paired-end reads from the reads mapped to single-isoform genes. Given that they give the transcript abundances in terms of FPKM, where a fragment comes from two reads that are paired, I believe they would disregard improperly paired reads.

Also if you are using TopHat to align the reads, I think you can use "--no-discordant" option to remove reads where the mates are aligned to different chromosomes.

ADD COMMENTlink written 5.3 years ago by vj390
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1416 users visited in the last hour