LeetCode-PyBites-PythonMorsels-OtherProgrammingChallenges/bite_262.py at main · jv199768/LeetCode-PyBites-PythonMorsels-OtherProgrammingChallenges · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90

# Online Python - IDE, Editor, Compiler, Interpreter

from collections import Counter

DNA_SEQUENCES = [
    (
        (  # from https://www.ncbi.nlm.nih.gov/nuccore/NC_000913.3
            "tagtgaaagatattcatttcgaaggccttcagcgtgtcgccgttggtgcggccctcctca"
            "gtatgccggtgcgcacaggcgacacggttaatgatgaagatatcagtaataccattcgcg"
            "ctctgtttgctaccggcaactttgaggatgttcgcgtccttcgtgatggtgatacccttc"
            "tggttcaggtaaaagaacgtccgaccattgccagcattactttctccggtaacaaatcgg"
            "tgaaagatgacatgctgaagcaaaacctcgaggcttctggtgtgcgtgtgggcgaatccc"
            "tcgatcgcaccaccattgccgatatcgagaaaggtctggaagacttctactacagcgtcg"
            "gtaaatatagcgccagcgtaaaagctgtcgtgaccccgctgccgcgcaaccgtgttgacc"
            "taaaactggtgttccaggaaggtgtgtcagctgaaatccagcaaattaacattgttggta"
            "accatgctttcaccaccgacgaactgatctctcatttccaactgcgtgacgaagtgccgt"
            "ggtggaacgtggtaggcgatcgtaaataccagaaacagaaactggcgggcgaccttgaaa"
            "ccctgcgcagctactatctggatcgcggttatgcccgtttcaacatcgactctacccagg"
            "tcagtctgacgccagataaaaaaggtatttacgtcacggtgaacatcaccgaaggcgatc"
            "agtacaagctttctggcgttgaagtgagcggcaaccttgccgggcactccgctgaaattg"
            "agcagctgactaagatcgagccgggtgagctgtataacggcaccaaagtgaccaagatgg"
            "aagatgacatcaaaaagcttctcggtcgctatggttatgcctatccgcgcgtacagtcga"
            "tgcccgaaattaacgatgccgacaaaaccgttaaattacgtgtgaacgttgatgcgggta"
            "accgtttctacgtgcgtaagatccgttttgaaggtaacgatacctcgaaagatgccgtcc"
            "tgcgtcgcgaaatgcgtcagatggaaggtgcatggctggggagcgatctggtcgatcagg"
            "gtaaggagcgtctgaatcgtctgggcttctttgaaactgtcgataccgatacccaacgtg"
            "ttccgggtagcccggaccaggttgatgtcgtctacaaggtaaaagagcgcaacaccggta"
            "gcttcaactttggtattggttacggtactgaaagtggcgtgagcttccaggctggtgtgc"),
        51.35,
    ),
    (
        (  # from https://www.ncbi.nlm.nih.gov/nuccore/CP010052.1
            "gttaatatttatgattcctgaaaagaaatcaatcgcaatcatgaaagaactaagcattgg"
            "aaatacaaagcaaatgctgatgattaatggagttgacgtgaaaaatccattgctgctttt"
            "tttacatggcgggccgggaacgccgcaaatcggatatgttagacattatcaaaaagagct"
            "ggaacagtattttacagtagttcattgggatcagagaggatcggggctttcttattctaa"
            "gcgaatttcgcatcactctatgacaataaatcacttcattaaagatacaatccaagtcac"
            "tcaatggcttttagctcatttttcaaaatcaaaactttacctagccggtcattcttgggg"
            "atcaatactggcgcttcatgtgctgcagcagcgtcctgatttattttacacgtattatgg"
            "aatcagccaggttgttaacccgcaagatgaagaatcaactgcttatcaacatattcgtga"
            "aatttccgaatcaaaaaaagccagcatattatctttccttacacgtttcattggtgctcc"
            "gccttggaagcaggatatccagcaccttatctatcggttttgtgttgagctaaccagggg"
            "aggattcactcaccgtcatcgtcaatctctcgctgtattatttcaaatgcttactggcaa"
            "tgagtatggagtgcggaacatgcacagcttccttaatggattgcgcttcagtaaaaaaca"
            "tttaactgatgagttgtaccggtttaatgcttttacatcagttccttctattaaagtacc"
            "gtgtgttttcatttcagggaaacatgacttaattgttcctgcagaaatatcgaaacagta"
            "ttatcaagaacttgaggcacctgaaaagcgctggtttcaatttgagaattcagctcacac"
            "cccgcatattgaggagccatcattattcgcgaacacattaagtcggcatgcacgacacca"
            "tttatgatagatccttgataaataagaaaaacccctgtataataaaaaaagtgtgcaaat"
            "catgcatattttaaataagtcttgcaacatgcgcctattttctgtataatggtgtatgtt"
            "ggtctttgactgcgatgaagtgagaggttgctgacacacccggccgctttgccatggcaa"
            "ggtgttcaggtttttctcacggagaactgtctaacgtgatgtaggcgaaaaggagggaaa"
            "ataatggcaaaacaaaaaattcgtattcgtttgaaagcatatgatcatagaatccttgat"),
        39.37,
    ),
    (
        (  # from https://www.ncbi.nlm.nih.gov/nuccore/CP010052.1
            # but all As removed, upper and lower case
            "gtttttttgttcctggtctcgctctggctgcttgg"
            "tcgctgctgtgtttgggttgcgtgtccttgctgctttt"
            "tttctggcgggccgggcgccgctcggttgttgctttcggct"
            "ggcgtttttcgtgttcttgggtcggggtcggggctttcttttct"
            "gcgtttcgctcctcttgcttccttcttgtctccgtcc"
            "tctggcttttgctctttttctcctttcctgccggtcttcttgggg"
            "tctctggcgcttctgtgctgcgcgcgtcctgtttttttccgttttgg"
            "tcgccggttgttcccgcgtGGtcctgctttcctttcgtg"
            "tttccgtcgccgcttttctttccttccgtttcttggtgctcc"
            "gccttgggcggttccgccctttcttcggttttgtgttggctccgggg"
            "ggttcctcccgtctcgtctctctcgctgttttttctgcttctggc"
            "tggttgggtgcggctgccgcttcctttggttgcgcttcgtc"
            "tttctgtggttgtccggttttgcttttctcgttccttctttgtcc"
            "gtgtgttttctttcgggCtgcttttgttcctgcgttcgcgt"),
        55.88,
    ),
]

def calculate_gc_content(sequence):
    """
    Receives a DNA sequence (A, G, C, or T)
    Returns the percentage of GC content (rounded to the last two digits)
    """
    c = Counter(sequence.upper())
    res = (c["G"] + c["C"]) / (c["A"] + c["G"] + c["C"] + c["T"]) * 100
    return round(res, 2)

if __name__ == "__main__":
    for dna, expected in DNA_SEQUENCES:
        result = calculate_gc_content(dna)
        print(result == expected)