Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,299	0	C	T	39.6%	38.1 / 50.2	53	intergenic (+490/+262)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (21/11);  new base T (8/13);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 8.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741166‑4741440                                                                                                                                           |                                                                                                                                               aTCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGcccc                                                                                                                                 >  1:953393/1‑150 (MQ=255)     gCCTCCTTGATGGTCGTAATGTGGAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac                                                                                                                             >  1:1466666/1‑150 (MQ=255)                 gTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGc                                                                                                                 <  2:1995979/150‑1 (MQ=255)                    gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc                                                                                                              <  1:303929/150‑1 (MQ=255)                    gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgc                                                                                                              <  2:2361992/150‑1 (MQ=255)                        tgtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGccc                                                                                                          >  2:925132/1‑149 (MQ=255)                                      tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGTCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                            <  2:365182/150‑1 (MQ=255)                                       cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                                                                           <  2:453509/150‑1 (MQ=255)                                         cgcgGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                                         >  1:2039666/1‑150 (MQ=255)                                          gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                        >  1:566690/1‑150 (MQ=255)                                             gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                                     <  2:2294784/150‑1 (MQ=255)                                             gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                                     <  2:2296754/150‑1 (MQ=255)                                              cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                    <  2:2071448/150‑1 (MQ=255)                                                tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                                                                           <  1:198566/141‑1 (MQ=255)                                                tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                                                                           >  2:198566/1‑141 (MQ=255)                                                  aCCGGCCATGGAGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTACTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                                >  2:618987/1‑150 (MQ=255)                                                   ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                      >  2:875895/1‑127 (MQ=255)                                                   ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                                                                                      <  1:875895/127‑1 (MQ=255)                                                   ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGaatc                                                                               >  1:626182/1‑147 (MQ=255)                                                     ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                             <  1:1485784/150‑1 (MQ=255)                                                       ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                    >  2:1161503/1‑141 (MQ=255)                                                       ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                    <  1:1161503/141‑1 (MQ=255)                                                       ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                           <  2:953393/150‑1 (MQ=255)                                                        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                                          >  2:686126/1‑150 (MQ=255)                                                                aTTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑cgcg                                                                 >  1:652594/1‑150 (MQ=255)                                                                 ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGa                                                                >  1:9883/1‑150 (MQ=255)                                                                 ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGa                                                                >  2:645115/1‑150 (MQ=255)                                                                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑AGCGATAGGc                                                           >  2:1968862/1‑149 (MQ=255)                                                                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑AGCGATAGGc                                                           <  1:1968862/149‑1 (MQ=255)                                                                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGcc                                                          >  2:821322/1‑150 (MQ=255)                                                                         ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCtt                                                        >  2:1514661/1‑150 (MQ=255)                                                                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                             >  1:1305310/1‑127 (MQ=255)                                                                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                             <  2:1305310/127‑1 (MQ=255)                                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACGTTGGCCTGCGTCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTAc                                                    <  1:686126/150‑1 (MQ=255)                                                                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAg                                                             >  1:653302/1‑139 (MQ=255)                                                                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACtt                                                  <  2:566690/150‑1 (MQ=255)                                                                               ggTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAg                                                             <  2:653302/139‑1 (MQ=255)                                                                                gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTg                                                 >  2:1339773/1‑150 (MQ=255)                                                                                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTgg                                                >  2:619118/1‑150 (MQ=255)                                                                                   cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGAc                                              >  1:852377/1‑150 (MQ=255)                                                                                    aGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGAcc                                             >  1:201337/1‑150 (MQ=255)                                                                                     gCTCTGGTCACCTTCGTACTTGAG‑CA‑GCTTG‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGAccc                                            <  1:766280/150‑1 (MQ=255)                                                                                       tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCgg                                          >  1:482137/1‑150 (MQ=255)                                                                                       tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCgg                                          <  2:482137/150‑1 (MQ=255)                                                                                        ctGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCATTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGc                                         <  1:2197836/150‑1 (MQ=255)                                                                                                                       tgcgacttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcg             >  1:941877/8‑150 (MQ=255)                                                                                                                       tgcgacttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTCGTCACCCACCGGCTCACCCTTGG‑CGGATAGTTGGACCCGGCCACCTCACTGAGGGGCACGTTGATGCcg             >  1:295260/8‑150 (MQ=255)                                                                                                                            cttggaTTCGTCAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc        <  2:1396263/145‑1 (MQ=255)                                                                                                                             ttGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcg       >  2:1067746/2‑150 (MQ=255)                                                                                                                                 aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcct   >  1:1629857/1‑150 (MQ=255)                                                                                                                                  ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc        >  2:2158271/1‑144 (MQ=255)                                                                                                                                  ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc        <  1:2158271/144‑1 (MQ=255)                                                                                                                                  ttCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt  <  2:9883/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                               ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAC‑CGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741166‑4741440

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: ‑