下記のコードでfastaファイルの生成を試みたのですが、ファイルが生成されないのですが。
3つのfastaフィルを作成しようとしています。
何か間違ってる箇所ありますか?
実行コマンドでは、何も特にエラーもでないし、コメントもでません。
また、もしfastaファイルが新しく生成されたなら、どこにそのファイルが出現するのでしょうか?
よろしくお願いいたします。
補足:初心者で、どこに生成するとか指定がわかりません。
インデントは、入力段階で入力してますが、質問欄で反映されてると良いのですが。
python S = SeqIO.read("vv.fasta", "fasta") S_S = S[20000:22000] def yield_S_t(): for var in S_S: var.seq = S_S.seq.translate() var.id = S_S.id var.description = S_S.description yield var S_t = yield_S_t() S_tt = S_t[10:1000] def generation(): with open('bb.fasta') as file: for list_record in list(SeqIO.parse(file, "fasta")): list_record.seq = list_record.seq[20000:22000] list_record.seq.translate() list_record.id = list_record.id yield list_record g_t = generation() g_tt = g_t[10:1000] seq1 = S_tt seq2 = g_tt All_Seq = [seq1, seq2] SeqIO.write(seq1,"EE.fasta", "fasta") SeqIO.write(seq2, "FF.fasta", "fasta") SeqIO.write(All_Seq,"GG.fasta","fasta") EE.fasta.close() FF.fasta.close() GG.fasta.close() ``` コード S = SeqIO.read("vv.fasta", "fasta") S_S = S[20000:22000] def yield_S_t(): for var in S_S: var.seq = S_S.seq.translate() var.id = S_S.id var.description = S_S.description yield var S_t = yield_S_t() S_tt = S_t[10:1000] def generation(): with open('bb.fasta') as file: for list_record in list(SeqIO.parse(file, "fasta")): list_record.seq = list_record.seq[20000:22000] list_record.seq.translate() list_record.id = list_record.id yield list_record g_t = generation() g_tt = g_t[10:1000] seq1 = S_tt seq2 = g_tt All_Seq = [seq1, seq2] SeqIO.write(seq1,"EE.fasta", "fasta") SeqIO.write(seq2, "FF.fasta", "fasta") SeqIO.write(All_Seq,"GG.fasta","fasta") EE.fasta.close() FF.fasta.close() GG.fasta.close()