我想找到相当于2个不同核苷酸的IUPAC。
示例:我有A和C,我想要M。 或者:我有R和T,我想要D。
在Biopython中是否有这样做的方法? (这听起来很简单,但我只找到了使用对齐方式来做到这一点的方法,这对我而言不合适。)
谢谢!
答案 0 :(得分:5)
我认为您正在寻找Bio.Data.IUPACData
中的字典,例如:
>>> from Bio.Data import IUPACData
>>> d = {v:k for k,v in IUPACData.ambiguous_dna_values.items()}
>>> d['AC']
'M'
请注意,此处反转字典映射意味着我们丢失了一个密钥,因为'X'
和'N'
都映射到'GATC'
。