1. contig N50 的定义

基因组的统计信息包含GC含量，N50等等，这里我们计算N50的算法：
N50是指一个基因组所有的contig，按照长度从大到小排列，一直长度加和，一直到某条contig达到该基因组总共长度的1/2时，那么这条contig的长度，就是N50的值。N50可以作为评价一个基因组拼接好坏的条件之一。

2. 脚本实现

1. N50 计算：

# 以sequence.fasta 为列，来计算N50的值
# 基因组拼接时，contig是按照长度从大到小排列的
f = open('sequence.fasta', 'r').readlines()
total_length = 0
for i in f:if not i.startswith('>'):seq_length = len(i.strip())total_length += seq_lengthvariable_length = 0
for i in f:if not i.startswith('>'):seq_length = len(i.strip())variable_length += seq_lengthif variable_length >= total_length / 2:print('N50 : %s'%seq_length)          #输出N50的值break

2. GC 含量计算

# 以sequence.fasta 为列，来计算N50的值f = open('sequence.fasta', 'r').readlines()
total_length = 0
total_seq = ''
for i in f:if not i.startswith('>'):seq_length = len(i.strip())total_length += seq_lengthtotal_seq += i.strip()GC_number = 0
for j in total_seq:if j == 'C' or j == 'G':GC_number += 1GC_percent = GC_number / total_length     # 输出为GC含量
print(f'GC percent : {GC_percent}')