Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,600,589	G→A	S129S (TCC→TCT)	dnaC ←	DNA biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,600,589	0	G	A	100.0%	133.6 / NA	42	S129S (TCC→TCT)	dnaC	DNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/19);  total (23/19)

TCTCATCGATCACCAGCAGATCAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACGGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGGCGCGCTTTGCTTA  >  NC_000913/4600465‑4600722                                                                                                                             |                                                                                                                                      tctcATCGATCACCAGCAGATCAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTag                                                                                                                         >  2:1661779/1‑139 (MQ=255)                    aTCAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGAtt                                                                                                      >  1:1613426/1‑139 (MQ=255)                    aTCAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGAtt                                                                                                      >  1:1617238/1‑139 (MQ=255)                      cAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTgc                                                                                                    <  1:772799/139‑1 (MQ=255)                                             aGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCgg                                                                             >  1:1220691/1‑139 (MQ=255)                                                tgttCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTc                                                                          >  1:1541556/1‑139 (MQ=255)                                                               ccGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGa                                                             >  2:1163851/1‑137 (MQ=255)                                                               ccGCTATTCCTGAAGGCATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGACTTACCAGa                                                             <  1:1163851/137‑1 (MQ=255)                                                                      tCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGa                                                                                                        <  2:410747/87‑1 (MQ=255)                                                                      tCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGa                                                                                                        >  1:410747/1‑87 (MQ=255)                                                                       ccTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTgccg                                                                                                  >  2:560068/1‑92 (MQ=255)                                                                       ccTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTgccg                                                                                                  <  1:560068/92‑1 (MQ=255)                                                                       ccTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTcc                                                                         <  2:227647/117‑1 (MQ=255)                                                                       ccTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTcc                                                                         >  1:227647/1‑117 (MQ=255)                                                                             ggTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACcgcg                                                                                                                            <  1:387145/60‑1 (MQ=255)                                                                             ggTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACcgcg                                                                                                                            >  2:387145/1‑60 (MQ=255)                                                                                 tCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGt                                         <  2:551975/139‑1 (MQ=255)                                                                                 tCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGt                                         <  1:1661779/139‑1 (MQ=255)                                                                                 tCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGt                                         >  2:1686624/1‑139 (MQ=255)                                                                                  cTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCTGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGtt                                        >  1:1025749/1‑139 (MQ=255)                                                                                  cTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGtt                                        >  2:86857/1‑139 (MQ=255)                                                                                         cGCCGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCa                                 <  2:1613426/139‑1 (MQ=255)                                                                                           ccGACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATATTGCCGTCaaa                               >  1:1126367/1‑139 (MQ=255)                                                                                             gACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                <  1:1402091/88‑1 (MQ=255)                                                                                             gACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                >  2:1402091/1‑88 (MQ=255)                                                                                             gACATAATATCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAAct                             <  2:1220691/139‑1 (MQ=255)                                                                                                   atatCGGCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGa                       >  1:832462/1‑139 (MQ=255)                                                                                                         gCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGtt                                                                           <  1:1421916/81‑1 (MQ=255)                                                                                                         gCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGtt                                                                           >  2:1421916/1‑81 (MQ=255)                                                                                                         gCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACt                 >  2:230631/1‑139 (MQ=255)                                                                                                         gCCACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACt                 >  2:1302273/1‑139 (MQ=255)                                                                                                          ccACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTgg                                                                      <  2:699621/85‑1 (MQ=255)                                                                                                          ccACGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTgg                                                                      >  1:699621/1‑85 (MQ=255)                                                                                                            aCGGTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGgc              <  1:86857/139‑1 (MQ=255)                                                                                                              ggTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                >  1:953124/1‑71 (MQ=255)                                                                                                              ggTGATGATCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                <  2:953124/71‑1 (MQ=255)                                                                                                                 gatgatCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCTGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGGCGCGCt         <  2:1025749/139‑1 (MQ=255)                                                                                                                    gatCAATACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGGCGCGCTTTg      >  1:984469/1‑139 (MQ=255)                                                                                                                        aaTACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                <  2:783284/61‑1 (MQ=255)                                                                                                                        aaTACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGc                                                                                >  1:783284/1‑61 (MQ=255)                                                                                                                        aaTACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGGCGCGCTTTGCTTa  <  2:984469/139‑1 (MQ=255)                                                                                                                          tACAGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGttttt                                                                                  <  1:1709670/57‑1 (MQ=255)                                                                                                                             |                                                                                                                                      TCTCATCGATCACCAGCAGATCAACGTTGCTCAGATCGTTAAGCAGTTGTTCTTCGCTGGTACCGCTATTCCTGAAGGTATCTTTCATCGCCGACATAATATCGGCCACGGTGATGATCAATACGGATTTACCGCGTAGCAGCAGCTCGTTGCAGATTGCCGCCGCCAGATGGTTTTTGCCGGTTCCTGGCTTACCAGAAAAGATGAAGCTGGCGATGTTGCCGTCAAACTCTTCGACATACTGGCGCGCTTTGCTTA  >  NC_000913/4600465‑4600722

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