Hi, I want to ask how to get the raw counts from the output of cufflinks. One way to do this is to use the fpkm.
raw counts = FPKM * (length of that transcript/1000) * (# of mapped reads / 1e6)
The FPKM and length of transcript are in the cufflinks FPKM Tracking Files. But how about the # of mapped reads?
For instance, we have a foo.bam. samtools view -c (-f|-F) flag foo.bam can do this job but I am not quite which flag should I set when it's single-end or paired-end.