Question: diffrent cufflinks output results with and without a Gtf feference file
gravatar for hana
4.8 years ago by
hana170 wrote:


I am currently working on RNA-seq using cufflinks, I am confused about the result of cufflinks when I used  –G and –g option and run cufflinks without gtf 

I found different number of genes and transcripts from my output results:

for the -G mode, I found 17,876 gens and 32,568 transcripts.

for the -g mode,I found 25,554 gens and 33,453 transcripts

without using gtf references file, I found 18,342 gees and 34,871 transcripts

Also I observed several new transcripts not listed in the result when I used -G mode where only used the known transcript. I also saw that some of known transcript from -G mode is not in the result using -g mode.

which result is correct for reporting the number of genes and transcripts?

Many thanks for your suggestions

rna-seq • 1.2k views
ADD COMMENTlink modified 4.8 years ago by Devon Ryan90k • written 4.8 years ago by hana170
gravatar for Devon Ryan
4.8 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

I have no idea why you find it confusing that these different methods produce different results, that's the expectation. Simply thinking about how this works would suggest that the -g output is more likely to match what's biologically happening.

ADD COMMENTlink written 4.8 years ago by Devon Ryan90k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1044 users visited in the last hour