Subset a gff3 file

This can be accomplished with the package "majiq_tools", which can be installed from here: https://bitbucket.org/biociphers/majiq-tools/src/master/

import majiq_tools as mt

gff3 = mt.gff3.load_gff3("path/to/gff3/file")

gene_ids = ["gene:ENSG00000006652", "gene:ENSG00000087460"]

gff3_subsetted = mt.gff3.subset_genes(df_gff3=gff3,gene_id=gene_ids)

mt.gff3.save_gff3(gff3_subsetted, "path/to/new/gff3/file")