Question: how to change gtf to fasta file of cufflink, gffread
gravatar for worarado.kan
13 days ago by
worarado.kan10 wrote:

Dear all,

I would like to change gtf file to fasta file. In gffread, if I need only TCONS_00033653 and gene_id "XLOC_028018 ID. what is command I should use?

My gtf file is

v.corymbosum_GDV_reftransV1_0039461 Cufflinks   exon    1171    1309    .-  .   gene_id "XLOC_028018"; transcript_id "TCONS_00033653"; exon_number "2"; oId "CUFF.29927.1"; tss_id "TSS29077";

Actually, I used gffread -w ECvsEL_transcript.fasta -g v.corymbosum_GDV_reftransV1.fasta EC_vs_EL.gtf

But I get only TCONS_00033653 ID

Thank you so much kanjana

ADD COMMENTlink modified 11 days ago by kristoffer.vittingseerup2.4k • written 13 days ago by worarado.kan10

If I understand this, are you referring to FASTA file header when you mentioned

But I get only TCONS_00033653 ID

ADD REPLYlink modified 13 days ago • written 13 days ago by Sej Modha4.5k

Yes I mean head of fasta file because next step I have to joint fasta file with genes.exp.cuffdiff file, so I need the same ID

Could you help how to use command? sorry Im new for this analysis

Thank you so much Kanjana

ADD REPLYlink written 12 days ago by worarado.kan10
gravatar for kristoffer.vittingseerup
11 days ago by
European Union
kristoffer.vittingseerup2.4k wrote:

For most genes you have quantified multiple transcripts. The expression of those transcripts have then been summed to get the gene expression - which is then used to o the gene differential expression. This also means that there is no single sequence which will represent the gene you have been analysing meaning it cannot be done.

What do you need the gene-sequence for? Maybe we can figure out an alternative approach if you tell us why you want to use it for :-)

Cheers Kristoffer

ADD COMMENTlink written 11 days ago by kristoffer.vittingseerup2.4k
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: 1535 users visited in the last hour