Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,828	A→G	12.2%	intergenic (‑67/+64)	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,828	0	A	G	12.2%	118.1 / 8.3	49	intergenic (‑67/+64)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base A (20/23);  new base G (4/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

CCCAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGG  >  NC_002947/1598680‑1598968                                                                                                                                                                  |                                                                                                                                             cccAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTa                                                                                                                                              >  4:96180/1‑149 (MQ=255)    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                                                                              <  4:147000/146‑1 (MQ=255)                                                                    aaaGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                              >  3:147000/1‑146 (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                 >  6:141317/1‑135 (MQ=25)                                                                                                                                                         ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                 <  5:141317/135‑1 (MQ=25)                                                                                                                                                           tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg  >  4:128821/1‑148 (MQ=255)                                                                                                                                                                  |                                                                                                                                             CCCAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGG  >  NC_002947/1598680‑1598968

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