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%	122.1 / NA	39	*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 (25/14);  total (25/14)

TGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966055‑2966323                                                                                                                                      |                                                                                                                                        tGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATaa                                                                                                                                    >  1:736914/1‑139 (MQ=255)           cTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGa                                                                                                                          <  2:510182/139‑1 (MQ=255)                agagAACAGCGATCGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCGGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCAAATCAg                                                                                                                     >  2:1129186/1‑139 (MQ=255)                  agaACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAgcc                                                                                                                   <  2:736914/139‑1 (MQ=255)                                         cgcAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATa                                                                                            >  2:824996/1‑139 (MQ=255)                                                   aGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCaa                                                                                  >  2:730453/1‑139 (MQ=255)                                                     gcgcTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAAcc                                                                                <  1:730453/139‑1 (MQ=255)                                                         tCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   <  2:660903/116‑1 (MQ=255)                                                         tCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   >  1:660903/1‑116 (MQ=255)                                                          ccAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATg                                                                           >  1:783206/1‑139 (MQ=255)                                                                   cGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTc                                                                  >  2:185567/1‑139 (MQ=255)                                                                             cacaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGt                                                        >  2:404027/1‑139 (MQ=255)                                                                                  ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  <  2:57071/92‑1 (MQ=255)                                                                                  ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  >  1:57071/1‑92 (MQ=255)                                                                                    ggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAgcccgcc                                                                                                               >  2:411759/1‑77 (MQ=255)                                                                                    ggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAgcccgcc                                                                                                               <  1:411759/77‑1 (MQ=255)                                                                                    ggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 >  2:1205559/1‑139 (MQ=255)                                                                                                                 aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  2:783206/139‑1 (MQ=255)                                                                                                                     tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              <  1:129696/77‑1 (MQ=255)                                                                                                                     tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              >  2:129696/1‑77 (MQ=255)                                                                                                                     tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCgg                >  1:576500/1‑139 (MQ=255)                                                                                                                       tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   <  1:264478/54‑1 (MQ=255)                                                                                                                       tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   >  2:264478/1‑54 (MQ=255)                                                                                                                       tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGt                                                        >  1:1140364/1‑97 (MQ=255)                                                                                                                       tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGt                                                        <  2:1140364/97‑1 (MQ=255)                                                                                                                       tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAAATTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCa              >  2:648545/1‑139 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcgcg                                                                                                      >  2:41941/1‑44 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcgcg                                                                                                      <  1:41941/44‑1 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGCACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTAGGCATAATCaa       >  1:937446/1‑139 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc               <  1:909745/131‑1 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc               >  2:909745/1‑131 (MQ=255)                                                                                                                              gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCaa       >  1:582295/1‑139 (MQ=255)                                                                                                                                tCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATg     >  2:921201/1‑139 (MQ=255)                                                                                                                                 ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  2:291965/1‑89 (MQ=255)                                                                                                                                 ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  1:291965/89‑1 (MQ=255)                                                                                                                                 ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgc    >  2:57167/1‑139 (MQ=255)                                                                                                                                  cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              >  2:464212/1‑64 (MQ=255)                                                                                                                                  cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTg                                                                              <  1:464212/64‑1 (MQ=255)                                                                                                                                  cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg   >  1:948145/1‑139 (MQ=255)                                                                                                                                   gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAAcc                                                                                >  2:719116/1‑61 (MQ=255)                                                                                                                                   gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACtt                                                                   <  1:790105/74‑1 (MQ=255)                                                                                                                                   gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCCCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTACACAAGATATTCTGCTTCGGCATAATCAATgcgc  >  1:1161970/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        TGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966055‑2966323

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