Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,274	A→C	32.3%	intergenic (+465/+287)	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,274	0	A	C	32.3%	44.0 / 24.5	40	intergenic (+465/+287)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base C (8/5);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741172‑4741405                                                                                                         |                                                                                                                                         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                                                        >  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                                                     <  2:304988/118‑1 (MQ=255)                                                                   ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                     >  1:304988/1‑118 (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)                                                                                                         |                                                                                                                                         CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACC  >  NC_002947/4741172‑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: ‑