Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,171,892	0	C	G	43.3%	38.9 / 10.6	30	intergenic (+136/‑165)	rrfB/murB	5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (13/4);  new base G (0/13);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 2.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTT  >  NC_000913/4171743‑4172039                                                                                                                                                       |                                                                                                                                                    gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGAGGc                                                                                                                                                     >  2:566798/1‑150 (MQ=14)  aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCgg                                                                                                                                                   <  2:688963/151‑1 (MQ=31)       ccAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAgg                                                                                                                                              <  2:535162/151‑1 (MQ=31)              tCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGAcgccc                                                                                                                                        >  1:601575/1‑150 (MQ=31)                      aaCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                                >  1:338695/1‑150 (MQ=31)                            aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                                           >  1:650567/1‑149 (MQ=31)                             ggCTCAGTCGAAAGACTGGGCTTTTCTTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGAAAAAACCCGCGGGGGGCGGATTTGAAGGTTGCGAAGCAACGGCCG‑GGAGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAgg                                                                                                                        <  2:539780/151‑1 (MQ=11)                                tCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAt                                                                                                                     >  2:472171/1‑151 (MQ=31)                                                                                cGCTCCCCTGGTTAGAACAATCCCGCGGGGAGTGTTTTCAAAGTTTGCGAAGCACCGGCGG‑GGGGGGTGGGGGGCAGGGCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTa                                                                     <  2:212866/151‑1 (MQ=2)                                                                                    cgccTGAGGAGGACAAACCCGCGGGGGGCGGATTTGAACTTTGCGAAGCAACGGCCG‑GGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa                                                                   <  2:205211/147‑1 (MQ=11)                                                                                              ggAAAAACCCGCGGGGGGGGAATTGGAAGTTTGCAAAGCAACGGCCG‑GGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCtt                                                            <  2:775043/146‑1 (MQ=14)                                                                                              gaaaaaccccgcggGGGGGGGATTGGAGCGTTGCGAAGCAGCGGCGG‑GGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTg                                                        <  1:19666/142‑1 (MQ=14)                                                                                                      ccgccgGGGGGGTATTTGAACGTTGGGAAGCACGGGCGG‑GGGGGGGGGGGGGCGGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCatat                                               <  1:195384/151‑1 (MQ=17)                                                                                                       cgccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATc                                              >  2:20204/1‑151 (MQ=255)                                                                                                        gccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCt                                             >  1:509281/1‑151 (MQ=255)                                                                                                           gggggggattttGACCGTTGCGAAGCAACGGCCG‑GGGGGGGGGGGGGAGGGAGGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc                                           <  2:681109/141‑1 (MQ=17)                                                                                                                cGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGcc                                      >  2:365500/1‑150 (MQ=255)                                                                                                                 ggATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAt                                    >  2:747391/1‑151 (MQ=255)                                                                                                                      ttgactttagcaaGCCACCGCCCGGGGGGGGGGGGGGCGGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAt                                    <  1:273793/136‑1 (MQ=18)                                                                                                                      tGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccccc                               >  2:204589/1‑151 (MQ=255)                                                                                                                        aaCGTTGCGAAGCAGCGGCCG‑GGGGGGGGGGGGGCGGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCca                             <  2:214551/151‑1 (MQ=21)                                                                                                                          cGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac                            >  1:256348/1‑150 (MQ=255)                                                                                                                            ttgagaggcaggggCCG‑GGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAg                          <  2:755930/138‑1 (MQ=21)                                                                                                                             tGCAAAGCACCGGCCG‑GGGGGGTGGCGGGCAGCACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGa                         <  1:336542/150‑1 (MQ=255)                                                                                                                                gAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACg                     >  2:401112/1‑151 (MQ=255)                                                                                                                                 aaGCACGGGCCG‑GGGGGGGGGGGGGCAGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACg                     <  1:157005/150‑1 (MQ=21)                                                                                                                                 aaGCAACGGCCG‑GGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACg                     <  2:861233/150‑1 (MQ=255)                                                                                                                                  acgaggggCCG‑GGAGGGGGGGGGGGGGGGCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACgg                    <  1:692629/144‑1 (MQ=21)                                                                                                                                                 ggggggCGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGtt     <  2:453372/151‑1 (MQ=255)                                                                                                                                                    tGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTCAACTAGCCTCGttttt  >  1:657156/1‑151 (MQ=255)                                                                                                                                                    tGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGttttt  >  1:242182/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTT  >  NC_000913/4171743‑4172039

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