GFF3提取fasta序列

时间 2020-02-12 标签 gff3 gff 3 提取 fasta 序列

awkweb

awk '{print ">"$9"\n"$10}' file.gff3

Pythonsvg

fr = open('./file.gff3', 'r')
for line in fr.readlines():
    name = '>' + line.strip().split('\t')[8]
    seq = line.strip().split('\t')[9]
    print(name)
    print(seq)