def normalized_gene_length(gene_length, read_length):
    """
    Compute the normalized gene length with the average read length if set
    Report in reads per kilobase
    """
    if read_length < 1:
        read_length = 1
    return (abs(gene_length - read_length) + 1) / 1000.0
