GFF3提取fasta序列

awk

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

Python

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)

猜你喜欢

转载自blog.csdn.net/weixin_40099163/article/details/84282786