fastq文件质量评估python代码

以下是一个基本的fastq文件质量评估的Python代码：

from Bio import SeqIO
import numpy as np

def quality_scores(file_path):
    # 从fastq文件读取序列和质量分数
    sequences = []
    quality_scores = []
    for record in SeqIO.parse(file_path, "fastq"):
        sequences.append(record.seq)
        quality_scores.append(record.letter_annotations["phred_quality"])
    
    # 计算平均质量分数、标准差和中位数
    mean_quality = np.mean([np.mean(scores) for scores in quality_scores])
    std_quality = np.mean([np.std(scores) for scores in quality_scores])
    median_quality = np.median([np.median(scores) for scores in quality_scores])
    
    # 返回结果
    return mean_quality, std_quality, median_quality

这个代码使用了Biopython库来读取fastq文件，并计算了平均质量分数、标准差和中位数。你可以根据自己的需要和数据格式进行修改和扩展。