怎么用Python對gff3格式進行處理

這篇文章主要介紹“怎么用Python對gff3格式進行處理”的相關(guān)知識，小編通過實際案例向大家展示操作過程，操作方法簡單快捷，實用性強，希望這篇“怎么用Python對gff3格式進行處理”文章能幫助大家解決問題。

1.平常從下載的植物基因組注釋文件，以gff3格式為例

      從JGI上下載的玉米基因組的注釋文件（gff3格式），第一例通常為基因的定位信息。1則表示位于玉米的1號染色體上，第二例表示注釋的版本信息，第三列通常為gene，mRNA，CDS等信息，同時一個基因可能對應(yīng)多個mRNA，對生物有些了解的也知道，一個mRNA即是一個轉(zhuǎn)錄信息， 這個和注釋過后的序列文件也是一一對應(yīng)的。第四列和第五列分別為基因的第三列信息在染色體上的物理位置。第七列則表示基因位于正鏈上還是負鏈上。第八列是相位信息。第九列則是基因注釋的一些ID信息。而本次數(shù)據(jù)處理主要是提取第一列、第三列、第四列、第五列和第九列的信息。

2.利用程序處理后的結(jié)果

      處理后的格式一共有五列，第一列為基因所在的染色體上，第二列是利用gene的起始位點和終止位點進行排序的邏輯順序而生成的新的ID信息，第三列和第四列就是基因的起始位點和終止位點了。第五列就是從原始的注釋信息第九列提取出來的，必須個序列文件一一對應(yīng)才行。下面直接上代碼。

3.代碼信息

#!usr/bin/python
import re,io
from operator import itemgetter
input_file = io.open('Zmays_284_Ensembl-18_2010-01-MaizeSequence.gene.gff3','r',encoding='UTF-8')
# 基因的注釋信息,GFF3格式的文件
out_file = open('Zm.newid.gff', 'w', encoding='UTF-8')
# 輸出文件的名字
list_two = []
chr_name = []
de_list = ('#','M','P','s') # 需要修改
for line in input_file:
    if line.startswith(de_list):
        continue
    list_one = line.strip().split()
    if list_one[2] == 'mRNA':
        # gene_id = list_one[8].split(';')[2] # 需要修改
        gene_id = list_one[8]
        gene_id = ''.join(re.findall(r'pacid=(.+?);longest',gene_id)) # 需要修改
		# 獲取gene的id信息
        list_one[0] = re.sub(r'\D',"",list_one[0])
        # list_two.append(gene_na_st_end)
        list_two.append((int(list_one[0]), int(list_one[3]), int(list_one[4]), int(gene_id)))
        chr_name.append(int(list_one[0]))
        # print (gene_id)
    else:
        continue
chr_name = list(set(chr_name))
chr_name.sort()

number = 0
list_thrre = sorted(list_two,key = itemgetter(0,1,2))
next_chr = 0
for i in list_thrre:
    new_i = "\t".join('%s' %id for id in i)
    # print (new_i)
    lp = str(new_i).strip().split()
    # chr_id = re.sub('\[',"",lp[0])
    if str(lp[0])== str(chr_name[next_chr]):
        number = number + 1
    else:
        number = 1
        next_chr = next_chr + 1
    # newid = "Zm"+'%02'%lp[0]+'G'+'%05'%number
    newid = "Zm"+str(lp[0]).zfill(2)+"G"+str(number).zfill(5) # 需要修改 
    print (newid)
    out_file.write('Zm'+str(lp[0])+"\t"+newid+"\t"+str(lp[1])+"\t"+str(lp[2])+"\t"+str(lp[3])+'\n')
input_file.close()
# make by ligaojie from North China University of Technology

關(guān)于“怎么用Python對gff3格式進行處理”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識，可以關(guān)注億速云行業(yè)資訊頻道，小編每天都會為大家更新不同的知識點。