Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,301	A→C	35.7%	intergenic (+492/+260)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,301	0	A	C	35.7%	35.2 / 35.6	42	intergenic (+492/+260)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base C (8/7);  total (22/20)
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.74e-01

CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741172‑4741440                                                                                                                                         |                                                                                                                                            ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                                                          <  2:1312539/150‑1 (MQ=255)       gATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc                                                                                                                    <  2:1701087/150‑1 (MQ=255)                aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgcgc                                                                                                           >  1:1273859/1‑150 (MQ=255)                        gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTATGGTAACCATCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCATCCTGCAGCGCCTCGACCt                                                                                                   >  1:1258930/1‑150 (MQ=255)                         caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                                                  >  1:88472/1‑150 (MQ=255)                          aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                                 <  1:284298/150‑1 (MQ=255)                                    gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCATGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                       >  1:1720677/1‑150 (MQ=255)                                    gcggcgTCACCGGCCATGGCGTATTGAAAAGGGTATAAGTCAAGCGCTGAACACCGTCGTGCGGGATC‑AG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                       >  2:1677040/1‑150 (MQ=255)                                     cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                                                      <  1:1345472/150‑1 (MQ=255)                                        cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                   <  2:135185/150‑1 (MQ=255)                                        cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                   <  1:178722/150‑1 (MQ=255)                                          tCACCGGCCTTGGCGTATTGAACTGGGTATAGTTGCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                                 <  2:1633123/150‑1 (MQ=255)                                          tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                                 <  2:1080947/150‑1 (MQ=255)                                          tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                                 <  1:1154349/150‑1 (MQ=255)                                            aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                               >  1:568224/1‑150 (MQ=255)                                                gCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                                           <  1:172594/150‑1 (MQ=255)                                                    tGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                                       >  2:1854476/1‑150 (MQ=255)                                                           ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                                                                >  2:1698015/1‑150 (MQ=255)                                                           ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                                                                >  2:277398/1‑150 (MQ=255)                                                             gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                                                              <  1:1681363/150‑1 (MQ=255)                                                               aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                                                       <  2:484355/91‑1 (MQ=255)                                                               aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                                                       >  1:484355/1‑91 (MQ=255)                                                                 tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                           <  1:1041744/117‑1 (MQ=255)                                                                 tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                           >  2:1041744/1‑117 (MQ=255)                                                                   ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                                        >  1:304988/1‑118 (MQ=255)                                                                   ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                                        <  2:304988/118‑1 (MQ=255)                                                                    gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                                                       <  2:695506/150‑1 (MQ=255)                                                                      ataGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa                                                     >  2:1485546/1‑150 (MQ=255)                                                                         ggTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                                                  <  2:88472/150‑1 (MQ=255)                                                                          gTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                 >  1:296576/1‑150 (MQ=255)                                                                          gTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                 >  2:474073/1‑150 (MQ=255)                                                                           tCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg                                                >  2:1256241/1‑150 (MQ=255)                                                                           tCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg                                                >  1:283226/1‑150 (MQ=255)                                                                           tCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg                                                >  1:772618/1‑150 (MQ=255)                                                                            ccAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa                                               >  1:1333123/1‑150 (MQ=255)                                                                                      tCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt                                                      <  2:1402680/133‑1 (MQ=255)                                                                                      tCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt                                                      >  1:1402680/1‑133 (MQ=255)                                                                                      tCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                           <  2:1497504/144‑1 (MQ=255)                                                                                      tCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc                                     >  1:1079106/1‑150 (MQ=255)                                                                                      tCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                           >  1:1497504/1‑144 (MQ=255)                                                                                                                        cttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc        <  1:1677040/148‑1 (MQ=255)                                                                                                                              ttCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt  <  2:1720677/150‑1 (MQ=255)                                                                                                                                         |                                                                                                                                            CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741172‑4741440

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