Tolong beri petunjuk Anda - Saya telah menghabiskan banyak waktu untuk kode ini tetapi belum dapat menyelesaikannya. Saya memiliki kode yang memberikan output dalam format yang saya inginkan namun hasil yang ditampilkan tidak 100% ...

1
Apex 9 April 2021, 17:55

1 menjawab

Jawaban Terbaik

Sepertinya ada bug dalam kode Anda. Ubah fungsi same_seq sebagai berikut:

def same_seq(a_record, brecord):
    window_size=6
    for j in range(len(a_record.seq)- window_size):
        for k in (a_record.seq.reverse_complement()[j: j + 6].split()):
            if brecord.seq.find(k) != -1:
                    return 1
    return 0

Semoga berhasil!

1
hamid ghasemi 9 April 2021, 19:41