Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,090,057	C→G	9.2%	S36S (TCC→TCG) ‡	gudP →	glucarate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,090,057	0	C	G	9.2%	145.5 / 5.4	53	S36S (TCC→TCG) ‡	gudP	glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base G (2/3);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCA  >  NC_002947/4089918‑4090205                                                                                                                                            |                                                                                                                                                      gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcga                                                                                                                                              <  7:217541/149‑3 (MQ=255) gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcga                                                                                                                                              <  8:158747/149‑3 (MQ=255) gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTg                                                                                                                                                 >  8:104837/1‑146 (MQ=255) gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTg                                                                                                                                                 >  4:255370/1‑146 (MQ=255)   acaAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcgata                                                                                                                                            <  6:244386/149‑5 (MQ=255)    caAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaaa                                                                                                                                           <  7:18336/149‑1 (MQ=255)    caAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaaa                                                                                                                                           <  2:80359/149‑1 (MQ=255)    caAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGa                                                                                                                                             >  1:185894/1‑147 (MQ=255)           aaCTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCg                                                                                                                                   <  7:330014/149‑1 (MQ=255)              tAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCa                                                                                                                                >  7:266439/1‑149 (MQ=255)                gTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCAtt                                                                                                                              >  5:109821/1‑149 (MQ=255)                   ggAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGAc                                                                                                                           >  2:23787/1‑149 (MQ=255)                      gTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGcc                                                                                                                        <  7:104837/149‑1 (MQ=255)                      gTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGcc                                                                                                                        <  7:131888/149‑1 (MQ=255)                         ttGCCAGCAGGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGt                                                                                                                      <  6:292332/148‑1 (MQ=255)                                 aTGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCgg                                                                                                             <  1:23787/149‑1 (MQ=255)                                 aTGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCgg                                                                                                             <  4:238341/149‑1 (MQ=255)                                                       aTGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCgg                                                                                        <  1:320096/149‑1 (MQ=255)                                                          tgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATg                                                                                     <  1:289876/148‑1 (MQ=255)                                                                  cTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCg                                                                                         <  3:34975/136‑1 (MQ=255)                                                                  cTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCg                                                                                         >  4:34975/1‑136 (MQ=255)                                                                     aTCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggc                                                                         <  3:255370/149‑1 (MQ=255)                                                                      tCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggct                                                                        >  6:281601/1‑149 (MQ=255)                                                                      tCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggct                                                                        >  2:60335/1‑149 (MQ=255)                                                                        cTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctg                                                                       <  2:185894/148‑1 (MQ=255)                                                                          gTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCAAGACTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                     >  5:306773/1‑148 (MQ=255)                                                                              aTGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATgg                                                                                    <  8:237799/129‑1 (MQ=255)                                                                                  tGTTTCTGGTGACCACGATAAAATATGCCGACCGCGCGACCATTTCTATCGCAGGCTACAGCCTGCAGAAAGA‑CATCGGCATTGACGACGTAACACTCGGTTATATCTTTTCCTCCTTCGGATGGGCTTACGTGTCTCGCCAGATTCCt                                                            >  7:42673/1‑149 (MQ=255)                                                                                    tttCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCt                                                            <  6:109821/147‑1 (MQ=255)                                                                                          gTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggct                                                     <  8:230906/148‑1 (MQ=255)                                                                                          gTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCAAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCtggtgg                                                       <  3:158642/146‑1 (MQ=255)                                                                                             aCCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTagctggct                                                 <  2:27657/149‑1 (MQ=255)                                                                                              ccACTATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGTGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTa                                                <  8:42673/149‑1 (MQ=255)                                                                                                    tCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGAt                                          <  7:292994/149‑1 (MQ=255)                                                                                                       aCTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGc                                       >  4:82946/1‑149 (MQ=255)                                                                                                       aCTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGc                                       >  6:178686/1‑149 (MQ=255)                                                                                                           tGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTc                                    <  5:281601/148‑1 (MQ=255)                                                                                                           tGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCg                                   >  6:120570/1‑149 (MQ=255)                                                                                                            gCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCgg                                  >  1:314736/1‑149 (MQ=255)                                                                                                                      cgACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAAc                        <  6:306773/149‑1 (MQ=255)                                                                                                                        aCCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACt                                              >  6:410914/1‑125 (MQ=255)                                                                                                                        aCCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACt                                              <  5:410914/125‑1 (MQ=255)                                                                                                                           aTTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTa                   >  3:50093/1‑149 (MQ=255)                                                                                                                           aTTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCGGGCTACTCGATGGCTTCGGGTCCATGAACGTCTa                   >  2:129627/1‑149 (MQ=255)                                                                                                                            tttcTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGTATGGGCTTACGTGGCTGGCCTGATTCCTGGTGGCTGGCTACTCGTTCGCTTCGGtt                                <  2:296376/132‑1 (MQ=255)                                                                                                                            ttGCTATCGCAGGCTCCAGCCTTCAGAAAGA‑CCTCGGCAATGCCGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGAGGCTGGCTACTCGATAGCATCGCtt                                >  1:296376/1‑135 (MQ=255)                                                                                                                              ctttctgCAGGCTGGAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGc                >  4:47163/7‑149 (MQ=255)                                                                                                                               cTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGc                <  3:82946/148‑1 (MQ=255)                                                                                                                                  tCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCt              <  5:120570/147‑1 (MQ=255)                                                                                                                                  tCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTc            <  2:314736/149‑1 (MQ=255)                                                                                                                                  tCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCAAAGAACGTCTACGCCtt             <  1:60335/148‑1 (MQ=255)                                                                                                                                     cAGGCTGGAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTAAGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAg          >  8:58250/1‑148 (MQ=255)                                                                                                                                        gCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTATTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATc      >  5:402589/1‑149 (MQ=255)                                                                                                                                            cAGCCTGCAGAAAGA‑CCTCGCCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGACTACGGATGGGATTACGTGGCTGGCAAGATTCATGGTCGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACTCCTACAGCATcataa  >  1:162112/1‑145 (MQ=255)                                                                                                                                            |                                                                                                                                                      GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCA  >  NC_002947/4089918‑4090205

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