Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,239	0	A	G	40.4%	39.0 / 37.9	52	intergenic (+476/+277)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (22/9);  new base G (10/11);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCT  >  NC_002947/4741169‑4741340                                                                         |                                                                                                        cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGcgcct                                                  >  1:4023240/1‑124 (MQ=255)           ggCGTATGGAACTGGGTATGGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGcgcctcga                                               <  2:4023120/118‑5 (MQ=255)             cGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCa                                      >  1:4023363/1‑125 (MQ=255)                  tGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                 <  2:4023126/125‑1 (MQ=255)                   gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                >  2:4023369/1‑125 (MQ=255)                   gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                >  1:4023302/1‑125 (MQ=255)                     aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                              >  1:4023183/1‑125 (MQ=255)                     aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                              >  2:4023237/1‑125 (MQ=255)                      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                             >  2:4023474/1‑125 (MQ=255)                      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                             <  1:4023128/125‑1 (MQ=255)                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                            >  1:4023396/1‑125 (MQ=255)                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                            >  2:4023307/1‑125 (MQ=255)                        gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                           >  1:4023359/1‑125 (MQ=255)                        gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                           >  1:4023472/1‑125 (MQ=255)                         ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                          >  1:4023352/1‑125 (MQ=255)                           tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                        >  2:4023304/1‑125 (MQ=255)                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAATCCAAGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                       >  1:4023344/1‑125 (MQ=255)                            ataGGTCCAGCTCTGGTCAACTTCGTGCTTGAA‑CA‑GCTTCTGCGACTGGCATTCGTGAA‑TCCAGGCACAGGACCCGGACACGTACTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCATTGATCa                       <  2:4023130/125‑1 (MQ=255)                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCTACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                      >  1:4023413/1‑125 (MQ=255)                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                      >  1:4023338/1‑125 (MQ=255)                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                      <  2:4023131/125‑1 (MQ=255)                              aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                     >  1:4023322/1‑125 (MQ=255)                              aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                     <  2:4023132/125‑1 (MQ=255)                              aGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                     <  2:4023133/125‑1 (MQ=255)                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                    >  2:4023238/1‑125 (MQ=255)                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                    <  1:4023134/125‑1 (MQ=255)                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                    <  1:4023135/125‑1 (MQ=255)                                gTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                   >  2:4023388/1‑125 (MQ=255)                                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                  >  2:4023255/1‑125 (MQ=255)                                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                  >  1:4023244/1‑125 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                 <  1:4023147/125‑1 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                 <  2:4023193/125‑1 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                 <  2:4023138/125‑1 (MQ=255)                                   cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                >  1:4023347/1‑125 (MQ=255)                                    aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca               >  1:4023418/1‑125 (MQ=255)                                    aGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca               <  2:4023198/125‑1 (MQ=255)                                     gCTCTGGTCACCTTCGTGCTTGGG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              <  2:4023199/125‑1 (MQ=255)                                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              <  1:4023169/125‑1 (MQ=255)                                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              >  2:4023471/1‑125 (MQ=255)                                      ctctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc             <  1:4023170/125‑1 (MQ=255)                                      ctctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc             >  1:4023329/1‑125 (MQ=255)                                       tctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg            >  1:4023502/1‑125 (MQ=255)                                         tGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg          >  2:4023427/1‑125 (MQ=255)                                          ggTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa         <  2:4023188/125‑1 (MQ=255)                                          ggTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa         <  1:4023154/125‑1 (MQ=255)                                           gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt        <  1:4023151/125‑1 (MQ=255)                                           gTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt        >  1:4023362/1‑125 (MQ=255)                                            tCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa       <  2:4023194/125‑1 (MQ=255)                                              aCCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg     >  2:4023250/1‑125 (MQ=255)                                               ccTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc    >  2:4023475/1‑125 (MQ=255)                                                cTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc   >  1:4023847/1‑125 (MQ=255)                                                 ttCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt  >  2:4023424/1‑125 (MQ=255)                                                                         |                                                                                                        CACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCT  >  NC_002947/4741169‑4741340

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