Example #1
0
def test_junction_overhang_3():
    """
    Test multiple junction overhangs with softclipped read
    :return:
    """
    cigar = '1S53M56N72M3261N8M'
    min_overhang = 10
    left, right = fr.get_junction_overhangs(cigar, min_overhang)
    assert left == 53 and right == 8
Example #2
0
def test_junction_overhang_1():
    """
    test large junction overhangs on each side
    :return:
    """
    cigar = '1S51M47N90M'
    min_overhang = 10
    left, right = fr.get_junction_overhangs(cigar, min_overhang)
    assert left == 51 and right == 90
Example #3
0
def test_junction_overhang_2():
    """
    Test no junction overhangs
    :return:
    """
    cigar = '255M'
    min_overhang = 10
    left, right = fr.get_junction_overhangs(cigar, min_overhang)
    assert left == -1 and right == -1
Example #4
0
def test_junction_overhang_5():
    """
    Test junction overhang with hardclipped reads
    :return:
    """
    cigar = '66H35M100N100M'
    min_overhang = 10
    left, right = fr.get_junction_overhangs(cigar, min_overhang)
    assert left == 35 and right == 100