def load_plot_regions(fn_regions):
    if fn_regions == None: return []
    plot_regions = []
    for line in open(fn_regions,'r').readlines():
        if line[0] == "#": continue
        print line
        sline = line.split()
        uID = "%s:"%(sline[1])
        uID += ":".join(sline[2:5])
        plot_regions.append(uID)
        print uID

    return plot_regions