利用python求一段DNA序列的互補序列

代碼如下：

complement = {'A':'T','G':'C','C':'G','T':'A'}
rev_seq = ''
with open(r'D:\Rosalind\haha.txt','w+') as f1:
    with open(r'D:\Rosalind\rosalind_revc.txt','r') as f:
        dna_seq = f.read()
        dna_seq = list(dna_seq.strip())
        for i in dna_seq:
            rev_seq += complement[i]
            
    rev_seq = rev_seq[::-1]
            
            
    print (rev_seq,file = f1)

這是將從一個txt文件中導入序列，然后將互補后的結果輸出到另外一個文件中。

如果一個段序列不長，直接中python交互式界面完成感覺更方便

先定義的一個字典：    complement = {'C': 'G', 'G': 'C', 'T': 'A', 'A': 'T'}  

然后 for i in list（seq）：

　　　　rev_dna += complement[i]

　　　　rev_dna = rev_dna[::-1]

　　print （rev_dna）

直接得到互補后的結果。