在數(shù)據(jù)分析中，數(shù)據(jù)以臟形式傳送給我們是很常見的，其中存在與轉(zhuǎn)錄或下載方式相關(guān)的錯(cuò)誤。由于我們知道DNA的任何序列都必須由四個(gè)堿基組成'a'，'g'因此DNA中出現(xiàn)'t'的'c'任何其他字母一定是錯(cuò)誤的。編寫一個(gè)函數(shù)clean(dna)，返回一個(gè)新的 DNA 字符串，其中除 A、C、G 或 T 之外的每個(gè)字符都被替換為 N。例如，clean('goat')應(yīng)返回字符串'gnat'。您可以假設(shè) dna 全部為小寫，但不要假設(shè)任何有關(guān)錯(cuò)誤字符的性質(zhì)（例如，它們甚至可能被意外轉(zhuǎn)錄為數(shù)字）。clean('') → ''clean('agct7ttczttctgactgcaacgggcaatatgtctctxtgtggattaaaaaaagagtgtcygatagcagcttctgaactggttacctgcc') → 'agctnttcnttctgactgcaacgggcaatatgtctctntgtggattaaaaaaagagtgtcngatagcagcttctgaactggttacctgcc'clean('gtgagtaaattaaaattttnttgacttaggtcactaaptactttaaccaatataggbatagcgcacagacagataaaaattacagagtac') → 'gtgagtaaattaaaattttnttgacttaggtcactaantactttaaccaatataggnatagcgcacagacagataaaaattacagagtac'使用for循環(huán)
2 回答
HUH函數(shù)
TA貢獻(xiàn)1836條經(jīng)驗(yàn) 獲得超4個(gè)贊
我希望我沒有替你做功課。我看到已經(jīng)使用 .sub 發(fā)布了答案，但您只要求使用循環(huán)
def clean(text):
cleaned_text=""
for i in range(0, len(text)):
if text[i] in "agtc":
cleaned_text=cleaned_text+text[i]
else:
cleaned_text=cleaned_text+"n"
return cleaned_text
print(clean("agct7ttczttctgactgcaacgggcaatatgtctctxtgtggattaaaaaaagagtgtcygatagcagcttctgaactggttacctgcc"))
# returns agctnttcnttctgactgcaacgggcaatatgtctctntgtggattaaaaaaagagtgtcngatagcagcttctgaactggttacctgcc
HUX布斯
TA貢獻(xiàn)1876條經(jīng)驗(yàn) 獲得超6個(gè)贊
使用re.sub：
Import re
dna = 'gtgagtaaattaaaattttnttgacttaggtcactaaptactttaaccaatataggbatagcgcacagacagataaaaattacagagtac'
dna = re.sub(r'[^ACTG]','N',dna.upper())
print(dna)
添加回答
舉報(bào)
0/150
提交
取消