def test_map_sequence_subsequence(): source_seq = 'ACTGATGC' dest_seq = 'TGAT' expected_positions = [3, 4, 5, 6] positions = map_sequence(source_seq, dest_seq) assert expected_positions == positions
def test_map_sequence_insertion(): source_seq = 'ACTGATGC' dest_seq = 'ACTGCATGC' expected_positions = [1, 2, 3, 4, None, 5, 6, 7, 8] positions = map_sequence(source_seq, dest_seq) assert expected_positions == positions
def test_map_sequence_low_quality(): source_seq = 'AGAGCGAACCGATTC' dest_seq = 'NNNNNNNNNNNATTC' expected_positions = 11 * [None] + [12, 13, 14, 15] positions = map_sequence(source_seq, dest_seq) assert expected_positions == positions
def test_map_sequence_identical(): source_seq = 'ACTG' dest_seq = 'ACTG' expected_positions = [1, 2, 3, 4] positions = map_sequence(source_seq, dest_seq) assert expected_positions == positions