Question: How Does Cufflinks Compute The Fpkm For Reads That Do Not Map In The Proper Pair?
3
gravatar for le456900
5.7 years ago by
le45690040
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.5k views
ADD COMMENTlink modified 5.6 years ago by vj400 • written 5.7 years ago by le45690040
0
gravatar for Phil S.
5.7 years ago by
Phil S.660
Stuttgart, Germany
Phil S.660 wrote:

hi,

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.7 years ago by Phil S.660
0
gravatar for vj
5.6 years ago by
vj400
UK
vj400 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.6 years ago by vj400
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 2.3.0
Traffic: 1191 users visited in the last hour