Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	934,532	G→A	14.8%	E437K (GAA→AAA)	ftsK →	DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,532	0	G	A	14.8%	71.1 / 5.1	27	E437K (GAA→AAA)	ftsK	DNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base A (2/2);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAACCCGTTGAACAGCAGCCTGTTGTGGAGCCTGAACCCGTT  >  NC_000913/934390‑934669                                                                                                                                               |                                                                                                                                          tGCAGTGCAATATAATGGGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGagca                                                                                                                                    <  2:15008/150‑1 (MQ=255)        cAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATcc                                                                                                                             <  2:155019/150‑1 (MQ=255)         aaTATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCa                                                                                                                            >  2:40891/1‑150 (MQ=255)           tataATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACt                                                                                                                          >  2:174899/1‑150 (MQ=255)           tataATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCAAAGAGCAGCAATCCACt                                                                                                                          >  1:110131/1‑150 (MQ=255)           tataATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCAAAGAGCAGCAATCCACt                                                                                                                          >  1:110118/1‑150 (MQ=255)               aTGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTg                                                                                                                      <  2:207771/150‑1 (MQ=255)               aTGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTg                                                                                                                      <  1:213566/150‑1 (MQ=255)                  aGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTc                                                                                                                   <  1:194859/150‑1 (MQ=255)                       cTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAg                                                                                                              >  2:96537/1‑150 (MQ=255)                                         ccGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGTGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAg                                                                                            <  1:226653/150‑1 (MQ=255)                                                  ccGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTAt                                                                                   >  1:18135/1‑150 (MQ=255)                                                                  cAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCATGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTc                                                                   <  2:70808/150‑1 (MQ=255)                                                                  cAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTc                                                                   <  2:70810/150‑1 (MQ=255)                                                                      tGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGa                                                               <  1:246084/150‑1 (MQ=255)                                                                      tGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGa                                                               <  1:689/150‑1 (MQ=255)                                                                      tGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGa                                                               <  1:702/150‑1 (MQ=255)                                                                           aacCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGt                                                          <  1:246056/150‑1 (MQ=255)                                                                             ccTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTg                                                        >  1:72047/1‑150 (MQ=255)                                                                               tGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTa                                                      >  2:159163/1‑150 (MQ=255)                                                                                   cAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCACTCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCaa                                                  <  1:171322/150‑1 (MQ=255)                                                                                   cAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCaa                                                  <  1:45794/150‑1 (MQ=255)                                                                                               tatGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCAAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAAccc                                      <  2:110131/150‑1 (MQ=255)                                                                                               tatGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCAAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAAccc                                      <  2:110118/150‑1 (MQ=255)                                                                                                                        tGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAACCCGTTGAACAGCAGCCTGTTGTGGAGc             <  1:174899/150‑1 (MQ=255)                                                                                                                                   gCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAACCCGTTGAACAGCAGCCTGTTGTGGAGCCTGAACCCGtt  >  1:236161/1‑150 (MQ=255)                                                                                                                                   gCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAACCCGTTGAACAGCAGCCTGTTGTGGAGCCTGAACCCGtt  >  1:236162/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                          TGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCCTGGCAAGCCGAAGAGCAGCAATCCACTTTTGCTCCACAGTCTACATACCAGACTGAGCAAACTTATCAGCAGCCAGCCGCTCAGGAGCCGTTGTACCAACAGCCGCAACCCGTTGAACAGCAGCCTGTTGTGGAGCCTGAACCCGTT  >  NC_000913/934390‑934669

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