读取 fastq 文件并计算有多少碱基具有 Phred 碱基读取质量-解网

问：

我需要帮助用 Python3 解决以下问题：
编写一个读取 fastq 文件的 Python 程序，并计算有多少碱基的 Phred 碱基读取质量为零，介于 1 和 10（含）、11 和 20、21 和 30、31 和 40 之间，以及 40 以上。

我从以下方面开始：

def decode(c):
    return ord(c) - 33

letters = "II93882$%@%%@"

values = map(decode, letters)
values = list(values)

print (values)

我正在尝试找到一种方法来使我的最终结果如下所示：

Base quality zero: XXX bases

Base quality 1-10: YYY bases

Base quality 11-20: ZZZ bases

Base quality 21-30: QQQ bases

Base quality 31-40: RRR bases

Base quality above 40: SSS bases

Python FastQ

读取 fastq 文件并计算有多少碱基具有 Phred 碱基读取质量

Reading a fastq file and calculating how many bases have Phred base read quality

评论