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%	200.4 / NA	49	coding (462/543 nt)	zapC	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/27);  total (22/27)

GGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTG  >  NC_000913/1006279‑1006545                                                                                                                                       |                                                                                                                                     ggCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                                  >  1:1611034/1‑136 (MQ=255)   cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACg                                                                                                                               >  2:1649286/1‑139 (MQ=255)   cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACg                                                                                                                               <  2:271286/139‑1 (MQ=255)      tGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                            >  1:1022929/1‑139 (MQ=255)      tGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                            >  1:1202884/1‑139 (MQ=255)                      ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                            <  2:1202884/139‑1 (MQ=255)                      ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCAGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                            <  2:1652451/139‑1 (MQ=255)                         aGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAg                                                                                                         >  1:19736/1‑139 (MQ=255)                              tACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGa                                                                                                    <  2:111476/139‑1 (MQ=255)                              tACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGa                                                                                                    <  2:1022929/139‑1 (MQ=255)                                   aCGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               <  2:1195522/139‑1 (MQ=255)                                        cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                          >  2:1117986/1‑139 (MQ=255)                                           gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               <  2:386343/131‑1 (MQ=255)                                           gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               >  1:386343/1‑131 (MQ=255)                                           gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTaa                                                                                       >  2:562455/1‑139 (MQ=255)                                                 ttGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                                 <  2:19736/139‑1 (MQ=255)                                                                     cGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg                                                             <  1:1649286/139‑1 (MQ=255)                                                                         aaTGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAg                                                          <  2:1665426/138‑1 (MQ=255)                                                                         aaTGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAg                                                          >  1:1665426/1‑138 (MQ=255)                                                                              cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAc                                                    >  2:1605678/1‑139 (MQ=255)                                                                                  cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACgcgc                                                <  1:1767235/139‑1 (MQ=255)                                                                                   tATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                               <  2:1450563/139‑1 (MQ=255)                                                                                     tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                                  >  1:1132464/1‑52 (MQ=255)                                                                                     tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                                  <  2:1132464/55‑4 (MQ=255)                                                                                     tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACgcgc                                                >  1:659050/1‑136 (MQ=255)                                                                                     tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACgcgc                                                <  2:659050/136‑1 (MQ=255)                                                                                           ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                                <  1:887050/98‑1 (MQ=255)                                                                                           ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                                >  2:887050/1‑98 (MQ=255)                                                                                               cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCt                                                                            >  2:95723/1‑98 (MQ=255)                                                                                               cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCt                                                                            <  1:95723/98‑1 (MQ=255)                                                                                               cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  <  2:381448/124‑1 (MQ=255)                                                                                               cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  >  1:381448/1‑124 (MQ=255)                                                                                               cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  2:552072/1‑139 (MQ=255)                                                                                                 caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                                 <  1:1117986/139‑1 (MQ=255)                                                                                                    caacgcagggttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                                 >  2:724373/10‑136 (MQ=255)                                                                                                    caacgcagggttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                                 <  1:724373/127‑1 (MQ=255)                                                                                                              ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   <  2:423204/124‑1 (MQ=255)                                                                                                              ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  1:423204/1‑124 (MQ=255)                                                                                                                 ttGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:106999/139‑1 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                           >  1:259240/1‑60 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                           <  2:259240/60‑1 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                               >  1:40473/1‑104 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                               <  2:40473/104‑1 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:1078754/134‑1 (MQ=255)                                                                                                                      ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 >  1:1078754/1‑134 (MQ=255)                                                                                                                        tCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTac                        >  1:690738/1‑125 (MQ=255)                                                                                                                        tCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTac                        <  2:690738/125‑1 (MQ=255)                                                                                                                        tCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCaaa          >  2:1892271/1‑139 (MQ=255)                                                                                                                           cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               <  1:1572304/51‑1 (MQ=255)                                                                                                                           cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               >  2:1572304/1‑51 (MQ=255)                                                                                                                               aTGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGt   <  1:552072/139‑1 (MQ=255)                                                                                                                                tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg  <  2:1611034/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     GGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTG  >  NC_000913/1006279‑1006545

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