评估基因组组装计算一些用于评估基因组组装的统计数据。依赖项：numpy matplotlib

用法：


import sys

from matplotlib import pyplot

from stats import AssemblyStatistics



# the input contig file, in FASTA format.

inputFile = sys.argv[1]

out = AssemblyStatistics(inputFile)



# L50 of the assembly

l50 = out.L50()



# N50 of the assembly

n50 = out.N50()



# size of the assembly

# ...