想通过python2.7 写一段程序,实现 DNA反向互补,如下列字符串:由文本文档1.txt输入(由文档输入是关键)
ccccgtatcgccaagggcaaggtcatcatcggcaacggcggtgcagagttcggcgcccgtggctttgtgacggcgtatgacgctgaaacgggcaagatggactggcgcttcttcaccgttccgaaccctgacaacaagccggacggcgcagcgtctgacgacgtgctgatgtccaaggcttatccgacatggggcaagggcggcgcatggaagcagcagggcggtggcggaaccgtctgggattcgctgatctatgaccccgtaacggatcttgtttaccttggcgtcggcaacggctcaccctggaactacaagttccgctcggaaggaaaa
ntComplement = {'A':'t', 'C':'g', 'T':'a', 'G':'c', 'a':'t', 'c':'g', 't':'a', 'g':'c'}
反向互补后的结果输出到文档 1rv.txt 如下
ttttccttccgagcggaacttgtagttccagggtgagccgttgccgacgccaaggtaaacaagatccgttacggggtcatagatcagcgaatcccagacggttccgccaccgccctgctgcttccatgcgccgcccttgccccatgtcggataagccttggacatcagcacgtcgtcagacgctgcgccgtccggcttgttgtcagggttcggaacggtgaagaagcgccagtccatcttgcccgtttcagcgtcatacgccgtcacaaagccacgggcgccgaactctgcaccgccgttgccgatgatgaccttgcccttggcgatacgggg
希望有高手写一段完整的代码,能运行的,一定采纳
不要从网上copy来不能运行的,多谢啦
https://ask.csdn.net/questions/688856
程序已经给你
程序非网上粘贴,而且调试通过
有问题请追问。
实在不好意思,我不是学生物的。之前的程序有错误。
import sys
def DNA_complement(sequence):
sequence = sequence.lower()
sequence = sequence.replace('a', 'T')
sequence = sequence.replace('t', 'A')
sequence = sequence.replace('c', 'G')
sequence = sequence.replace('g', 'C')
return sequence.lower()
def DNA_reverse(sequence):
sequence = sequence.lower()
return sequence[::-1]
file_name = '1.txt'
f = open(file_name)
f1 = open(r"1rv.txt", "a+")
line = f.readline()
s = ""
while line:
s = DNA_reverse(DNA_complement(line))
print (s)
f1.write(s)
line = f.readline()
f.close()
f1.close()
ttttccttccgagcggaacttgtagttccagggtgagccgttgccgacgccaaggtaaacaagatccgttacggggtcatagatcagcgaatcccagacggttccgccaccgccctgctgcttccatgcgccgcccttgccccatgtcggataagccttggacatcagcacgtcgtcagacgctgcgccgtccggcttgttgtcagggttcggaacggtgaagaagcgccagtccatcttgcccgtttcagcgtcatacgccgtcacaaagccacgggcgccgaactctgcaccgccgttgccgatgatgaccttgcccttggcgatacgggg
·如果还有问题,请在这里追问,不要再发新问题了。抱歉·
您好,结尾处 aaggaaaa 反向互补后应该变成结果序列的开头部分,因该是 ttttcctt ,不是AAAACCAA,好像 a 都没有反补,
您看是不是,完成后这40个币您也收了吧
这回好用了,多谢您了