Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,298	1	.	T	43.9%	60.2 / 84.6	41	intergenic (+489/+263)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (15/8);  new base T (8/10);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741163‑4741440                                                                                                                                              |                                                                                                                                               gCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGATCTTTTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGc                                                                                                                                   >  1:1051814/1‑150 (MQ=255)       cGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGa                                                                                                                             >  1:491204/1‑150 (MQ=255)           tcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACg                                                                                                                         >  1:512475/1‑150 (MQ=255)                    gTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGc                                                                                                                >  1:704125/1‑150 (MQ=255)                     tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCa                                                                                                               >  2:1444164/1‑150 (MQ=255)                        tAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgcg                                                                                                            >  1:226866/1‑150 (MQ=255)                                   aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAAGTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGTCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                                 >  2:985474/1‑150 (MQ=255)                                   aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                                 <  1:1450573/150‑1 (MQ=255)                                    cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGTTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                                                <  1:3519/150‑1 (MQ=255)                                        aTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAg                                                                                            <  1:1412465/150‑1 (MQ=255)                                         tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                           <  2:512475/150‑1 (MQ=255)                                           aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                                                         >  1:371212/1‑150 (MQ=255)                                             gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                       >  2:959950/1‑150 (MQ=255)                                                   tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                                 >  1:1397940/1‑150 (MQ=255)                                                     aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                     >  2:1133511/1‑128 (MQ=255)                                                     aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                     <  1:1133511/128‑1 (MQ=255)                                                      ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                              >  1:1197011/1‑150 (MQ=255)                                                       cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                             <  2:1397940/150‑1 (MQ=255)                                                       cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCT‑GCGACT‑TGCATAGGGGAA‑TCCAGGCACAGGCCCCGGCCACGTGCTCCTGCAGCGCCTCGACCTTGGCAAGCGGCGCCTTGATCAc                                                                             >  1:543852/1‑150 (MQ=255)                                                        ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            <  2:1051814/150‑1 (MQ=255)                                                        ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            <  1:54567/150‑1 (MQ=255)                                                          ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGTTGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                          <  2:376882/150‑1 (MQ=255)                                                            aTGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                        <  2:371212/150‑1 (MQ=255)                                                                    ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                                                                >  2:913295/1‑150 (MQ=255)                                                                         cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTTAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                                   >  2:1054084/1‑142 (MQ=255)                                                                         cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTTAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                                   <  1:1054084/142‑1 (MQ=255)                                                                          tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                                          <  1:913295/150‑1 (MQ=255)                                                                               ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa                                                     >  1:210208/1‑150 (MQ=255)                                                                                  ggTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑GGCGACTTGCATTCGTGAA‑TCCAGGAATACGACAGCTACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                                                  <  1:295105/150‑1 (MQ=255)                                                                                   gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                 <  2:491204/150‑1 (MQ=255)                                                                                        gCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                     <  1:46271/93‑1 (MQ=255)                                                                                        gCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                     >  2:46271/1‑93 (MQ=255)                                                                                         ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                           >  2:1292441/1‑150 (MQ=255)                                                                                         ctctGGTCACCTTCGTGCTTGAG‑CA‑GCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                           >  1:765187/1‑150 (MQ=255)                                                                                                cACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATGCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCt                                    >  2:1236690/1‑150 (MQ=255)                                                                                                                               acttGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCGTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct         >  2:1611680/4‑150 (MQ=255)                                                                                                                                    cATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcc    <  1:985474/150‑1 (MQ=255)                                                                                                                                     aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACAcc                      <  2:1019661/131‑1 (MQ=255)                                                                                                                                     aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACAcc                      >  1:1019661/1‑131 (MQ=255)                                                                                                                                     aTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcct   >  1:117727/1‑150 (MQ=255)                                                                                                                                      ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt  <  2:1197011/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               GCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741163‑4741440

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