python 中 将fastq数据的碱基质量值转换为数值

 

001、

def fastq_quality_to_phred(fastq_file):
    with open(fastq_file, 'r') as f:
        for line in f:
            if line.startswith('+'):
                continue  # 跳过加号行
            quality_line = next(f)  # 下一行为质量值行
            phred_scores = [ord(q) - 33 for q in quality_line.rstrip()]  # 将ASCII转换为Phred分数
            print(phred_scores)
 
# 使用示例
fastq_quality_to_phred('example.fastq')