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%	189.3 / NA	52	S697L (TCG→TTG)	metL	fused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/26);  total (26/26)

GACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATCGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGCAGATG  >  NC_000913/4131848‑4132000                                                                             |                                                                             gACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTg                                                                          >  1:770052/1‑81 (MQ=255) gACGTGATGCGCAAGCTGGTGATTATGGCGAGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTg                                                                          <  2:769979/81‑1 (MQ=255)  aCGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTgg                                                                         <  2:769980/81‑1 (MQ=255)   cGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGt                                                                        >  1:770046/1‑81 (MQ=255)     tGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGc                                                                      >  1:770035/1‑81 (MQ=255)       aTGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCt                                                                    <  2:769981/81‑1 (MQ=255)          cgcAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCt                                                                 <  2:769982/81‑1 (MQ=255)           gcAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCt                                                                 >  1:770044/1‑80 (MQ=255)             aaGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGCGTGGAATTGCTGGTGCCTGCTCAt                                                              >  1:770039/1‑81 (MQ=255)              aGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCAtt                                                             <  2:769983/81‑1 (MQ=255)                   gTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGa                                                         <  1:769985/80‑1 (MQ=255)                    tGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAg                                                       <  2:769985/81‑1 (MQ=255)                       ttCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAggcg                                                    <  1:769987/81‑1 (MQ=255)                       ttCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAggcg                                                    >  2:770061/1‑81 (MQ=255)                            gcgcGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgcagc                                               <  1:769989/81‑1 (MQ=255)                            gcgcGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGgcag                                                <  1:769988/80‑1 (MQ=255)                               cgTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATc                                            >  1:770054/1‑81 (MQ=255)                                gTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATc                                            >  2:770065/1‑80 (MQ=255)                                gTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCg                                           <  1:769990/81‑1 (MQ=255)                                 tGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGa                                          <  2:769990/81‑1 (MQ=255)                                  gAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGAc                                         >  1:770048/1‑81 (MQ=255)                                  gAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGAc                                         <  2:769991/81‑1 (MQ=255)                                   aaGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGAc                                         <  2:769992/80‑1 (MQ=255)                                    aGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCa                                       >  2:770057/1‑81 (MQ=255)                                     gCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAt                                      <  1:769994/81‑1 (MQ=255)                                     gCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAt                                      <  2:769994/81‑1 (MQ=255)                                      cAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAtt                                     >  1:770063/1‑81 (MQ=255)                                       aGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAttt                                    >  1:770047/1‑81 (MQ=255)                                       aGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCAttt                                    <  2:769995/81‑1 (MQ=255)                                            acaacaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCttt                                >  1:770077/1‑80 (MQ=255)                                            acaacaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCttt                                >  1:770055/1‑80 (MQ=255)                                             caacaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTg                               >  2:770050/1‑80 (MQ=255)                                              aacaTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGaa                             <  2:769996/81‑1 (MQ=255)                                                 aTCGAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAAt                          <  1:769998/81‑1 (MQ=255)                                                    gAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGc                       >  2:770075/1‑81 (MQ=255)                                                    gAACCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGc                       <  2:769998/81‑1 (MQ=255)                                                     aaCCGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGc                       >  2:770068/1‑80 (MQ=255)                                                       ccGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAt                    >  2:770063/1‑81 (MQ=255)                                                       ccGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAt                    >  1:770050/1‑81 (MQ=255)                                                       ccGGATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAt                    >  2:770071/1‑81 (MQ=255)                                                         ggATCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtg                   >  1:770078/1‑80 (MQ=255)                                                           aTCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaac                <  1:770000/81‑1 (MQ=255)                                                            tCAGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaact               >  2:770088/1‑81 (MQ=255)                                                              aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactga             <  2:770000/81‑1 (MQ=255)                                                              aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactg              >  2:770101/1‑80 (MQ=255)                                                              aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactg              <  2:770001/80‑1 (MQ=255)                                                              aGGTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactg              <  1:770003/80‑1 (MQ=255)                                                               ggTACGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGAtgaactgaa            <  2:770003/81‑1 (MQ=255)                                                                   cGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAg        >  2:770066/1‑81 (MQ=255)                                                                   cGTGTGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAg        >  1:770084/1‑81 (MQ=255)                                                                    gtgtGGAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGc       <  1:770005/81‑1 (MQ=255)                                                                         gAATTGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGCAGATg  >  1:770068/1‑81 (MQ=255)                                                                             |                                                                             GACGTGATGCGCAAGCTGGTGATTCTGGCGCGTGAAGCAGGTTACAACATCGAACCGGATCAGGTACGTGTGGAATCGCTGGTGCCTGCTCATTGCGAAGGCGGCAGCATCGACCATTTCTTTGAAAATGGCGATGAACTGAACGAGCAGATG  >  NC_000913/4131848‑4132000

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: ‑