def parse_read_id(read_id):
    """
    Parse a read_id into a sequence, strand, position tuple.

    :param read_id: A read id (e.g. "NC_000913.3_-_3061809")
    :return: (sequence, strand, position) tuple e.g. ('NC_000913.3','-','3061809')
    """
    parts = read_id.split(":")
    parts1 = parts[-1].split("-")
    return (":".join(parts[:-1]), "+", int(parts1[0]))