Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,642,689	0	T	G	56.2%	5.4 / 13.8	16	A52A (GCA→GCC)	rodZ	MreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (9/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-02

GAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTC  >  NC_000913/2642539‑2642835                                                                                                                                                       |                                                                                                                                                     gAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGt                                                                                                                                                      <  2:263133/151‑1 (MQ=255)       aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCtt                                                                                                                                                 >  2:55783/1‑150 (MQ=255)                 aTCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTcttcttt                                                                                                                                        >  1:243849/1‑148 (MQ=255)                                                    cTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATCTTCTTTAATGTTGCGTACCGTGGAAAACTTCAGGCaaaaacg                                                                                                    >  2:332836/1‑146 (MQ=255)                                                       ttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTTTTTTTTTATTTCGCGGACCGCGGGAAACCTCCGGCaaaa                                                                                                       >  2:218339/1‑143 (MQ=255)                                                                tGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTATTTTCTTTCATTTTGCGTACCGTGGAAAACCTCCGGCAAAATCGCT‑C‑TGCAACGGcc                                                                                     >  1:123593/1‑151 (MQ=255)                                                                      cagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCGTATTTTTTTCTATGTGTCGTGCCCTGGGAACATTCCGGGAAAAACGCT‑C‑TGCCACCGcc                                                                                     >  1:393552/1‑145 (MQ=255)                                                                              cttctGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGGCTTTTAATTTTCAATTTCGCCTAACGTGGGAACAATTCGGCAAAGACGCT‑ccccaaacgcc                                                                                      >  1:370100/1‑127 (MQ=255)                                                                                                                gATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTATTTTCTTCAATGTCCGGCACCGAGGGAACCTCCCGGCAAAGACGCT‑CGGGCACCGGCC‑GCCGACCAAGTCCTAGTTTTTCGCGAGCATTTCGGAGGCGAggggcg                                     >  1:15139/1‑146 (MQ=255)                                                                                                                    tatCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTTTTTTTTTTAAATTCGCGGACCGTGGGAAACCTCAGGCAAAAACGCTCC‑GCCAAGGGCCGCCGGATAAAGTCCT‑GTTTTTTGCGGGCATTTCGCCGGGGGGGGCCggg                                   >  2:259464/1‑148 (MQ=255)                                                                                                                    tatCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGGGCCTTTTTTTTTTTCATTTCGCGTACCCCGGAAAACATTACGCAAAAACGCT‑C‑GGCAACGGCCTGCCGACTAAGTCCTAGTTTTTTTCGCGCATTATGCAGGCGGGCGCCggggg                                 >  2:203983/1‑148 (MQ=255)                                                                                                                         gcgcAGGGAAGTTGAAGCTAGATCGGCGGGTGCCGTATCTTCTTCGATGTTGCGTACCGTGGAAACCTGCAGGCAACGTCGCT‑C‑GGCAGCGGCCTGCTGAGTAAGTCCGAGTTGTTCGCGAGCATTACGCAGGGGAGCGCCGGTAGTAAg                             <  1:235762/150‑1 (MQ=255)                                                                                                                            cAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCt                         <  2:185418/151‑1 (MQ=255)                                                                                                                                   gTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCAttt                   <  2:120049/150‑1 (MQ=255)                                                                                                                                     tGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGGACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCATCGGCCTGATGGCTAAGACCTGGTTGTTCGCGGGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTGTATTTTg                 <  1:253223/150‑1 (MQ=255)                                                                                                                                      gAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTg                 <  2:221231/149‑1 (MQ=255)                                                                                                                                                   cggGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc  <  2:257284/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCT‑C‑GGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTC  >  NC_000913/2642539‑2642835

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