Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,223,121	A→G	V14A (GTA→GCA)	lysC ←	lysine‑sensitive aspartokinase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,223,121	0	A	G	100.0%	145.7 / NA	39	V14A (GTA→GCA)	lysC	lysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/20);  total (19/20)

TCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTACGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCGCTACTATTT  >  NZ_CP009273/4223004‑4223252                                                                                                                      |                                                                                                                                    tCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTg                                                                                                                 >  2:322872/1‑138 (MQ=255)             aCCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAAttt                                                                                                    >  1:448119/1‑139 (MQ=255)               cagcagATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCa                                                                                                  <  2:142315/139‑1 (MQ=255)                 gcagATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGa                                                                                                <  1:95433/139‑1 (MQ=255)                   agATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACa                                                                                              <  2:448119/139‑1 (MQ=255)                      ttAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATaa                                                                                           <  2:375826/139‑1 (MQ=255)                        aGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAAc                                                                                          <  1:431420/138‑1 (MQ=255)                            aTACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCt                                                                                     >  2:323840/1‑139 (MQ=255)                            aTACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCt                                                                                     >  2:38130/1‑139 (MQ=255)                            aTACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCt                                                                                     >  2:531288/1‑139 (MQ=255)                               ccAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGa                                                                                                <  2:251301/125‑1 (MQ=255)                               ccAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGa                                                                                                >  1:251301/1‑125 (MQ=255)                                         ccGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGAt                                                                        <  1:471824/139‑1 (MQ=255)                                                            cACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGcac                                                     <  1:292114/139‑1 (MQ=255)                                                               gTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAg                                                  >  2:312540/1‑139 (MQ=255)                                                                      tCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAagag                                           <  1:309661/139‑1 (MQ=255)                                                                        aGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCg                                         <  2:113238/139‑1 (MQ=255)                                                                         gAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCg                                         >  1:508321/1‑138 (MQ=255)                                                                             gCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGa                                    <  1:57271/139‑1 (MQ=255)                                                                              cacaATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGa                                                                                                <  1:450798/78‑1 (MQ=255)                                                                              cacaATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGa                                                                                                >  2:450798/1‑78 (MQ=255)                                                                                  atatCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGt                               <  1:533421/139‑1 (MQ=255)                                                                                   tatCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCt                                                           >  1:257966/1‑110 (MQ=255)                                                                                   tatCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCt                                                           <  2:257966/110‑1 (MQ=255)                                                                                        gcgcTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATc                                                                       >  1:458697/1‑93 (MQ=255)                                                                                        gcgcTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATc                                                                       <  2:458697/93‑1 (MQ=255)                                                                                        gcgcTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGgcg                          >  2:159985/1‑138 (MQ=255)                                                                                        gcgcTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGgcg                          >  1:444743/1‑138 (MQ=255)                                                                                         cgcTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGgcg                          >  2:89520/1‑137 (MQ=255)                                                                                           cTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCaga                      <  1:322872/139‑1 (MQ=255)                                                                                           cTGCGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCa                        >  2:371144/1‑137 (MQ=255)                                                                                              cGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATc                                                                       >  2:239349/1‑87 (MQ=255)                                                                                              cGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATc                                                                       <  1:239349/87‑1 (MQ=255)                                                                                              cGGTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGc                    >  2:438127/1‑138 (MQ=255)                                                                                               ggTTCATGGCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGa                  <  2:508321/139‑1 (MQ=255)                                                                                                       gCGTCAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCg           <  1:312540/138‑1 (MQ=255)                                                                                                           cAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCGctac       >  1:358632/1‑138 (MQ=255)                                                                                                           cAAAATCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCGctac       >  2:305714/1‑138 (MQ=255)                                                                                                               aTCAGCTGCGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCGCTACTAttt  <  1:424764/139‑1 (MQ=255)                                                                                                                      |                                                                                                                                    TCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGACAACTAAACGCACGTTGGCATCAGAAAGCACAATATCAGCGCTGCGGTTCATGGCGTCAAAATCAGCTACGCTGGTACCGCCAAATTTGGAGACAACAATTTCAGACATAACTACCTCGTGTCAGGGGATCCATTTTCAGCCTTGGCACAAGGGAAGAGCGGAAGACGGGTGGGCGCAGAGCGATACTTCGCTACTATTT  >  NZ_CP009273/4223004‑4223252

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