Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,129,362	G→T	V231L (GTG→TTG)	metB →	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,129,362	0	G	T	100.0%	163.8 / NA	45	V231L (GTG→TTG)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/15);  total (30/15)

ACGTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCGTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAG  >  NC_000913/4129285‑4129439                                                                              |                                                                              aCGTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGa                                                                            >  1:837022/1‑81 (MQ=255)  cGTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGAc                                                                           <  1:836947/81‑1 (MQ=255)  cGTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGAc                                                                           <  2:836948/81‑1 (MQ=255)   gTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACg                                                                          >  1:836997/1‑81 (MQ=255)      gTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGggc                                                                       >  1:836991/1‑81 (MQ=255)        ggCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGggcgg                                                                     >  2:837021/1‑81 (MQ=255)        ggCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGggcgg                                                                     >  1:837036/1‑81 (MQ=255)            ggCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTTGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGcgcg                                                                 >  2:837037/1‑81 (MQ=255)            ggCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGcgcg                                                                 >  1:837011/1‑81 (MQ=255)              cGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGtt                                                               >  1:837049/1‑81 (MQ=255)                  gtgATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGa                                                            >  1:837010/1‑80 (MQ=255)                    gATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAg                                                         >  1:837024/1‑81 (MQ=255)                     aTTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGc                                                        <  2:836949/81‑1 (MQ=255)                       tGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTa                                                      >  2:837006/1‑81 (MQ=255)                         cTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACTATATTGGCTTGACGGGCAGCGCGTTTGACAGCTAt                                                     <  2:836951/80‑1 (MQ=255)                         cTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTAt                                                     <  2:836950/80‑1 (MQ=255)                           aaaGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATctg                                                  >  2:837013/1‑81 (MQ=255)                              gACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATctgctg                                               <  2:836952/81‑1 (MQ=255)                              gACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATctgctg                                               >  1:837044/1‑81 (MQ=255)                               aCCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCtgctgc                                              >  1:837026/1‑81 (MQ=255)                               aCCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCtgctgc                                              >  1:837009/1‑81 (MQ=255)                               aCCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCtgctgc                                              >  1:837002/1‑81 (MQ=255)                                 ccGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTa                                            <  1:836953/81‑1 (MQ=255)                                  cGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTAc                                           <  1:836955/81‑1 (MQ=255)                                  cGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTAc                                           >  2:837018/1‑81 (MQ=255)                                  cGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTAc                                           <  2:836954/81‑1 (MQ=255)                                  cGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTAc                                           <  2:836956/81‑1 (MQ=255)                                   ggACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACg                                          >  2:837027/1‑81 (MQ=255)                                             aCTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGa                                <  2:836957/81‑1 (MQ=255)                                              cTGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGaa                               >  2:837017/1‑81 (MQ=255)                                               tGAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAAc                              >  1:837062/1‑81 (MQ=255)                                                gAACTGGCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACg                             >  1:837030/1‑81 (MQ=255)                                                      gCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGt                        <  2:836958/80‑1 (MQ=255)                                                      gCCTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTg                       >  2:837028/1‑81 (MQ=255)                                                       ccTGGTGGGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGc                      <  2:836959/81‑1 (MQ=255)                                                          ggtggGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCcgc                   >  2:837032/1‑81 (MQ=255)                                                          ggtggGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCcgc                   >  1:837034/1‑81 (MQ=255)                                                            tggGCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGt                 <  1:836960/81‑1 (MQ=255)                                                             gggCAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTa                >  1:837035/1‑81 (MQ=255)                                                                cAAACAATATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATg              >  1:837043/1‑80 (MQ=255)                                                                     aaTATTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTg        >  1:837047/1‑81 (MQ=255)                                                                       tatTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGgc      <  2:836961/81‑1 (MQ=255)                                                                       tatTGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGgc      >  1:837046/1‑81 (MQ=255)                                                                          tGGCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGgcg     >  1:837031/1‑79 (MQ=255)                                                                           ggCTTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGgcgcag  >  2:837041/1‑81 (MQ=255)                                                                              |                                                                              ACGTAGTGGCCGGCGTGGTGATTGCTAAAGACCCGGACGTTGTCACTGAACTGGCCTGGTGGGCAAACAATATTGGCGTGACGGGCGGCGCGTTTGACAGCTATCTGCTGCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAG  >  NC_000913/4129285‑4129439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑