https://github.com/fjruizruano/ngs-protocols
Raw File
Tip revision: 39a091d1fa569a7fc717ac73c4b3de07f0a1204d authored by fjruizruano on 03 August 2023, 11:48:27 UTC
adding gfa2fas.py and extract_gfa.py
Tip revision: 39a091d
gff_creator.py
#!/usr/bin/python

from Bio import SeqIO

input = raw_input("Introduce FASTA file: ")

secu = SeqIO.parse(input, "fasta")

output = open(input+".gff", "w")

for sq in secu:
	output.write("%s\tprotein_coding\texon\t1\t%s\t.\t+\t.\tgene_id \042%s\042\n" % (sq.id, str(len(sq.seq)), sq.id))

output.close()

back to top