Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,298:1	+T	60.9%	intergenic (+489/+263)	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,298	1	.	T	60.9%	95.4 / 61.6	46	intergenic (+489/+263)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (7/11);  new base T (16/12);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741153‑4741440                                                                                                                                                       |                                                                                                                                               gggTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGaatcc                                                                                                                                              >  2:148387/1‑147 (MQ=255)                   cctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGaca                                                                                                                           >  2:963894/1‑150 (MQ=255)                       ctTGATGGTCGTAATGTGCAGCACGGCATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                                                          <  2:182499/147‑1 (MQ=255)                       ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtc                                                                                                                       >  1:205311/1‑150 (MQ=255)                       ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GAGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac                                                                                                                          >  1:182499/1‑147 (MQ=255)                           aTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc                                                                                                                   >  1:939143/1‑150 (MQ=255)                               tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCa                                                                                                               <  2:864023/150‑1 (MQ=255)                                                     aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCAGCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                                                         >  2:597558/1‑150 (MQ=255)                                                     aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                                                         >  2:568052/1‑150 (MQ=255)                                                           cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                   <  1:701747/150‑1 (MQ=255)                                                            gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGcc                                                                                                          <  2:990596/126‑1 (MQ=255)                                                            gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGcc                                                                                                          >  1:990596/1‑126 (MQ=255)                                                             tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                                 <  2:939143/150‑1 (MQ=255)                                                              cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                                <  2:787408/150‑1 (MQ=255)                                                               aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                               >  1:559961/1‑150 (MQ=255)                                                               aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTGCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                               >  2:157928/1‑150 (MQ=255)                                                                ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                                                        <  1:141991/124‑1 (MQ=255)                                                                ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                              >  2:939419/1‑150 (MQ=255)                                                                ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                                                        >  2:141991/1‑124 (MQ=255)                                                                     cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                                                  >  1:241250/1‑125 (MQ=255)                                                                     cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                                                  <  2:241250/125‑1 (MQ=255)                                                                      aTGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                        <  1:568052/150‑1 (MQ=255)                                                                         gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                                     >  2:503655/1‑150 (MQ=255)                                                                              ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            >  1:526508/1‑138 (MQ=255)                                                                              ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            <  2:526508/138‑1 (MQ=255)                                                                                   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                             >  2:683097/1‑132 (MQ=255)                                                                                   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                             <  1:683097/132‑1 (MQ=255)                                                                                      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                                                        <  2:772298/150‑1 (MQ=255)                                                                                      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                                                        <  1:1068396/150‑1 (MQ=255)                                                                                      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                                                        >  1:448193/1‑150 (MQ=255)                                                                                        tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcc                                                                                                                      >  1:164578/1‑86 (MQ=255)                                                                                        tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcc                                                                                                                      <  2:164578/86‑1 (MQ=255)                                                                                           aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACt                                                   >  2:1070604/1‑150 (MQ=255)                                                                                            ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                                                  <  1:963894/150‑1 (MQ=255)                                                                                            ggTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                                                  <  2:82814/150‑1 (MQ=255)                                                                                             gTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                 <  1:567795/150‑1 (MQ=255)                                                                                                    tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg                                          <  2:205311/150‑1 (MQ=255)                                                                                                    tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg                                          >  1:1087533/1‑150 (MQ=255)                                                                                                    tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg                                          >  1:66283/1‑150 (MQ=255)                                                                                                        gTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                                                        >  1:240311/1‑84 (MQ=255)                                                                                                        gTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                                                        <  2:240311/84‑1 (MQ=255)                                                                                                        gTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACA‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            <  1:801230/112‑1 (MQ=255)                                                                                                        gTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                            >  2:801230/1‑112 (MQ=255)                                                                                                         tCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc                                     <  1:157928/150‑1 (MQ=255)                                                                                                                                              aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGa                              >  2:1020699/1‑123 (MQ=255)                                                                                                                                              aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGa                              <  1:1020699/123‑1 (MQ=255)                                                                                                                                               ttCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt  <  1:148387/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               GGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT  >  NC_002947/4741153‑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: ‑