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.
Thanks!
3 answers
A workaround could be to get the number reads that overlap with an interval via bedtools or bedops something like:
bedtools intersect -c ...
and use the reads file and the transcript interval file. The -c flag can give you counts that overlap.
A workaround could also be to get the number of reads for all genes/transcripts specified in a GTF/GFF file using htseq-count
htseq-count [options] <sam_file> <gff_file>
Log in to answer this question.
-
FPKM to raw count conversion
written by Dominik •
-
filtered tab delimited file with awk
written by yaghoub.amraei •
-
Differences between Cufflinks results and home-made FPKM values
written by nicolas.hipp •
-
Effective length in normalisation for reads counts
written by Xiaokang •
-
Why cufflinks does not report raw read count anymore? Which version of it do report?
written by bibilong1111 •
-
How Does Cufflinks Compute The Fpkm For Reads That Do Not Map In The Proper Pair?
written by le456900 •
-
Which Expression Units To Use, Fpkm Or Rpkm ?
written by biorepine
-
Draw Heatmap Or Do Pca Analysis With Raw Read Counts?
written by camelbbs
-
Stranded Paired-End Data With Rsem
written by Ge •
-
How To Calculate Fpkm (Fragments Per Kilobase Of Exon Per Million Fragments Mapped)
written by User 5402 •
