diff --git a/make_mCherry_transgene.py b/make_mCherry_transgene.py index 7cca3308fbb84a24fc6b8c7786cfae9ca5b39fa8..165c3e9c9ac95b15bfa6caeb4669d3182bef2b1d 100755 --- a/make_mCherry_transgene.py +++ b/make_mCherry_transgene.py @@ -50,7 +50,7 @@ def main(): with fs.open("mCherry_transcript.fa", "w") as transcript_file: transcript_file.write(f">mCherry\n{mCherry_seq}\n") with fs.open("mCherry.gtf", "w") as transcript_gtf_file: - annotations = "gene_id \"mCherry\"; transcript_id \"mCherry\"; gene_biotype \"mCherry\";" + annotations = "gene_id \"mCherry\"; transcript_id \"mCherry\"; gene_biotype \"protein_coding\";" transcript_gtf_file.write("\t".join([ "mCherry", "local", "transcript", "1", str(len(mCherry_seq)), "0", "+", ".", annotations])) @@ -77,7 +77,7 @@ def main(): "mCherry_transgene", "local", "exon", str(start_in_transgene), str(end_in_transgene), "0", "+", ".", annotations])) transgene_gtf_file.write("\n") - annotations = f"exon_id \"mCherry.e{exon_number}\"; exon_number \"{exon_number}\"; gene_id \"mCherry\"; transcript_id \"mCherry\"; gene_biotype \"mCherry\";" + annotations = f"exon_id \"mCherry.e{exon_number}\"; exon_number \"{exon_number}\"; gene_id \"mCherry\"; transcript_id \"mCherry\"; gene_biotype \"protein_coding\";" transcript_gtf_file.write("\t".join([ "mCherry", "local", "exon", str(start_in_transcript), str(end_in_transcript), "0", "+", ".", annotations]))