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


用 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

 


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM