From f56c0607c311f07158c04ea1baf23e9904a4742a Mon Sep 17 00:00:00 2001 From: Minh Bui Date: Thu, 7 Jun 2018 17:18:51 +1000 Subject: [PATCH] support amino-acid O as the 22nd amino-acid (reported by Cuong) --- alignment/alignment.cpp | 2 ++ 1 file changed, 2 insertions(+) diff --git a/alignment/alignment.cpp b/alignment/alignment.cpp index 093e6807..99a4af3c 100644 --- a/alignment/alignment.cpp +++ b/alignment/alignment.cpp @@ -1113,6 +1113,7 @@ void Alignment::buildStateMap(char *map, SeqType seq_type) { map[(unsigned char)'J'] = 22; // I or L map[(unsigned char)'*'] = STATE_UNKNOWN; // stop codon map[(unsigned char)'U'] = STATE_UNKNOWN; // 21st amino acid + map[(unsigned char)'O'] = STATE_UNKNOWN; // 22nd amino acid return; case SEQ_MULTISTATE: @@ -1201,6 +1202,7 @@ char Alignment::convertState(char state, SeqType seq_type) { if (state == 'J') return 22; if (state == '*') return STATE_UNKNOWN; // stop codon if (state == 'U') return STATE_UNKNOWN; // 21st amino-acid + if (state == 'O') return STATE_UNKNOWN; // 22nd amino-acid loc = strchr(symbols_protein, state); if (!loc) return STATE_INVALID; // unrecognize character