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%	108.8 / NA	35	*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 (20/15);  total (20/15)

TTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966061‑2966323                                                                                                                                |                                                                                                                                        ttCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCATTATAACCCTcc                                                                                                                              >  1:996020/1‑139 (MQ=255) ttCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCAGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTcc                                                                                                                              <  1:1179438/139‑1 (MQ=255)                  gCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCa                                                                                                             >  1:534539/1‑139 (MQ=255)                   cGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCAt                                                                                                            >  2:788484/1‑139 (MQ=255)                            tGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   <  1:1165112/139‑1 (MQ=255)                                                  cTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTcc                                                                                                                              <  2:1227965/90‑1 (MQ=255)                                                  cTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTcc                                                                                                                              >  1:1227965/1‑90 (MQ=255)                                                             cGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTta                                                                                                                                       >  2:238983/1‑70 (MQ=255)                                                             cGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTc                                                                  >  2:1430010/1‑139 (MQ=255)                                                                         caAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:323927/139‑1 (MQ=255)                                                                         caAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:534539/139‑1 (MQ=255)                                                                            ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTgcgc                                                   >  1:1582863/1‑139 (MQ=255)                                                                            ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTgcgc                                                   <  1:788484/139‑1 (MQ=255)                                                                                               ggATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGa                                                                                                     <  1:49304/70‑1 (MQ=255)                                                                                               ggATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGa                                                                                                     >  2:49304/1‑70 (MQ=255)                                                                                                     cGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  1:1017122/111‑1 (MQ=255)                                                                                                     cGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  2:1017122/1‑111 (MQ=255)                                                                                                       ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 <  1:1370683/114‑1 (MQ=255)                                                                                                       ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 >  2:1370683/1‑114 (MQ=255)                                                                                                       ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTc                        >  1:1576807/1‑139 (MQ=255)                                                                                                       ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTc                        >  1:412536/1‑139 (MQ=255)                                                                                                           aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  2:1582863/139‑1 (MQ=255)                                                                                                           aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  1:655656/139‑1 (MQ=255)                                                                                                                         aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAAt      >  1:1019496/1‑139 (MQ=255)                                                                                                                          tCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  >  1:850264/1‑46 (MQ=255)                                                                                                                          tCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  <  2:850264/46‑1 (MQ=255)                                                                                                                          tCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                           >  1:778395/1‑53 (MQ=255)                                                                                                                          tCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                           <  2:778395/53‑1 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:1352915/89‑1 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  1:1352915/1‑89 (MQ=255)                                                                                                                           ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgc    >  2:1543706/1‑139 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              >  2:852785/1‑64 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              <  1:852785/64‑1 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg   >  1:637737/1‑139 (MQ=255)                                                                                                                            cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg   >  2:850973/1‑139 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  2:445043/1‑87 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  1:445043/87‑1 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc               <  2:1225092/126‑1 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc               >  1:1225092/1‑126 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  <  2:412536/139‑1 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  >  1:845251/1‑139 (MQ=255)                                                                                                                             gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  >  1:222791/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                        TTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966061‑2966323

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