SEQanswers (
-   RNA Sequencing (
-   -   Filter cuffcompare by FPKM (

SS00 12-08-2014 03:47 AM

Filter cuffcompare by FPKM
How do I filter out low cuffcompare FPKM values out of by combined.gtf file?

I am trying to use the lncrScan workflow for novel long non-coding RNA from total RNA samples. I assembled RABT assembly with cufflinks using combined gtf from UCSC, NCBI and emsembl and then ran cuffcompare on the resulting assembly.

However before I pull out my sequences from the gtf in FASTA format, I want t filter out all transcripts below FPKM of 2. I suppose there could be an awk script for this, but in my combined.gtf there is no FPKM value listed.

Any ideas?

Thanks for the help!

pengchy 04-26-2015 01:03 AM

please take a look at filtercuffcompare from seqscripts:

I have searched around found nobody has used this script here. I have tested this script, filtered plenty of the combined gtf. The rationality is under checking.

All times are GMT -8. The time now is 06:19 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.