使用Biopython读取fastq文件

使用Biopython读取fastq文件

读取fastq文件，输出碱基序列和预测的准确度

from Bio import SeqIO

with open("./data2/ERR000020_2.fastq") as handle:
    record = SeqIO.parse(handle,"fastq")
    for lin in record:
        print("line information \n")
        print(lin.id)
        print(lin.seq)
        print(len(lin.letter_annotations['phred_quality']))
        q = lin.letter_annotations['phred_quality']
        num = [1- 10**(qi/(-10.0)) for qi in q]
        print(num)

这里是具体代码的结果