Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,374	G→T	A680S (GCG→TCG)	thrA →	Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,374	0	G	T	100.0%	151.9 / NA	39	A680S (GCG→TCG)	thrA	Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/15);  total (24/15)

AGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATC  >  NC_000913/2228‑2518                                                                                                                                                   |                                                                                                                                                 aGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTcgcg                                                                                                                                              >  1:253775/1‑151 (MQ=255)      tcCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAAc                                                                                                                                         >  2:353736/1‑151 (MQ=255)       ccGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTaaa                                                                                                                                          <  2:408451/149‑1 (MQ=255)                tttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATtct                                                                                                                               >  1:145082/1‑151 (MQ=255)                  tctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTctcg                                                                                                                             >  2:504573/1‑151 (MQ=255)                   ctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTctcgc                                                                                                                            >  1:142515/1‑151 (MQ=255)                    tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTctcgct                                                                                                                           >  2:154483/1‑151 (MQ=255)                         cGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTg                                                                                                                       >  1:564830/1‑150 (MQ=255)                          gCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGt                                                                                                                        >  2:549590/1‑148 (MQ=255)                             ttCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACg                                                                                                                  >  2:189276/1‑151 (MQ=255)                                 tataTCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGAc                                                                                                              >  1:285180/1‑151 (MQ=255)                                  atatCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGAc                                                                                                              <  1:134223/150‑1 (MQ=255)                                    atCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTg                                                                                                           >  1:438475/1‑151 (MQ=255)                                          ggCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTgg                                                                                                     <  1:328308/151‑1 (MQ=255)                                             aaGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGc                                                                                                  >  1:217989/1‑151 (MQ=255)                                              aGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCt                                                                                                 >  2:263861/1‑151 (MQ=255)                                                  aGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCg                                                                                             <  2:357811/151‑1 (MQ=255)                                                      gAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGat                                                                                          >  1:500257/1‑150 (MQ=255)                                                       aaGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGatat                                                                                        >  1:547143/1‑151 (MQ=255)                                                            aTGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattgaa                                                                                    >  1:507450/1‑150 (MQ=255)                                                                      ccGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCtgt                                                                         <  2:569775/151‑1 (MQ=255)                                                                      ccGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCtg                                                                          <  1:120123/150‑1 (MQ=255)                                                                             gACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGccc                                                                  >  2:139478/1‑151 (MQ=255)                                                                                   gCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCag                                                              >  2:347658/1‑149 (MQ=255)                                                                                    cTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCagag                                                            <  1:246289/150‑1 (MQ=255)                                                                                              aaaTGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGa                                                 <  1:167682/151‑1 (MQ=255)                                                                                               aaTGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAg                                                <  2:438563/151‑1 (MQ=255)                                                                                                     ttATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGAt                                          >  1:197271/1‑151 (MQ=255)                                                                                                               ccGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCtt                                >  2:26026/1‑151 (MQ=255)                                                                                                                cGGACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCttt                               <  1:333598/151‑1 (MQ=255)                                                                                                                 ggACCCGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCtttt                              <  2:410148/151‑1 (MQ=255)                                                                                                                     ccGCGAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATgg                          <  2:300461/151‑1 (MQ=255)                                                                                                                      cgcgAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGGGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGc                         >  1:51114/1‑151 (MQ=255)                                                                                                                        cgAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCg                        >  1:156653/1‑150 (MQ=255)                                                                                                                        cgAGATGATCTTTCTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGa                       <  1:66586/151‑1 (MQ=255)                                                                                                                                     cTGGTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTcaca             <  1:283630/148‑1 (MQ=255)                                                                                                                                       ggTATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACt          <  2:6897/149‑1 (MQ=255)                                                                                                                                         tATGGATGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTcgac      >  2:144653/1‑151 (MQ=255)                                                                                                                                              aTGTGTCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGAtc  >  1:336751/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATC  >  NC_000913/2228‑2518

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