@@ -460,7 +460,7 @@ def main(args=None):
460
460
461
461
fastafile = generate_fasta (args .output ,queryfile ,resultfile )
462
462
smorf_number = int (predicted_smorf_count (fastafile )/ 2 )
463
- summary .append (f'{ smorf_number } smORFs aligned against GMSC in total.\n ' )
463
+ summary .append (f'{ smorf_number } smORFs are aligned against GMSC in total.\n ' )
464
464
465
465
if not args .noquality :
466
466
summary .append (f'# Quality' )
@@ -476,14 +476,14 @@ def main(args=None):
476
476
if not args .notaxonomy :
477
477
summary .append (f'# Taxonomy' )
478
478
annotated_number ,rank_number ,rank_percentage = taxonomy (args ,resultfile ,tmpdir )
479
- summary .append (str (annotated_number )+ '(' + str (round ((1 - rank_percentage ['no rank' ])* 100 ,2 ))+ '%) aligned smORFs have taxonomy annotation.' )
480
- summary .append (str (rank_number ['kingdom' ])+ '(' + str (round (rank_percentage ['kingdom' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on kingdom.' )
481
- summary .append (str (rank_number ['phylum' ])+ '(' + str (round (rank_percentage ['phylum' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on phylum.' )
482
- summary .append (str (rank_number ['class' ])+ '(' + str (round (rank_percentage ['class' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on class.' )
483
- summary .append (str (rank_number ['order' ])+ '(' + str (round (rank_percentage ['order' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on order.' )
484
- summary .append (str (rank_number ['family' ])+ '(' + str (round (rank_percentage ['family' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on family.' )
485
- summary .append (str (rank_number ['genus' ])+ '(' + str (round (rank_percentage ['genus' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on genus.' )
486
- summary .append (str (rank_number ['species' ])+ '(' + str (round (rank_percentage ['species' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on species.' )
479
+ summary .append (str (annotated_number )+ ' (' + str (round ((1 - rank_percentage ['no rank' ])* 100 ,2 ))+ '%) aligned smORFs have taxonomy annotation.' )
480
+ summary .append (str (rank_number ['kingdom' ])+ ' (' + str (round (rank_percentage ['kingdom' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on kingdom.' )
481
+ summary .append (str (rank_number ['phylum' ])+ ' (' + str (round (rank_percentage ['phylum' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on phylum.' )
482
+ summary .append (str (rank_number ['class' ])+ ' (' + str (round (rank_percentage ['class' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on class.' )
483
+ summary .append (str (rank_number ['order' ])+ ' (' + str (round (rank_percentage ['order' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on order.' )
484
+ summary .append (str (rank_number ['family' ])+ ' (' + str (round (rank_percentage ['family' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on family.' )
485
+ summary .append (str (rank_number ['genus' ])+ ' (' + str (round (rank_percentage ['genus' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on genus.' )
486
+ summary .append (str (rank_number ['species' ])+ ' (' + str (round (rank_percentage ['species' ]* 100 ,2 ))+ '%) aligned smORFs are annotated on species.' )
487
487
488
488
with atomic_write (f'{ args .output } /summary.txt' , overwrite = True ) as ofile :
489
489
for s in summary :
0 commit comments