def remove_outlier_lengths(seq_record_data, protein_analysis_data, max_length=2000):
    """
    Parameters:
    -----------
    seq_record_data : `list` of `Bio.SeqRecord.SeqRecord`
        List of sequence records.    
    protein_analysis_data : `list` of `Bio.SeqUtils.ProtParam.ProteinAnalysis`
        List of sequence records.
    max_length : `int`
        The maximum length of the sequence.        
    
    Returns:
    --------
    modified_data : `list` of `Bio.SeqRecord.SeqRecord`
        The data with outlier lengths removed.
    """    
    modified_data = [seq_record_data[i] for i, example in enumerate(protein_analysis_data) if example.length < max_length]
    print("Data removed: ", len(seq_record_data)-len(modified_data))
    return modified_data