Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,276	C→G	15.0%	intergenic (+467/+285)	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,276	0	C	G	15.0%	145.1 / 14.4	52	intergenic (+467/+285)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base G (3/5);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCA  >  NC_002947/4741155‑4741408                                                                                                                            |                                                                                                                                          gTCAGGCTGCCATCCGCCTCCTGGACGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CC‑GCGACT‑TGCATTCGTGaatccag                                                                                                            <  2:1196155/150‑6 (MQ=255)  tCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGaatccagg                                                                                                           <  1:2199534/150‑7 (MQ=255)          ccATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACATTCGTGCATGAG‑CA‑GC‑TGCG‑GCGATT‑TGCATTCGTTAA‑TCCAGGCACAGGcc                                                                                                  >  2:260724/1‑150 (MQ=255)                  ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATGCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                          >  1:1451634/1‑150 (MQ=255)                   tcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACg                                                                                         >  2:296201/1‑150 (MQ=255)                    cctTGATGGTCGTAATGTGCAGCACGGAATCACGAGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGt                                                                                        >  1:332049/1‑150 (MQ=255)                     ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                       <  1:1961742/150‑1 (MQ=255)                          tGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTcctcct                                                                                  >  2:2415303/1‑150 (MQ=255)                              cGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAg                                                                              >  1:269999/1‑150 (MQ=255)                               gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc                                                                             >  1:331965/1‑150 (MQ=255)                               gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc                                                                             <  2:331965/150‑1 (MQ=255)                                 aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgcgc                                                                           >  1:592224/1‑150 (MQ=255)                                      gcagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      <  2:269999/150‑1 (MQ=255)                                         gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCGAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAcat                                                                   >  1:1094181/1‑148 (MQ=255)                                          caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACA‑CGTTGTCTGGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                  <  2:1165220/150‑1 (MQ=255)                                           aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                 <  1:786799/150‑1 (MQ=255)                                            cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                <  1:2416616/150‑1 (MQ=255)                                                   aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                         >  2:120772/1‑150 (MQ=255)                                                     gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                       >  2:175499/1‑150 (MQ=255)                                                       ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                     <  1:798856/150‑1 (MQ=255)                                                          gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTAACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                  >  1:118490/1‑150 (MQ=255)                                                            cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                <  1:1293308/150‑1 (MQ=255)                                                              ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                              >  2:2407899/1‑150 (MQ=255)                                                              ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                              >  1:1951835/1‑150 (MQ=255)                                                                ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                            <  1:2116771/150‑1 (MQ=255)                                                                       gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                     <  2:592224/150‑1 (MQ=255)                                                                          tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                  >  2:2388218/1‑150 (MQ=255)                                                                             tGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCGTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                               <  2:29412/150‑1 (MQ=255)                                                                               aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                             <  1:2415303/150‑1 (MQ=255)                                                                                 cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                           <  2:2247397/150‑1 (MQ=255)                                                                                  tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                          >  1:1865906/1‑150 (MQ=255)                                                                                   gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                         <  1:766812/150‑1 (MQ=255)                                                                                    ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                        <  2:2214123/150‑1 (MQ=255)                                                                                      tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                           <  1:1547492/113‑1 (MQ=255)                                                                                      tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                           >  2:1547492/1‑113 (MQ=255)                                                                                      tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt                      >  2:1222072/1‑150 (MQ=255)                                                                                        taGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTAc                    <  1:260724/150‑1 (MQ=255)                                                                                          ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  >  1:79963/1‑150 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                        <  2:1794721/95‑1 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                        >  1:1794721/1‑95 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTCTG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                 >  1:665188/1‑150 (MQ=255)                                                                                            tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                <  2:1109011/102‑1 (MQ=255)                                                                                            tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                >  1:1109011/1‑102 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  >  2:1782486/1‑146 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  <  1:1782486/146‑1 (MQ=255)                                                                                               aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc             >  1:275194/1‑150 (MQ=255)                                                                                               aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc             >  2:306901/1‑150 (MQ=255)                                                                                                 ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      >  2:2377254/1‑91 (MQ=255)                                                                                                 ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      <  1:2377254/91‑1 (MQ=255)                                                                                                 ctctGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg           <  2:332049/150‑1 (MQ=255)                                                                                                       tCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc     >  1:2131376/1‑150 (MQ=255)                                                                                                          ccTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCa  <  2:118490/150‑1 (MQ=255)                                                                                                                            |                                                                                                                                          GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCA  >  NC_002947/4741155‑4741408

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