Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,782,970	G→A	A231V (GCG→GTG)	gpsA ←	glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,782,970	0	G	A	100.0%	77.2 / NA	24	A231V (GCG→GTG)	gpsA	glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/13);  total (11/13)

ATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCA  >  NC_000913/3782846‑3783099                                                                                                                             |                                                                                                                                  aTCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACa                                                                                                                    <  1:536597/140‑1 (MQ=255)  cccATGCCCTGACCGAGTATCACGCCAAAACGGCGGTTAAGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACAt                                                                                                                   <  2:248505/139‑1 (MQ=255)      tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGTAAGCACCAGAACGCCAAACCCCGCCCTGCCCATAAAAGGGGGAGGGGCGGCGCCCAGCGCCCCACCAAAACGCGCCAttttc                                                                                                               >  2:533353/1‑138 (MQ=255)        cccTGCCCGAGCATCATGCCAAACCGGGGGTTCCGCGACGGGTTGTGGGTACAGGTAAGCCCCAGACCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTcagc                                                                                                             <  2:139276/140‑1 (MQ=255)        cccTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTcagc                                                                                                             <  1:533353/140‑1 (MQ=255)           tGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTcagccag                                                                                                          <  2:117196/140‑1 (MQ=255)             aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCACCAGAACGCCCAGCCCCCCCCCGCCCCCAAAAGGGGGAGGGGCGGCACCCCCCGCCCCCCCAAGACACGACATTTCagccagcc                                                                                                        >  2:87793/1‑140 (MQ=255)               cGAGCATCATGCAAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGACCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                      <  2:573672/140‑1 (MQ=255)                             aaaCGGCGGTTACGCGACTGGGTGTCGGGACAGGTAAGCACCAAAACGCCAAACCCCCCCACGCCCCTAAAAGTGGCAGGGGCGGCACCCCGCGCCACACCAAAACGCGACATTTTAGCCAGCCCCCGGGGGATCAgcgc                                                                                        >  2:496412/1‑140 (MQ=255)                                   cggTTGCGCGACTGGTTGCCGTTACGGGTAAGCGCCATAGGTCTAGGCACCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcg                                                                                  <  1:115379/140‑1 (MQ=255)                                                       ggTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                              <  1:34304/140‑1 (MQ=255)                                                              gTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTc                                                       >  1:52578/1‑140 (MQ=255)                                                              gTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTc                                                       >  2:635887/1‑140 (MQ=255)                                                                 aGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                    >  2:636792/1‑140 (MQ=255)                                                                  cccccAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGAc                                                   <  2:44026/137‑1 (MQ=255)                                                                   cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGTTGATCAGCCCCGTACGGGCATTCGCACCAAAACCGATACCGTCGGACa                                                  >  2:280072/1‑140 (MQ=255)                                                                   cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACa                                                  >  2:6353/1‑140 (MQ=255)                                                                       aGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                              <  1:635887/140‑1 (MQ=255)                                                                        gATCGCCAAGCCCCGTCATGCCAATCAAGCTGGCAGGCTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATcccc                                             <  1:175554/140‑1 (MQ=255)                                                                        gATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATcccc                                             >  2:153755/1‑140 (MQ=255)                                                                         aTCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                            <  1:431627/140‑1 (MQ=255)                                                                                  ccccGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATgg                                   >  1:414431/1‑140 (MQ=255)                                                                                     cGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCaa                                >  1:331097/1‑140 (MQ=255)                                                                                               aTAAAGGTGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGttttt                      >  1:164877/1‑140 (MQ=255)                                                                                                      tGGCAGGGTCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg               >  2:619880/1‑140 (MQ=255)                                                                                                              tCGGCACCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAg       <  1:619880/140‑1 (MQ=255)                                                                                                                   aCCCAGCGCCACACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTGTTCACCGCGCCGCCAAGCTGCa  >  1:243276/1‑140 (MQ=255)                                                                                                                             |                                                                                                                                  ATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCA  >  NC_000913/3782846‑3783099

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