Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,966,553	C→G	A628P (GCC→CCC)	ptsP ←	phosphoenolpyruvate‑protein phosphotransferase PtsP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,966,553	0	C	G	100.0%	104.7 / NA	29	A628P (GCC→CCC)	ptsP	phosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (22/7);  total (22/7)

ATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAA  >  NC_000913/2966477‑2966628                                                                             |                                                                            atcGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAg                                                                          >  1:538484/1‑80 (MQ=255)     ccAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGt                                                                     <  2:538434/81‑1 (MQ=255)            tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGt                                                              >  1:538472/1‑81 (MQ=255)              gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGt                                                              >  1:538468/1‑79 (MQ=255)                  cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGt                                                              >  1:538455/1‑75 (MQ=255)                   aTTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCg                                                       >  2:538462/1‑81 (MQ=255)                                gACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCaa                                             >  1:538495/1‑78 (MQ=255)                                  cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACa                                           >  1:538493/1‑78 (MQ=255)                                    gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGAt                                      >  2:538467/1‑81 (MQ=255)                                     tCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATg                                     <  2:538435/81‑1 (MQ=255)                                       aTAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGa                                    <  1:538436/80‑1 (MQ=255)                                         aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAAt                                 >  2:538481/1‑81 (MQ=255)                                          aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATc                                >  2:538497/1‑81 (MQ=255)                                          aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATc                                >  1:538483/1‑81 (MQ=255)                                           aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATc                                >  2:538482/1‑80 (MQ=255)                                               tCGCCACCCGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGAccc                           >  2:538502/1‑81 (MQ=255)                                                     cccGGGTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTg                     >  1:538476/1‑81 (MQ=255)                                                         ggTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                 <  2:538438/81‑1 (MQ=255)                                                         ggTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                 <  1:538437/81‑1 (MQ=255)                                                          gTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGc                    >  1:538479/1‑77 (MQ=255)                                                          gTATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGc                    >  2:538499/1‑77 (MQ=255)                                                           tATTGTTGCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGAt               <  1:538439/81‑1 (MQ=255)                                                               gttgCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgg           >  2:538474/1‑81 (MQ=255)                                                                ttgCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgg           >  1:538475/1‑80 (MQ=255)                                                                 tgCGATCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgg           >  2:538490/1‑79 (MQ=255)                                                                     aTCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCAt     <  1:538440/81‑1 (MQ=255)                                                                      tCAACGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAAACGACCCGCTTTGCCAGATGCGGCAGCAt     >  1:538494/1‑80 (MQ=255)                                                                        aaCGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATaaa  >  1:538503/1‑81 (MQ=255)                                                                        aaCGGGCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATaaa  >  2:538486/1‑81 (MQ=255)                                                                             |                                                                            ATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAA  >  NC_000913/2966477‑2966628

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