How to calculate all transcripts from cufflinks
1
0
Entering edit mode
7.7 years ago
Paul ★ 1.4k

Dear all,

I have two mouses paired-end sequenced. I mapped with Tophat2 - make Cufflinks - and Cufdiff (one replicates - two conditions) - and CummeRbund. It is very nice graphs and outputs, but I would like to get just one number, which describe me total transcript expression for mouse1 and for mouse2 and compare it.

It is possible to take output form Cufflinks - transcripts.gtf and extract FPKM value for all genes and do sum of all FPKM's? Something like - awk 'OFS="\t"{sum+=$16}END {print sum}' transcripts.gtf.

Thank you so much for explanation and help,

Paul.

RNA-Seq FPKM cufflinks • 2.3k views
ADD COMMENT
1
Entering edit mode
7.6 years ago
Paul ★ 1.4k

Answer is pretty easy - just print the FPKM > 0 from isoforms.fpkm_tracking file.

So it could be:

awk '{if ($10>0) print $0}' isoforms.fpkm_tracking | awk 'OFS="\t" {print $1,$4,$6,$7,$8}' | awk 'END{print NR}'

This tell us how many isoforms we have.

ADD COMMENT

Login before adding your answer.

Traffic: 3284 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6