Esempio n. 1
0
def test_non_ag_mismatches_5():
    # One AG(+) mismatch, two non-ag mismatch (TC, TC)
    read_seq = 'CCCCCCCCCCGCCCCC'
    md = '10A3T0T10'
    sense = True
    assert fr.non_ag_mismatches(read_seq, md, sense) == 2
Esempio n. 2
0
def test_non_ag_mismatches_4():
    # One TC(-) mismatch, one non-ag mismatch (TA)
    read_seq = 'ATCGA'
    md = '2T1T'
    sense = False
    assert fr.non_ag_mismatches(read_seq, md, sense) == 1
Esempio n. 3
0
def test_non_ag_mismatches_2():
    # One non-ag mismatch
    read_seq = 'CCAAAATTCAGCCCGCGAAGGCATGACGTCAGCGCAAGGCAGTAGTTTCCAGAAGAACTCTGTCGTCTACCTTAATGCCTCAAATGCGAACCCGCTTCGGCCATCCTTCTCGCTCAGAGAATGGATTAGAGTTCTA'
    md = '35G99'
    sense = False
    assert fr.non_ag_mismatches(read_seq, md, sense) == 1
Esempio n. 4
0
def test_non_ag_mismatches_3():
    # One AG mismatch
    read_seq = 'ATCG'
    md = '2T1'
    sense = False
    assert fr.non_ag_mismatches(read_seq, md, sense) == 0
Esempio n. 5
0
def test_non_ag_mismatches_1():
    # No mismatches
    read_seq = 'ACGCCTCGGAGTCAACCCGGAGGAAGTTTTGGCGGATCTTCGTGCTCGTAATCAATTCCAATAAATATTCTTTGCCCTAAATACTTTAAATTATCCATCTGACAACTAAAATTTCGGTTCTTCTTGGCTTCTTCTATTTGTGAAATGGTTTATTT'
    md = '155'
    sense = True
    assert fr.non_ag_mismatches(read_seq, md, sense) == 0