Note: The output file for gene or transcript coordinates should always on 1-basis.
Script for manipulating GFF3 annotation file using pysam