Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,090,056	C→G	9.5%	S36C (TCC→TGC) ‡	gudP →	glucarate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,090,056	0	C	G	9.5%	145.4 / 7.5	53	S36C (TCC→TGC) ‡	gudP	glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (19/28);  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 = 1.95e-01

GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATC  >  NC_002947/4089918‑4090201                                                                                                                                           |                                                                                                                                                   gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcga                                                                                                                                          <  8:158747/149‑3 (MQ=255) gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTGcga                                                                                                                                          <  7:217541/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)                                                                                                                                           |                                                                                                                                                   GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGA‑CCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATC  >  NC_002947/4089918‑4090201

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