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')