Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,966,188	C→T	*749* (TAG→TAA)	ptsP ←	PEP‑protein phosphotransferase enzyme I, GAF domain containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,966,188	0	C	T	100.0%	96.5 / NA	31	*749* (TAG→TAA)	ptsP	PEP‑protein phosphotransferase enzyme I, GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/12);  total (19/12)

TTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCG  >  NC_000913/2966061‑2966322                                                                                                                                |                                                                                                                                       ttCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTcc                                                                                                                             >  1:610782/1‑139 (MQ=255) ttCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTcc                                                                                                                             >  1:265373/1‑139 (MQ=255)                      ttGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGaa                                                                                                                        >  2:1393890/1‑123 (MQ=255)                      ttGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGaa                                                                                                                        <  1:1393890/123‑1 (MQ=255)                       tGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcg                                                                                                       <  2:265373/139‑1 (MQ=255)                           aTGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGAc                                                                                                   <  2:539892/139‑1 (MQ=255)                               gcgcCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTc                                                                                               >  1:1221839/1‑139 (MQ=255)                                          ttCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                    >  1:1451034/1‑139 (MQ=255)                                                           tGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACt                                                                   >  2:430887/1‑139 (MQ=255)                                                              gAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCg                                                                <  1:1431676/139‑1 (MQ=255)                                                                       cacaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGt                                                       >  2:894356/1‑139 (MQ=255)                                                                       cacaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGt                                                       >  2:1442366/1‑139 (MQ=255)                                                                         caAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                     <  2:1451034/139‑1 (MQ=255)                                                                            ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                     <  1:158523/136‑1 (MQ=255)                                                                            ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                     >  2:158523/1‑136 (MQ=255)                                                                                          ttGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCt                                    >  1:77467/1‑139 (MQ=255)                                                                                           tGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTg                                   >  2:443791/1‑139 (MQ=255)                                                                                               ggATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCa                               >  1:1098291/1‑139 (MQ=255)                                                                                                           aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                   <  2:610782/139‑1 (MQ=255)                                                                                                           aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                   <  1:1442366/139‑1 (MQ=255)                                                                                                                       tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGt                                                              <  1:392257/84‑1 (MQ=255)                                                                                                                       tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGt                                                              >  2:392257/1‑84 (MQ=255)                                                                                                                       tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                   <  2:1167308/127‑1 (MQ=255)                                                                                                                       tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                   >  1:1167308/1‑127 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc              <  1:314688/128‑1 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc              >  2:314688/1‑128 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGATTGGGAATTATCAATgc   >  2:1109410/1‑139 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                          <  1:1085810/51‑1 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                          >  2:1085810/1‑51 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg  >  2:1615063/1‑139 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg  >  2:922374/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       TTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCG  >  NC_000913/2966061‑2966322

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