Hi folks, When trying to use pepinfo on sequences containing peptide ambiguity codes (e.g. BJXZ) the program gives an error. For example running pepinfo on UniProtKB:Q9DY04 gives:
Error: At position 83 in seq, couldn't find key X This is a reference to these amino acid codes not having corresponding entries in the data files used by pepinfo: - Eaa_properties.dat - Eaa_hydropathy.dat It would be good if this was handled more gracefully, either by providing "average" values in the data files, or by handling these cases with a warning and an appropriate "null" value. FWIW the same issue comes when handling sequences with unusual amino acids (i.e. 'O' and 'U'), although there should be actual figures for these. Does anyone have already fixed data files that could be used instead? All the best, Hamish -- ============================================================ Mr Hamish McWilliam, Web Production, European Bioinformatics Institute (EMBL-EBI), European Molecular Biology Laboratory, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD United Kingdom URL: http://www.ebi.ac.uk/ ============================================================ _______________________________________________ EMBOSS mailing list EMBOSS@lists.open-bio.org http://lists.open-bio.org/mailman/listinfo/emboss