Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,131,924	C→T	S697L (TCG→TTG)	metL →	fused aspartate kinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,131,924	0	C	T	100.0%	209.4 / NA	57	S697L (TCG→TTG)	metL	fused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/35);  total (22/35)

AAGACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATCGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGC  >  NC_000913/4131846‑4131995                                                                               |                                                                        aaGACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGc                                                                       >  1:837881/1‑81 (MQ=255)     cGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGt                                                                   >  2:837862/1‑81 (MQ=255)         aTGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCt                                                               <  1:837793/81‑1 (MQ=255)         aTGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCt                                                               <  1:837791/81‑1 (MQ=255)         aTGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCt                                                               <  2:837792/81‑1 (MQ=255)            cgcAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGc                                                             <  2:837794/80‑1 (MQ=255)             gcAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTc                                                           <  2:837796/81‑1 (MQ=255)             gcAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTc                                                           <  2:837795/81‑1 (MQ=255)               aaGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCAt                                                         <  2:837798/81‑1 (MQ=255)               aaGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCAt                                                         <  2:837797/81‑1 (MQ=255)                aGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCAtt                                                        <  1:837799/81‑1 (MQ=255)                  cttgTGATTCTGGCGCGTGAAGCAGGTTACAACATCTAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGc                                                      <  1:837801/78‑1 (MQ=255)                  cTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGc                                                      <  1:837800/81‑1 (MQ=255)                     gTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCg                                                     <  1:837803/79‑1 (MQ=255)                     gTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCg                                                     >  1:837880/1‑79 (MQ=255)                     gTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGaa                                                   <  1:837802/81‑1 (MQ=255)                     gTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGa                                                    <  2:837804/80‑1 (MQ=255)                      tggttCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAg                                                  <  1:837805/78‑1 (MQ=255)                        aTTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAggc                                                <  2:837806/81‑1 (MQ=255)                           cTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAggcgg                                              <  1:837807/80‑1 (MQ=255)                              gcgcGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAggcggc                                             <  2:837809/78‑1 (MQ=255)                              gcgcGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgcag                                           >  1:837891/1‑80 (MQ=255)                              gcgcGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgca                                            <  1:837808/79‑1 (MQ=255)                               cgcgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgcagca                                         >  1:837868/1‑81 (MQ=255)                               cgcgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgcagca                                         >  1:837876/1‑81 (MQ=255)                                 cgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATc                                       >  2:837899/1‑81 (MQ=255)                                 cgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATc                                       <  2:837811/81‑1 (MQ=255)                                 cgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATc                                       <  1:837810/81‑1 (MQ=255)                                   tGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGa                                     >  2:837895/1‑81 (MQ=255)                                   tGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGa                                     >  2:837877/1‑81 (MQ=255)                                     aaGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGAcc                                   <  2:837812/81‑1 (MQ=255)                                     aaGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGAcc                                   <  1:837813/81‑1 (MQ=255)                                      aGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCa                                  >  2:837870/1‑81 (MQ=255)                                         aGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAttt                               >  2:837879/1‑81 (MQ=255)                                            ttACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCtt                            >  1:837865/1‑81 (MQ=255)                                                aacaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGa                         >  2:837887/1‑80 (MQ=255)                                                 acaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGaaa                       <  1:837814/81‑1 (MQ=255)                                                 acaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGaaa                       <  1:837815/81‑1 (MQ=255)                                                 acaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGaa                        >  2:837897/1‑80 (MQ=255)                                                  caTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGaaaa                      >  1:837888/1‑81 (MQ=255)                                                      gAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGc                  <  2:837816/81‑1 (MQ=255)                                                        aCCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGa                <  1:837817/81‑1 (MQ=255)                                                        aCCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGa                >  2:837892/1‑81 (MQ=255)                                                          cGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAt               >  2:837912/1‑80 (MQ=255)                                                              tCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaact          <  1:837818/81‑1 (MQ=255)                                                               cAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaact          >  1:837890/1‑80 (MQ=255)                                                                aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactga        <  2:837822/81‑1 (MQ=255)                                                                aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactga        <  1:837819/81‑1 (MQ=255)                                                                aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactga        <  1:837821/81‑1 (MQ=255)                                                                aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactg         <  2:837820/80‑1 (MQ=255)                                                                   tACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACg     <  2:837824/81‑1 (MQ=255)                                                                   tACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACg     <  2:837823/81‑1 (MQ=255)                                                                    aCGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactgaac      >  2:837889/1‑79 (MQ=255)                                                                      gtgtGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAg   >  2:837893/1‑80 (MQ=255)                                                                      gtgtGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAg   >  2:837894/1‑80 (MQ=255)                                                                      gtgtGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGc  <  2:837825/81‑1 (MQ=255)                                                                        gtgGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAg   >  1:837883/1‑78 (MQ=255)                                                                               |                                                                        AAGACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATCGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGC  >  NC_000913/4131846‑4131995

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