如何得到给定序列的互补序列以及反向互补序列

用 python 实现如下:

 1 #!/usr/bin/python
 2 # Complementing DNA
 3
 4 my_dna = "ACTGATCGATTACGTATAGTATTTGCTATCATACATATATATCGATGCGTTCAT"
 5 # 由于python 区分大小写，所以先替换成小写可以有效避免后面重复替换的问题。
 6 complement_AT = my_dna.replace("A","t")
 7 complement_TA = complement_AT.replace("T","a")
 8 complement_GC = complement_TA.replace("G","c")
 9 complement_CG = complement_GC.replace("C","g")
10 # 把所以小写变为大写即可。
11 complement_DNA = complement_CG.upper()
12 reverse_complement_DNA = complement_DNA[::-1]  # str[::-1] used to reverse the string
13
14 print(complement_DNA)
15 print(reverse_complement_DNA)

运行结果：

ubuntu$ python3 complement_py3.py
TGACTAGCTAATGCATATCATAAACGATAGTATGTATATATAGCTACGCAAGTA
ATGAACGCATCGATATATATGTATGATAGCAAATACTATACGTAATCGATCAGT