Algorithms-In-Bioinformatics-CCA-3/test_dna_data_structures.py at main · siyaagarwal2005/Algorithms-In-Bioinformatics-CCA-3 · 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
from dna_data_structures import DNA, InvalidNucleotideError

""""
This file tests the DNA class functionality.
We will run three test cases-
1. where DNA sequence is valid
2. where DNA sequence is not valid
3. where DNA sequence is empty
"""

print ("Test Case 1- Valid DNA sequence:")
try:
    dna_sequence_1= DNA("ATTCggaCT")
    print(dna_sequence_1)
    print(f"Length: {dna_sequence_1.get_length()}")
    print(f"Counts: {dna_sequence_1.get_nucleotide_count()}")
    print(f"GC-Content: {dna_sequence_1.calculate_GC_content()}%")
    print(f"AT-Content: {dna_sequence_1.calculate_AT_content()}%\n")
except InvalidNucleotideError as e:
    print(f"Error: {e}\n")

print("Test Case 2- Invalid DNA sequence:")
try:
    dna_sequence_2= DNA("AttxyGBAC")
    print(dna_sequence_2)
except InvalidNucleotideError as e:
    print(f"Error: {e}\n")

print("Test Case 3- Empty DNA sequence:")
try:
    dna_sequence_3= DNA()
    print(dna_sequence_3)
    print(f"Length: {dna_sequence_3.get_length()}")
    print(f"GC-Content: {dna_sequence_3.calculate_GC_content()}%")
except InvalidNucleotideError as e:
    print(f"Error: {e}\n")