Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,528,340	C→T	22.3%	G213G (GGC→GGT)	ftsA →	cell division protein FtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,528,340	0	C	T	22.3%	22.9 / 5.3	18	G213G (GGC→GGT)	ftsA	cell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (2/2);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

GCAGTCAACGCCGCCCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATC  >  NC_002947/1528191‑1528478                                                                                                                                                       |                                                                                                                                           gcAGTCAACGCCGCCCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑ggc                                                                                                                                            <  1:41709/150‑1 (MQ=255)          gccgccCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCaccacc                                                                                                                                   >  1:790142/1‑150 (MQ=255)            cgccCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGa                                                                                                                                 <  1:567039/150‑1 (MQ=255)             ggccAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat                                                                                                                                <  2:219664/148‑2 (MQ=255)                cAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATc                                                                                                                             >  2:621596/1‑150 (MQ=255)                                                          gatgacATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATc                                                                                   <  2:759481/147‑1 (MQ=255)                                                                atcatcCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACACCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGAtt                                                                             <  2:790142/150‑1 (MQ=255)                                                                        cGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGa                                                                     <  1:46095/150‑1 (MQ=255)                                                                                    ttCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCaa                                                         >  1:368765/1‑150 (MQ=255)                                                                                       ggCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat                                                                                                                                <  1:632295/76‑2 (MQ=255)                                                                                       ggCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat                                                                                                                                >  2:632295/1‑75 (MQ=255)                                                                                                                 aaaaaGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCACGGAAATTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGa                            >  1:250726/1‑150 (MQ=255)                                                                                                                 aaaaaGAGCTGGGCGTGTGCCTGGTGGACATCCGTAGCCGGCA‑GTTCGAGCTGACCGGCTTCACCGAGGGGTGCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTTACGCACCACCTCGCCATGGCCCTGCGTACCCCGa                            >  1:337503/1‑150 (MQ=255)                                                                                                                  aaaaGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCTACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGAATGCAACGACATCGCCATGGCCCTGCGTACCCCGAc                           <  1:585651/150‑1 (MQ=255)                                                                                                                   aaaGAGCTGGGCGTTTGCCTGTTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGAcc                          <  2:337503/150‑1 (MQ=255)                                                                                                                         cTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTAc                    >  1:730802/1‑150 (MQ=255)                                                                                                                          tGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACg                   >  1:272352/1‑150 (MQ=255)                                                                                                                                   gCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGaaa          <  2:730802/150‑1 (MQ=255)                                                                                                                                           gatgtCGGT‑GGTGCCGCCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATc  >  1:298754/5‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                           GCAGTCAACGCCGCCCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATC  >  NC_002947/1528191‑1528478

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