Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,903	G→A	E523K (GAA→AAA)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,903	0	G	A	100.0%	103.1 / NA	32	E523K (GAA→AAA)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/16);  total (16/16)

GACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAA  >  NC_000913/1680‑2141                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               gACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCGGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCTTTAATCTGGAAAACGGGCAGAAAGAACTTGCGCAAGCCAAAGAGCCGt                                                                                                                                                                                                                     <  1:132784/251‑1 (MQ=255)    cACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGttt                                                                                                                                                                                                                   >  1:172661/1‑250 (MQ=255)    cACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGttt                                                                                                                                                                                                                   >  1:222285/1‑250 (MQ=255)                   tACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTa                                                                                                                                                                                                    <  1:295809/250‑1 (MQ=255)                          caGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTGTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTGAATGCGCCt                                                                                                                                                                                            <  1:346322/251‑1 (MQ=255)                                 tGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                                                                                                                                                                      >  1:49764/1‑250 (MQ=255)                                  gTTCAATACGGATCAGGTTAGCGAAGTGATTGTGATTGGAGTCGGGGGCGTTGGCGGAGCGCTGCTGGAGCAACTGAAGAGTCAGCAAAGCAGGCTGAAGAATAAAAATATCGACTTACGGGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAAATGGCGCAAGACAAAGAGACGGTTAATCGAGGGCGCTTAATTCGCCTCGTGaa                                                                                                                                                                                      >  2:1925/1‑249 (MQ=255)                                  gTTCAATACCGATCATGTTATCGAAGTGTTTGTGATTGGCGTCTGTGGCGTTGGCGGTTCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTTAAGAATAAACATAACGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCGGCTCACCAATGTACATGGCCTTAATCTGGAAAACGGGCAGAAAGAACTGGCACAAGCCAAAGAGCCGTTTGATCTCGGGCGCTTAAATCGCCCCGTGaaa                                                                                                                                                                                     <  2:675/250‑1 (MQ=255)                                               gaGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTGACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACGGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgcg                                                                                                                                                                       <  2:325229/250‑2 (MQ=255)                                                       cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGCGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCgg                                                                                                                                                               >  1:35258/1‑251 (MQ=255)                                                        gacgtgtTTGTGGTTGTCGCCGGTGGCGTTGACGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGGTGAAGAATATACATATCGGCTTACGTGTCTGCGGGGTTGCCAACTCGAAGGATCTGCTCGTCAATGTACATGGCATTAAGTTGGAAAACTGGCAGAAATAACTGGCTCAAGCCAAAGTTCCGTGTAATCTCGGGAGCTTAATTCGCCTCTTGAAAGAAGATCAGCTGCGGTACCCGGt                                                                                                                                                              <  2:166128/248‑1 (MQ=255)                                                                       ggCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCAc                                                                                                                                               >  2:175253/1‑251 (MQ=255)                                                                            cGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcca                                                                                                                                          <  2:327955/251‑1 (MQ=255)                                                                                                                          aGCTGGCCGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTAGAAGGCTCTGCTCACAAATGTACATGGCCTTACTCTGGAAAACTAGCAGAAAGACCTGGCGCAAGCCACAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGCAAGCATATCATCAGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCAAGGCAGTGGAGGATCAATATGACGACTTCCTGCGCGAAGGtat                                                                                            >  1:414707/1‑249 (MQ=255)                                                                                                                                        tAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                                                                               <  1:259877/234‑1 (MQ=255)                                                                                                                                        tAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                                                                               >  2:259877/1‑234 (MQ=255)                                                                                                                                         aaaCATATCGACTTACGTGGCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAAGATCATCTTCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCATGCAGTGGCGGATCAAGATGCCGACGTCCTGCGCGAAGGTTTCCGCGTTGTCACGcc                                                                             <  2:416943/251‑1 (MQ=255)                                                                                                                                              tatCGAATTACGTGACTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATAGCCTTAATATGGCAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCCGACTTCCTGCGCGACGGTTTACACGTTGTCACGCCGAAc                                                                         >  2:157496/1‑250 (MQ=255)                                                                                                                                                         gtgtCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAAATGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCCGCCAGGCAGTGGCGGATCTATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACCCCGAACAAAAAGGCCacc                                                             >  1:422100/1‑249 (MQ=255)                                                                                                                                                           gtCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAaca                                                            <  2:187201/250‑1 (MQ=255)                                                                                                                                                                 ggTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACAACtcgtc                                                     >  1:392125/1‑251 (MQ=255)                                                                                                                                                                  gtgtTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACCCCGAACAAAAAGGCCAACACCccgtcg                                                    >  2:60424/1‑251 (MQ=255)                                                                                                                                                                       gCCAACTCGAAGGCTCTGCTCAACAATGTTCATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGa                                               <  2:93842/251‑1 (MQ=255)                                                                                                                                                                                 aGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCAt                                     <  1:14858/251‑1 (MQ=255)                                                                                                                                                                                    ctctGCTCACCAATGTACATGGCCTTAATCTGGAAACCTGGCAGAAAGAACTGGCGCAAGCCAACGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCAAGGCAGTGGCGGATCAATATGCCGACTTCCTTCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATc                                    >  2:35711/1‑249 (MQ=255)                                                                                                                                                                                        gCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGAATTCCTGCGCGAAGGTTTACACGTTGTCACGCCGAACAAAAAGGCCAACACCACGACGATGGATtatcaccaacagtagc                              >  1:134109/1‑237 (MQ=255)                                                                                                                                                                                                  gTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATgcggc                    <  1:175253/251‑1 (MQ=255)                                                                                                                                                                                                    aCATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGCACAAAAAAGCCAACACCTCGGCGGTGGATTACTACCATCAGTTGCGTTATGcggcgg                  >  1:388233/1‑251 (MQ=255)                                                                                                                                                                                                                   ggAACACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTCCTACCATCAGTTGCGTTATGCGGCGGAAAAATCgcggc      <  2:297691/248‑1 (MQ=255)                                                                                                                                                                                                                   ggAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCgcggc      <  1:178455/248‑1 (MQ=255)                                                                                                                                                                                                                    gAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGGTATGCGGCGGAAAAATCGCGGCGTaa  <  2:49764/251‑1 (MQ=255)                                                                                                                                                                                                                     aaaaCTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCCGTTGCGTTATGCGGCGGAAAAATCGCGGCGTaa  >  1:230507/1‑250 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               GACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAA  >  NC_000913/1680‑2141

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