Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,831	C→G	12.8%	intergenic (‑70/+61)	dctB ← / ← bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,598,831	0	C	G	12.8%	121.9 / 10.9	47	intergenic (‑70/+61)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base G (4/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

AGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTG  >  NC_002947/1598683‑1598978                                                                                                                                                                  |                                                                                                                                                    aGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActc                                                                                                                                                     <  1:168410/149‑1 (MQ=255)  gCCGGCGGTGAGTGCTGTAACCAGCACGATCAGCAGGGTGATCAGGCCCAGAGGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActct                                                                                                                                                    <  6:271130/149‑1 (MQ=255)    cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                  <  1:124818/149‑1 (MQ=255)          tGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAcc                                                                                                                                            >  7:107041/1‑149 (MQ=255)                     cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTTGGAGAAGGGGATGGGGAGAAGTCGGGCCCTGCCTCTTCGCGGGTAAACCCGCGCATACAGGAACTCTccagca                                                                                                                                              >  1:271274/1‑149 (MQ=255)                        aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGttt                                                                                                                              >  2:168158/1‑149 (MQ=255)                                 agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcgcg                                                                                                                      <  3:146882/148‑1 (MQ=255)                                 agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGCAGGCGCGGTTTTACCcgcg                                                                                                                      >  4:13038/1‑148 (MQ=255)                                   cagGGTGATCAGCCCCAGACGCAGTGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                       <  8:148282/129‑1 (MQ=37)                                   cagGGTGATCAGCCCCAGACGCAGTGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                       >  7:148282/1‑129 (MQ=37)                                          aTCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                             <  8:107041/148‑1 (MQ=255)                                               ccccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCAGatcgg                                                                                                       <  4:286715/149‑5 (MQ=255)                                               ccccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTACGCGGGTAAACCCGCGCCTACAGGTAATCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAg                                                                                                       >  5:171561/1‑149 (MQ=255)                                                 ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                     >  7:6715/1‑149 (MQ=255)                                                     aCGCAGGGCGCGAAAGGAAATGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGt                                                                                                                                         <  8:307577/109‑1 (MQ=255)                                                     aCGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCAGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGt                                                                                                                                         >  7:307577/1‑109 (MQ=255)                                                          gggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                            <  1:324549/149‑1 (MQ=255)                                                           ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCCCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                     >  5:8504/1‑139 (MQ=255)                                                           ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                     <  6:8504/139‑1 (MQ=255)                                                            gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGCAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                          <  3:13038/149‑1 (MQ=255)                                                            gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGGCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                          <  7:109050/149‑1 (MQ=255)                                                                gAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGCCCTTTTTCAGTc                                                                                       >  7:39632/1‑148 (MQ=255)                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                                        >  3:147000/1‑146 (MQ=255)                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                                        <  4:147000/146‑1 (MQ=255)                                                                      aaaaGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCaa                                                                                >  6:28893/1‑149 (MQ=255)                                                                          gggCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCa                                                                            <  3:96180/149‑1 (MQ=255)                                                                            gCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAATCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                          <  2:84425/149‑1 (MQ=255)                                                                                 gggggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGctc                                                                      <  8:6715/148‑1 (MQ=255)                                                                                    ggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCCCCTGTAGGCGCGGTTTAACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                    >  4:217759/1‑131 (MQ=255)                                                                                    ggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTAACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                    <  3:217759/131‑1 (MQ=255)                                                                                                          g‑tgg‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTg                                                     <  1:168158/149‑1 (MQ=255)                                                                                                              g‑‑‑gATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTggg                                                   <  5:28893/148‑1 (MQ=255)                                                                                                                        ggggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                           >  8:15855/1‑149 (MQ=255)                                                                                                                         gggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGcc                                          <  7:15855/149‑1 (MQ=255)                                                                                                                          ggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTGTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGccc                                         >  1:205139/1‑149 (MQ=255)                                                                                                                             ctgcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCGTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAc                                      >  7:174771/1‑149 (MQ=255)                                                                                                                                cctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTg                                   >  6:1577/1‑149 (MQ=255)                                                                                                                                 ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                           >  4:299755/1‑140 (MQ=255)                                                                                                                                 ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                           <  3:299755/140‑1 (MQ=255)                                                                                                                                       gctggTAAAACCGCGCCTACAGGTACTCTGAAGCAACTGTATGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTTCTGCTGCACCGCCTTGGAATCCAGCCCTACCTGGACATTg                            <  2:271274/146‑1 (MQ=255)                                                                                                                                         gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                   >  4:195506/1‑76 (MQ=255)                                                                                                                                         gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                   <  3:195506/76‑1 (MQ=255)                                                                                                                                             aaaaCCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGATCTGCTGCTGCACCGCCTgg                                                    >  3:52106/1‑119 (MQ=18)                                                                                                                                              aaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                     >  5:326961/1‑149 (MQ=255)                                                                                                                                              aaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGCCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                     <  8:39632/149‑1 (MQ=255)                                                                                                                                                   gcgcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGt                >  4:316359/1‑149 (MQ=37)                                                                                                                                                      ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                           <  5:141317/135‑1 (MQ=25)                                                                                                                                                      ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                           >  6:141317/1‑135 (MQ=25)                                                                                                                                                        tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg            >  4:128821/1‑148 (MQ=255)                                                                                                                                                                 tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTg  >  2:34402/1‑149 (MQ=255)                                                                                                                                                                  |                                                                                                                                                    AGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTG  >  NC_002947/1598683‑1598978

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