Extract single transcript sequences from genome assembly

Question

[kaz0816]

Hi,

I want to extract one transcript per gene from the genome assembly and its GFF3 file registered in Refseq. And I want to extract a transcript that has a longest or a certain frame. Anyone know of an efficient method?

I am looking for this method to do ortholog comparisons.

Best regards,

[Bastien Hervé]

Python and biopython can help you on that

https://biopython.org/wiki/GFF_Parsing

https://pypi.org/project/gff3/

https://pythonhosted.org/gffutils/