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
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
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
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
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