Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,801	G→A	T66I (ACA→ATA)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,801	0	G	A	100.0%	107.3 / NA	34	T66I (ACA→ATA)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/11);  total (23/11)

CGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTT  >  NC_000913/3379666‑3379921                                                                                                                                        |                                                                                                                         cGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTatgt                                                                                                                       >  2:351698/1‑139 (MQ=255)   cagcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCt                                                                                                                     >  2:177394/1‑139 (MQ=255)          aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGtt                                                                                                              <  2:378689/139‑1 (MQ=255)              ggTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTAcc                                                                                                          <  2:371118/139‑1 (MQ=255)                tGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTAcccc                                                                                                        >  2:454137/1‑139 (MQ=255)                                gTTTCCGCTCTCCCGGCAGTCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACCCCGCCCCACATTTATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCcg                                                                                        <  1:748336/139‑1 (MQ=255)                                         ctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAg                                                                               >  2:714198/1‑139 (MQ=255)                                                    cTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTcc                                                                    >  1:303651/1‑139 (MQ=255)                                                      ttGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTg                                                                  >  1:40573/1‑139 (MQ=255)                                                               aaGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCtggt                                                         >  2:399109/1‑139 (MQ=255)                                                                   tCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCg                                                     <  2:319589/139‑1 (MQ=255)                                                                     cATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTg                                                   <  1:678686/139‑1 (MQ=255)                                                                             aGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATtgat                                           <  1:635581/139‑1 (MQ=255)                                                                             aGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATtgat                                           >  1:88211/1‑139 (MQ=255)                                                                                                    ccccACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCa                    >  2:793027/1‑139 (MQ=255)                                                                                                     cccACATATATAAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAg                   >  2:553991/1‑139 (MQ=255)                                                                                                             tataAGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGa           <  2:303651/139‑1 (MQ=255)                                                                                                                aaGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTa                                                                                <  2:390400/67‑1 (MQ=255)                                                                                                                aaGCCACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTa                                                                                >  1:390400/1‑67 (MQ=255)                                                                                                                   ccACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCct     >  1:727089/1‑139 (MQ=255)                                                                                                                   ccACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCct     >  1:242303/1‑139 (MQ=255)                                                                                                                    cACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCctc    >  2:505640/1‑139 (MQ=255)                                                                                                                    cACGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCctc    >  1:786598/1‑139 (MQ=255)                                                                                                                     aCGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCctct   >  2:96041/1‑139 (MQ=255)                                                                                                                     aCGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCctct   >  2:16062/1‑139 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCcgc                                                                                       >  2:12443/1‑54 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCcgc                                                                                       <  1:12443/54‑1 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCAc                                                                          <  2:453455/67‑1 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCAc                                                                          >  1:453455/1‑67 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCtt  >  2:322135/1‑139 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCtt  >  2:151732/1‑139 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCtt  >  2:115911/1‑139 (MQ=255)                                                                                                                      cGCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCtt  >  1:793338/1‑139 (MQ=255)                                                                                                                       gCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCca                                                                                                       >  2:276869/1‑37 (MQ=255)                                                                                                                       gCACTGGCGTATGCTTTATGTCTTCGCGTTTACCCca                                                                                                       <  1:276869/37‑1 (MQ=255)                                                                                                                                        |                                                                                                                         CGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTT  >  NC_000913/3379666‑3379921

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: ‑