Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,006,415	Δ1 bp	coding (464/543 nt)	zapC →	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,006,413	0	G	.	100.0%	191.7 / NA	45	coding (462/543 nt)	zapC	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (25/20);  total (25/20)

GCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTGCC  >  NC_000913/1006280‑1006547                                                                                                                                      |                                                                                                                                       gCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcgac                                                                                                                                   >  1:282194/1‑135 (MQ=255)  cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACg                                                                                                                                 <  1:365269/139‑1 (MQ=255)        ggCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTa                                                                                                                           <  2:282194/139‑1 (MQ=255)                   ggTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                                    >  1:620268/1‑117 (MQ=255)                   ggTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                                    <  2:620268/120‑4 (MQ=255)                        aGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAg                                                                                                           >  2:726808/1‑139 (MQ=255)                                 cACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAAcc                                                                                                  <  1:726808/139‑1 (MQ=255)                                  aCGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                                 >  1:1252706/1‑139 (MQ=255)                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAATTGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                            >  1:403207/1‑139 (MQ=255)                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                            >  1:1127978/1‑139 (MQ=255)                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                            >  2:246490/1‑139 (MQ=255)                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                            >  2:570582/1‑139 (MQ=255)                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                              <  1:229595/118‑1 (MQ=255)                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                              >  2:229595/1‑118 (MQ=255)                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTaa                                                                                         >  2:860658/1‑139 (MQ=255)                                                    tGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGc                                                                               >  1:446032/1‑139 (MQ=255)                                                            gAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCt                                                                       <  1:860658/139‑1 (MQ=255)                                                             aaTCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCAt                                                                                                                    >  2:1149988/1‑93 (MQ=255)                                                             aaTCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCAt                                                                                                                    <  1:1149988/93‑1 (MQ=255)                                                                                  tATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                                 <  2:446032/139‑1 (MQ=255)                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAg                                         >  1:224998/1‑139 (MQ=255)                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAg                                                                >  2:869541/1‑114 (MQ=255)                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAg                                                                <  1:869541/114‑1 (MQ=255)                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAAcc                                                                                                  <  1:1169409/69‑1 (MQ=255)                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAAcc                                                                                                  >  2:1169409/1‑69 (MQ=255)                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                            >  1:1081272/1‑139 (MQ=255)                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                            >  2:923537/1‑139 (MQ=255)                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTa                           <  2:224998/139‑1 (MQ=255)                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTa                           <  2:229237/139‑1 (MQ=255)                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTa                           >  2:473538/1‑139 (MQ=255)                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTa                           <  1:494302/134‑1 (MQ=255)                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTa                           >  2:494302/1‑134 (MQ=255)                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                      >  1:371225/1‑139 (MQ=255)                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                      >  1:1265746/1‑139 (MQ=255)                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                      >  1:214739/1‑139 (MQ=255)                                                                                                               aTTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCt                    >  1:1196426/1‑139 (MQ=255)                                                                                                               aTTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCt                    <  1:246490/139‑1 (MQ=255)                                                                                                                     ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                             <  1:642458/60‑1 (MQ=255)                                                                                                                     ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                             >  2:642458/1‑60 (MQ=255)                                                                                                                     ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCa              >  1:188562/1‑139 (MQ=255)                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                                 <  1:833234/51‑1 (MQ=255)                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                                 >  2:833234/1‑51 (MQ=255)                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa         >  2:332722/1‑139 (MQ=255)                                                                                                                               tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg    <  2:371225/139‑1 (MQ=255)                                                                                                                               tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg    <  1:473538/139‑1 (MQ=255)                                                                                                                               tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg    <  2:214739/139‑1 (MQ=255)                                                                                                                               tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg    <  2:188562/139‑1 (MQ=255)                                                                                                                                 cAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTGcc  <  2:1265746/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTGCC  >  NC_000913/1006280‑1006547

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