Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,273	1	.	C	73.3%	71.2 / 26.1	30	intergenic (+464/+288)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base C (12/10);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741129‑4741405                                                                                                                                                  |                                                                                                                                          aaGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTTAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagct                                                                                                                                       >  1:194792/1‑145 (MQ=255)  aGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagctt                                                                                                                                      >  1:201042/1‑144 (MQ=255)                    gTTACGGGTCAGGCTGCCATCCGCCTCCTTGCTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCT‑GCGACT‑TGCATTCGTg                                                                                                                >  1:841331/1‑150 (MQ=255)                                          gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac                                                                                         <  1:4356/150‑1 (MQ=255)                                            ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                       >  1:1143653/1‑150 (MQ=255)                                                                   gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCt                                                                >  1:814817/1‑150 (MQ=255)                                                                     aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                              <  1:288279/150‑1 (MQ=255)                                                                     aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                              >  2:1256829/1‑150 (MQ=255)                                                                           tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGt                                                                                     <  1:1142808/121‑1 (MQ=255)                                                                           tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGt                                                                                     >  2:1142808/1‑121 (MQ=255)                                                                           tCACGCGGCGTCACCGGACATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                        >  1:1305858/1‑150 (MQ=255)                                                                               gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                       >  2:80580/1‑115 (MQ=255)                                                                               gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                       <  1:80580/115‑1 (MQ=255)                                                                                        ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                           >  2:1227388/1‑150 (MQ=255)                                                                                             cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac                                                                                         >  1:1144394/1‑99 (MQ=255)                                                                                             cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGac                                                                                         <  2:1144394/99‑1 (MQ=255)                                                                                                    tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                               >  2:1210089/1‑150 (MQ=255)                                                                                                      ttGACCTGGGTATAGGTCCAGCTCTCGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGAAACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                             <  2:194792/150‑1 (MQ=255)                                                                                                      ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                             >  1:1121690/1‑150 (MQ=255)                                                                                                        gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                           <  2:1121690/150‑1 (MQ=255)                                                                                                         aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                          <  2:586104/150‑1 (MQ=255)                                                                                                           cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                        <  1:1210089/150‑1 (MQ=255)                                                                                                           cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac                                                                                       >  1:462130/1‑87 (MQ=255)                                                                                                           cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac                                                                                       <  2:462130/87‑1 (MQ=255)                                                                                                            tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                       >  2:306835/1‑150 (MQ=255)                                                                                                                 ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa                  >  1:400982/1‑150 (MQ=255)                                                                                                                     gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg              >  1:1267895/1‑150 (MQ=255)                                                                                                                      tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg             <  2:964917/150‑1 (MQ=255)                                                                                                                          gCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc         <  2:734782/150‑1 (MQ=255)                                                                                                                                 tCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc  >  1:262180/1‑150 (MQ=255)                                                                                                                                  cACCTTCGTGCTTGAGCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                           <  1:338550/108‑1 (MQ=255)                                                                                                                                  cACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                           >  2:338550/1‑108 (MQ=255)                                                                                                                                                  |                                                                                                                                          AAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741129‑4741405

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