def check_find_known_in_range(query_junc):
txid = query_junc.chrom
minus_range,plus_range = find_match_range(query_junc,seqs,20)
found_known = find_known_in_range(query_junc,minus_range,plus_range,all_known_juncs)
expected = set([str(X) for X in known_juncs.get(txid,[])])
found_known = set([str(X) for X in found_known])
assert_set_equal(expected,found_known,("Known junction not found for %s. "+\
"Query: %s\n Expected %s\n " +\
"Got %s\n match range: (%s,%s)") % (txid,
str(query_junc),
expected,
found_known,
minus_range,
plus_range))
def check_find_known_in_range(query_junc):
txid = query_junc.chrom
minus_range, plus_range = find_match_range(query_junc, seqs, 20)
found_known = find_known_in_range(query_junc, minus_range, plus_range,
all_known_juncs)
expected = set([str(X) for X in known_juncs.get(txid, [])])
found_known = set([str(X) for X in found_known])
assert_set_equal(expected,found_known,("Known junction not found for %s. "+\
"Query: %s\n Expected %s\n " +\
"Got %s\n match range: (%s,%s)") % (txid,
str(query_junc),
expected,
found_known,
minus_range,
plus_range))
def check_find_none_when_none_known_in_range(query_junc):
minus_range, plus_range = find_match_range(query_junc, seqs, 20)
found_known = find_known_in_range(query_junc, minus_range, plus_range,
all_known_juncs)
found_known = set([str(X) for X in found_known])
assert_set_equal(found_known, set([]))
def check_find_none_when_none_known_in_range(query_junc):
minus_range, plus_range = find_match_range(query_junc,seqs,20)
found_known = find_known_in_range(query_junc,minus_range,plus_range,all_known_juncs)
found_known = set([str(X) for X in found_known])
assert_set_equal(found_known,set([]))