Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,073,292	C→A	V134F (GTT→TTT)	epd ←	D‑erythrose‑4‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,073,292	0	C	A	100.0%	126.2 / NA	36	V134F (GTT→TTT)	epd	D‑erythrose‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/17);  total (19/17)

CAATTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAACGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCAT  >  NC_000913/3073216‑3073369                                                                             |                                                                              cAATTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTc                                                                           <  1:533609/81‑1 (MQ=255)  aaTTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTcg                                                                          <  2:533611/81‑1 (MQ=255)  aaTTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCC‑CACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTcgc                                                                         <  2:533610/81‑1 (MQ=255)   aTTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTcgc                                                                         >  1:533659/1‑81 (MQ=255)    ttCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTcgcg                                                                        <  1:533612/81‑1 (MQ=255)      cGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGt                                                                       <  2:533613/80‑1 (MQ=255)       gTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTc                                                                      >  1:533656/1‑80 (MQ=255)          gTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCgag                                                                   <  1:533614/80‑1 (MQ=255)                 aaGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTAc                                                           >  2:533678/1‑81 (MQ=255)                     gTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGcc                                                       <  2:533615/81‑1 (MQ=255)                        aGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAg                                                     >  1:533667/1‑80 (MQ=255)                        aGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAg                                                     >  2:533662/1‑80 (MQ=255)                          aaaCGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGAt                                                  >  1:533647/1‑81 (MQ=255)                             cGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGAtgtg                                               <  2:533616/81‑1 (MQ=255)                              gATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGAtgtg                                               >  2:533658/1‑80 (MQ=255)                              gATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGAtgtg                                               >  1:533673/1‑80 (MQ=255)                               aTGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGaa                                             <  1:533617/81‑1 (MQ=255)                                  cGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGaa                                             >  2:533665/1‑78 (MQ=255)                                    gtgtTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAagag                                        <  1:533618/81‑1 (MQ=255)                                         ccGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGAGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACtt                                   >  2:533645/1‑81 (MQ=255)                                          cGCACGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACttt                                  <  1:533619/81‑1 (MQ=255)                                             aCGAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACtttttt                               >  1:533668/1‑81 (MQ=255)                                               gAAGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTgg                             >  1:533681/1‑81 (MQ=255)                                                aaGTTGATCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGc                            >  2:533674/1‑81 (MQ=255)                                                       tCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGCTTTTAAAAGTTCACTTTTTTTGCCCcggcg                     >  1:864502/1‑81 (MQ=255)                                                       tCCTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCcggcg                     >  2:533663/1‑81 (MQ=255)                                                         cTGATTGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCggcggc                   <  2:533620/81‑1 (MQ=255)                                                             ttGACGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAata               <  2:533621/81‑1 (MQ=255)                                                                aCGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATgc            <  1:533622/81‑1 (MQ=255)                                                                 cGCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATgcg           <  1:533623/81‑1 (MQ=255)                                                                  gCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATgcgc          >  2:533672/1‑81 (MQ=255)                                                                  gCCGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGAACTTTTTTGGCCCCGGCGGCAATATgcgc          >  1:533669/1‑81 (MQ=255)                                                                    cGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCtt        <  2:533625/81‑1 (MQ=255)                                                                    cGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCtt        <  2:533624/81‑1 (MQ=255)                                                                    cGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCtt        >  1:533684/1‑81 (MQ=255)                                                                    cGTAAACAAAGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCtt        >  2:533685/1‑81 (MQ=255)                                                                        aacaaaGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGGGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGc     <  2:533626/80‑1 (MQ=255)                                                                        aacaaaGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGc     <  1:533627/80‑1 (MQ=255)                                                                          caaaGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCAt  <  2:533628/81‑1 (MQ=255)                                                                          caaaGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCAt  >  2:533660/1‑81 (MQ=255)                                                                          caaaGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCAt  <  1:533629/81‑1 (MQ=255)                                                                             |                                                                              CAATTCGTGGTACAGGAAGCGTTAGAAACGATGCGGTGTTCCGCACGAAGTTGATCCTGATTGACGCCGTAAACAACGGTCGCGTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCAT  >  NC_000913/3073216‑3073369

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