Skip to content
Snippets Groups Projects
Verified Commit bd5b4ace authored by Laurent Modolo's avatar Laurent Modolo
Browse files

kb: change t2g.py to use classical transcript and gene gtf format

parent 192fef09
No related branches found
No related tags found
No related merge requests found
......@@ -16,7 +16,7 @@ if __name__ == "__main__":
help="gtf file", metavar="FILE")
args = parser.parse_args()
db = gffutils.create_db(args.gtf, dbfn = ":memory:", force = True, merge_strategy="merge", disable_infer_transcripts=False, disable_infer_genes=False)
db = gffutils.create_db(args.gtf, dbfn = ":memory:", force = True, merge_strategy="merge", disable_infer_transcripts=True, disable_infer_genes=True)
with open("t2g.txt", "w") as t2g:
for gene in db.all_features():
for transcript in db.children(gene, featuretype='transcript', order_by='start'):
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment