Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,273:1	+C	65.0%	intergenic (+464/+288)	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,273	1	.	C	65.0%	86.0 / 45.6	40	intergenic (+464/+288)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (9/5);  new base C (13/13);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741126‑4741405                                                                                                                                                      |                                                                                                                                          ggTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTga‑gca                                                                                                                                          <  1:6147/150‑3 (MQ=255)      gTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCCTCCGCCTCTTTGCTGGTGGTATTGTGGAGCAGCGGATCACCTTGTGTGACCGGCCATGGCGTATTTTTCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTga‑gcagcttc                                                                                                                                     <  2:1189744/150‑8 (MQ=255)                                               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                                                                                    >  1:484355/1‑91 (MQ=255)                                                                                                             aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                    <  2:484355/91‑1 (MQ=255)                                                                                                               tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                        >  2:1041744/1‑117 (MQ=255)                                                                                                               tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                        <  1:1041744/117‑1 (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)                                                                                                                                                      |                                                                                                                                          GGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741126‑4741405

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