def test_find_start():
assert DNA('ATGGG').find_first_start_site == 0
assert DNA('CCCCATG').find_first_start_site == 4
with pytest.raises(TypeError):
DNA('GGG').find_first_start_site
assert DNA('GTC').compliment == DNA('CAG')
assert DNA('ATC').compliment == DNA('TAG')
def test_complimentary_sequence_works():
assert DNA('GTC').complimentary_sequence == DNA('CAG')
assert DNA('ATC').complimentary_sequence == DNA('TAG')
assert DNA('GTC').complimentary_sequence == DNA('CAG')
assert DNA('ATC').complimentary_sequence == DNA('TAG')
def test_transcribe():
assert DNA('GTC').transcribe() == 'GAC'
assert DNA('ATC').transcribe() == 'GAU'
def test_bad_sequence_raises_error():
with pytest.raises(ValueError):
DNA('ATB')
def test_find_start_codons():
"""New test to test the function to find start codons"""
assert DNA('ATGGTACATGCGA').find_start_codons() == [0, 7]
def test_triplets():
assert DNA('AAA').split_DNA_triplets == ['AAA']
assert DNA('AAATTTGGG').split_DNA_triplets == ['AAA', 'TTT', 'GGG']
assert DNA('AAAT').split_DNA_triplets == ['AAA', 'T']
def test_gc_content():
assert DNA('ATTTATGGCC').gc_content == 0.4
assert DNA('AGGTATGGCC').gc_content == 0.6
assert DNA('ATAT').gc_content == 0
def test_is_gc_rich():
assert DNA('GTGT').gc_content() == 0.5
def test_find_start_codons():
assert DNA('ATGGTACATGCGA').find_start_codons() == [0, 7]
def test_find_start():
assert DNA('ATGGG').find_first_start_site == 0
assert DNA('CCCCATG').find_first_start_site == 4
with pytest.raises(TypeError):
DNA('GGG').find_first_start_site()
def get_reverse_complement(seq, reverse=True):
""" Take a DNA sequence and returns its complement"""
complement = str(DNA(seq).compliment())
if reverse:
complement = complement[::-1]
click.echo(complement)