Implement a stand-alone Python program to calculate N50 from the FASTA format contigs file.

Usage:

python myquast.py test_data/contigs_1.fasta \ -o quast_test_output (optinal for output directory) -h (showing usage)

Output:

report.txt

The report should include:

Number of contigs and total length of them.

Length of the largest contig.

N50 (length of a contig, such that all the contigs of at least the same length together cover at least 50% of the assembly).

Recommend you to use the BioPython library to read the FASTA format file, but you can implement the section by yourself.