1

FastqGeneralIterator を使用していますが、fastq ファイルの 1 行目から @ が削除され、3 行目の情報も削除されることがわかりました (3 行目全体が削除されます)。次の方法で 1 行目に @ を追加しました。

for line in open("prova_FiltraN_CE_filt.fastq"):
fout.write(line.replace('SEQ', '@SEQ'))

+ で始まり、その後に何もない 3 行目も追加したいと思います。例えば：

    @SEQILMN0
    TCATCGTA....
    +
    #<BBBFFF.....

誰かが私を助けることができますか？

4

1 に答える 1

2

使用できます、文字列の書式設定操作 %

from Bio.SeqIO.QualityIO import FastqGeneralIterator
with open("prova_FiltraN_CE_filt.fastq", "rU") as handle:
    for (title, sequence, quality) in FastqGeneralIterator(handle):
        print("@%s\n%s\n+\n%s" % (title, sequence, quality))

fastqを使用して、印刷形式を取得しますFastqGeneralIterator

@SEQILMN0
TCATCGTA....
+
#<BBBFFF....
于 2015-10-13T18:09:58.800 に答える