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

t2g.py: fool proof output format

parent fda70fd4
No related branches found
No related tags found
No related merge requests found
...@@ -3,6 +3,7 @@ import os ...@@ -3,6 +3,7 @@ import os
import gffutils import gffutils
import argparse import argparse
def validate_file(f): def validate_file(f):
if not os.path.exists(f): if not os.path.exists(f):
# Argparse uses the ArgumentTypeError to give a rejection message like: # Argparse uses the ArgumentTypeError to give a rejection message like:
...@@ -32,10 +33,15 @@ if __name__ == "__main__": ...@@ -32,10 +33,15 @@ if __name__ == "__main__":
with open("t2g.txt", "w") as t2g: with open("t2g.txt", "w") as t2g:
for gene in db.all_features(): for gene in db.all_features():
for transcript in db.children( for transcript in db.children(
gene, featuretype='transcript', order_by='start'): gene, featuretype='transcript', order_by='start'
):
t2g_line = str(transcript["transcript_id"][0]) + \
"\t" + \
str(gene["gene_id"][0])
t2g_line = t2g_line.split("\t")
t2g.write( t2g.write(
str(transcript["transcript_id"][0]) + str(t2g_line[0].split(".")[0]) +
"\t" + "\t" +
str(gene["gene_id"][0]) + str(t2g_line[1].split(".")[0]) +
"\n" "\n"
) )
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment