Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	544,106	0	G	A	37.5%	9.2 / 12.5	16	K781K (AAG→AAA)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base A (3/3);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATCCGTCAGCTGGCCGGTATGCGTGGTCTGATGGCCAAGCCGGACGGTTCCATCATCGAGACGCCGATCACCGCGAACTTCCGTGAAGGTCTGAGCGTACTTCAGTACTTCATCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTGGTGATGACGC  >  NC_002947/543971‑544233                                                                                                                                        |                                                                                                                                gATCCGTCAGCTGGCCGGTATGCGTGGTCTGATGGCCAAGCCGGACGGTTCCATCATCGAGACGCCGATCACCGCGAACTTCCGTGAAGGTCTGAGCGTACTTCAGTACTTCATCTCGACTCACGGTGCTCGTAAg                                                                                                                                 <  2:633000/136‑1 (MQ=255)             tgctggTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCCCGGTGCTCGTAAAGGTCTGGCGGAt                                                                                                                     >  1:598951/5‑136 (MQ=255)             tgctggTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGAt                                                                                                                     >  1:292590/5‑136 (MQ=255)             tgctggTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGAt                                                                                                                     >  1:482964/5‑136 (MQ=255)                                                  tCCATCATCGAGACGCCGATCACCGCGAACTTCCGTGAAGGTCTGAGCGTACTTCAGTACTTCATCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACtc                                                                                <  1:381946/136‑1 (MQ=255)                                                                          gcgAACTTCCGTGAAGGTCTGAGCGTACTTCAGTACTTCATCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAgg                                                        <  1:290475/136‑1 (MQ=255)                                                                                                       tcagtacttcaTCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCAGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGc                           >  2:112774/1‑136 (MQ=255)                                                                                                            acttcaTCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGa                      <  1:243237/136‑1 (MQ=255)                                                                                                                caTCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCaa                  >  1:506881/1‑135 (MQ=255)                                                                                                                   ctcGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGa                                  <  1:401962/117‑1 (MQ=255)                                                                                                                   ctcGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGa                                  >  2:401962/1‑117 (MQ=255)                                                                                                                        cTCACGCTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTggg          >  2:235479/1‑135 (MQ=255)                                                                                                                                gCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACtc                                                                                >  1:146803/1‑58 (MQ=255)                                                                                                                                gCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACtc                                                                                <  2:146803/58‑1 (MQ=255)                                                                                                                                gCTCGTAAAGGTCTGGCGTATACCGCACTGAATACTGCGAACTCCTGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTGGTGATGACGc  <  2:292590/136‑1 (MQ=255)                                                                                                                                gCTCGTAAAGGTCTGGCGGATACCGCACTGAAATCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTGGTGATGACGc  <  2:598951/136‑1 (MQ=255)                                                                                                                                gCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTGGTGATGACGc  <  2:482964/136‑1 (MQ=255)                                                                                                                                        |                                                                                                                                GATCCGTCAGCTGGCCGGTATGCGTGGTCTGATGGCCAAGCCGGACGGTTCCATCATCGAGACGCCGATCACCGCGAACTTCCGTGAAGGTCTGAGCGTACTTCAGTACTTCATCTCGACTCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACCGCTAACTCCGGTTACCTGACTCGTCGTCTGGTGGACGTTGCCCAGGATCTGGTAGTGACCGAGATCGACTGCGGCACCGATCAGGGCCTGGTGATGACGC  >  NC_002947/543971‑544233

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