使用Cufflinks里面的工具gffread
執(zhí)行命令：?gffread **.gff?-T -o **.gtf? ?
??????????????????? gffread merged.gtf -o- > merged.gff3
cufflinks詳細(xì)用法參考：陳連福的生信博客
問題挺好：
junction：剪接接頭
博主師兄您好捡硅， 我剛開始學(xué)習(xí)RNA-seq數(shù)據(jù)分析措近，樓主的文章給我收獲很多批幌，謝謝翎猛。 有幾個問題想請教拘荡。
我根據(jù)12年Nature Protocol 的文章臼节，用Tophat和cufflinks處理數(shù)據(jù)。
1. 我的數(shù)據(jù)是細(xì)菌珊皿，基因組是剛測序完的草圖网缝。自己根據(jù)網(wǎng)上的格式介紹嘗試做了一個.gtf 的文件，做Tophat時蟋定，發(fā)現(xiàn)警告“Warning: TopHat did not find any junctions in GTF file”粉臊。不知是.gtf 弄錯了，還是細(xì)菌本來就沒有junction.
我的gtf是這樣的：
scaffold1majorCDS933254.–0gene_id “orf00001”; transcript_id “transcript_orf00001”;
scaffold1majorexon933254.–0gene_id “orf00001”; transcript_id “transcript_orf00001”;
scaffold1majorCDS933254.–0gene_id “orf00001”; transcript_id “transcript_orf00001”;
scaffold1majorstart_codon9395.–0gene_id “orf00001”; transcript_id “transcript_orf00001”;
scaffold1majorstop_codon32523254.–0gene_id “orf00001”; transcript_id “transcript_orf00001”;
2. Cufflinks的結(jié)果里只有FPKM驶兜，如果想要RPKM用什么軟件做呢扼仲？還是cufflinks能算RPKM 只是我沒找到。我看了manual抄淑，找不到北屠凶。
3. 我每個條件只有一組RNAseq， 沒做生物重復(fù)肆资。cuffdiff 能做出結(jié)果矗愧，但是這樣的結(jié)果能用嗎? 有沒有其他軟件能做？ （你文中說“結(jié)果中沒有差異表達(dá)基因”）
1.細(xì)菌本來就沒有junction郑原。junction是對于真核生物含有內(nèi)含子的基因來說的唉韭。
2.Cufflinks應(yīng)該是不能用于計算RPKM的夜涕，現(xiàn)在公認(rèn)的是FPKM，RPKM算是過時的属愤。
3.沒有生物重復(fù)女器，則需要使用cufflinks 2.0版本的。高版本貌似必須有重復(fù)春塌。
