Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,832	T→C	12.7%	intergenic (‑71/+60)	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,832	0	T	C	12.7%	117.7 / 7.7	47	intergenic (‑71/+60)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base C (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 = 9.33e-01

GCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGC  >  NC_002947/1598684‑1598979                                                                                                                                                                  |                                                                                                                                                    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                                                                                                    <  3:195506/76‑1 (MQ=255)                                                                                                                                        gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                    >  4:195506/1‑76 (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                            >  6:141317/1‑135 (MQ=25)                                                                                                                                                     ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                            <  5:141317/135‑1 (MQ=25)                                                                                                                                                       tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg             >  4:128821/1‑148 (MQ=255)                                                                                                                                                                tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTg   >  2:34402/1‑149 (MQ=255)                                                                                                                                                                  tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                           <  6:36999/59‑1 (MQ=255)                                                                                                                                                                  tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                           >  5:36999/1‑59 (MQ=255)                                                                                                                                                                  tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc  >  6:185268/1‑148 (MQ=255)                                                                                                                                                                  |                                                                                                                                                    GCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGC  >  NC_002947/1598684‑1598979

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