GFF與GTF格式怎么轉(zhuǎn)換

這篇“GFF與GTF格式怎么轉(zhuǎn)換”文章的知識(shí)點(diǎn)大部分人都不太理解，所以小編給大家總結(jié)了以下內(nèi)容，內(nèi)容詳細(xì)，步驟清晰，具有一定的借鑒價(jià)值，希望大家閱讀完這篇文章能有所收獲，下面我們一起來(lái)看看這篇“GFF與GTF格式怎么轉(zhuǎn)換”文章吧。

GFF主要是用來(lái)注釋基因組，格式如下：

J15     glean   gene    25308430        25309140        .       +       .       ID=Gglean067903; status=novel;
J15     glean   mRNA    25308430        25309140        0.976879        +       .       ID=Gglean067903-TA; Parent=Gglean067903; status=novel;
J15     glean   CDS     25308430        25308501        .       +       0       Parent=Gglean067903-TA;
J15     glean   CDS     25308646        25309140        .       +       0       Parent=Gglean067903-TA;
J15     glean   gene    126763  129003  .       +       .       ID=Gglean075841; status=novel;
J15     glean   mRNA    126763  129003  1       +       .       ID=Gglean075841-TA; Parent=Gglean075841; status=novel;
J15     glean   CDS     126763  126973  .       +       0       Parent=Gglean075841-TA;
J15     glean   CDS     127285  127628  .       +       2       Parent=Gglean075841-TA;
J15     glean   CDS     127719  127854  .       +       0       Parent=Gglean075841-TA;
J15     glean   CDS     128049  128185  .       +       2       Parent=Gglean075841-TA;

而GTF主要是用來(lái)對(duì)基因進(jìn)行注釋，格式如下：

J01     glean   CDS     6976    7317    .       +       0       transcript_id "Gglean025939-TA"; gene_id "Gglean025939";
J01     glean   CDS     7912    8162    .       +       0       transcript_id "Gglean025939-TA"; gene_id "Gglean025939";
J01     glean   CDS     8245    8413    .       +       1       transcript_id "Gglean025939-TA"; gene_id "Gglean025939";
J01     glean   CDS     8479    8790    .       +       0       transcript_id "Gglean025939-TA"; gene_id "Gglean025939";
J01     glean   CDS     9444    9708    .       -       2       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";
J01     glean   CDS     9778    9935    .       -       0       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";
J01     glean   CDS     10012   10216   .       -       2       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";
J01     glean   CDS     10299   10754   .       -       2       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";
J01     glean   CDS     10838   10926   .       -       0       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";
J01     glean   CDS     11015   11082   .       -       1       transcript_id "Gglean025954-TA"; gene_id "Gglean025954";

轉(zhuǎn)換：

Cufflinks里面的工具gffread可以直接實(shí)現(xiàn)GFF與GTF的互相轉(zhuǎn)換。

#gff2gtfgffread my.gff3 -T -o my.gtf#gtf2gffgffread merged.gtf -o- > merged.gff3

以上就是關(guān)于“GFF與GTF格式怎么轉(zhuǎn)換”這篇文章的內(nèi)容，相信大家都有了一定的了解，希望小編分享的內(nèi)容對(duì)大家有幫助，若想了解更多相關(guān)的知識(shí)內(nèi)容，請(qǐng)關(guān)注億速云行業(yè)資訊頻道。