-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathREVC-Complementing a Strand of DNA.py
18 lines (15 loc) · 1.16 KB
/
REVC-Complementing a Strand of DNA.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
'''
Rosalind Problems: [REVC] Complementing a Strand of DNA
'''
seq = 'TCAGAAGTATGTGGGGCCCCCGTCCGGCCTGGAATTTAGCGCTATGAGTAACGTGGCGTGGATCGCCATGCATGTAGGTCACAGTTCACACTTAGCCTCATGCCATTCTCTCGGCAAATCGCCGGTGCTTTAACATGATCGGAATCGGTGCGGAGCTCAAACTCAGCTAGGCAACACCTGTGTGTGCAGGGGGGCCACTGATAAAGTGCGGAGACCTCAGTACACTTTGAGGGCGCTTACCCGCTGCCTGATACGCAAGTGTCGTGAACCAGTCTTGGCCGGCCGGGTTACGTGCATCAACTCAGGTGACTGTACGTTTACTTCGTGGCAGAACGAATATTTTCTGGGTTAAAGTTCCGGAGCCTGTGACTCCAGTGTTCACACAGGGGACTTGACTACATCCCTTGATCAGCCAGAGTGGTCTGTTGTGCGCTCGGCATTAAGAATTCTGCTATGCTGCGATATGCAATATGATCATGGGCGCAGTTTCTTTAGGAGAGTCAATTCAGCACTACGGGTGCAGACTTCTCCGCAGCGCGCCTGCGGGGACAGGCACCCGTTAGGCGGTCGATATTTCATGATGGCGGGTCACAACCGGTAGACCTGACGTCTATATATGAGAGGCTGAAGAATTCAAACCTATGATGATTAAGCAAGACCCTTGCCTGTGAAATAGGTTTAGATGCAGTGGAGACCCCTGCGGTCAGCCCTATCACATCTCGAGTTTGTAGGTGGAGGGAAGGAGGCGAACACCATTCACAGCATCGCTCGATGGAAGGACAGGAACAGCTCCTGTCGTGCTCCCCCTTGATATCCCACCTACGTCCCCGAGGGTAACTTTGTCACCG'
seq_c = []
complement = {'C':'G', 'G':'C', 'A':'T', 'T':'A'}
for nt in seq:
seq_c.append(complement[nt])
#notice that the sequence needs to be reversed
seq_c.reverse()
seq_cr = "".join(seq_c)
out = ''
for i in seq_cr:
out = out + i
print(out)