Skip to content

Commit

Permalink
Update mitohifi.py
Browse files Browse the repository at this point in the history
  • Loading branch information
marcelauliano authored Nov 7, 2023
1 parent 9a76a7d commit 2cfc1ef
Showing 1 changed file with 5 additions and 3 deletions.
8 changes: 5 additions & 3 deletions src/mitohifi.py
Original file line number Diff line number Diff line change
Expand Up @@ -410,7 +410,9 @@ def main():
"length(bp)", "number_of_genes", "was_circular\n"]))
f.write("\t".join(["final_mitogenome", all_frameshifts, final_annotation,
str(contig_len), str(num_genes), str(is_circ)+"\n"]))

f.write("\t".join([repr_contig_id, all_frameshifts, final_annotation,
str(contig_len), str(num_genes), str(is_circ)+"\n"]))

# compare list of genes from related mito and final mitogenome
related_mito_annotation = args.g
related_mito_genes = get_genes_list(related_mito_annotation, "genbank")
Expand All @@ -434,9 +436,9 @@ def main():
contigs_stats_files = []
for curr_file in os.listdir('.'):
if curr_file.endswith('.individual.stats'):
# don't skips addition of representative contig, which is the
# skips addition of representative contig, which is the
# same as the final_mitogenome
if curr_file.split('.individual.stats')[0] == repr_contig_id:
if curr_file.split('.individual.stats')[0] != repr_contig_id:
contigs_stats_files.append(curr_file)

with open("contigs_stats.tsv", "a") as outfile:
Expand Down

0 comments on commit 2cfc1ef

Please sign in to comment.