Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,090,057	C→G	11.0%	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	11.0%	149.4 / 8.0	53	S36S (TCC→TCG) ‡	gudP	glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (2/4);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCA  >  NC_002947/4089918‑4090205                                                                                                                                              |                                                                                                                                                     gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTGGAGCCTGcga                                                                                                                                             <  4:136584/149‑3 (MQ=255) gCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTGGAGCCTGcga                                                                                                                                             <  1:220615/149‑3 (MQ=255)    caAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTGGAGCCTGcgatag                                                                                                                                          <  8:41755/149‑6 (MQ=255)     aaaaaTAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTGGAGCCTGcgata                                                                                                                                           >  4:104786/1‑143 (MQ=255)     aaaaaTAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTGGAGCCTGcgata                                                                                                                                           <  3:104786/147‑5 (MQ=255)         aTAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCAGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCt                                                                                                                                     <  4:253206/149‑1 (MQ=255)            ccTAGTGGGAGTATTGCACTCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCgg                                                                                                                                  <  4:243053/148‑1 (MQ=255)            aCTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCg                                                                                                                                   <  5:127657/148‑1 (MQ=255)              tAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAACGACCTCGGCa                                                                                                                                >  7:267652/1‑149 (MQ=255)                gTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCAtt                                                                                                                              <  7:6678/149‑1 (MQ=255)                      gTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATCGACGcc                                                                                                                        >  2:290532/1‑149 (MQ=255)                            ccAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGc                                                                                                                                               >  2:63650/1‑120 (MQ=255)                            ccAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGc                                                                                                                                               <  1:63650/120‑1 (MQ=255)                              aGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCt                                                                                                                >  1:277068/1‑149 (MQ=255)                               gCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTc                                                                                                               >  7:51795/1‑149 (MQ=255)                                 aTGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCgg                                                                                                             <  3:103170/149‑1 (MQ=255)                                    cAAGCGACTAAGAAAACGCAGGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGCTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGtt                                                                                                           <  6:59087/148‑1 (MQ=255)                                      aGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTtata                                                                                                        <  5:207652/149‑1 (MQ=255)                                       gCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTtata                                                                                                        >  5:285098/1‑148 (MQ=255)                                            tAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACAACGATCAACAATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCttat                                                                                                  >  1:209867/1‑147 (MQ=255)                                                aaaaCGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGc                                                                                              >  2:34704/1‑149 (MQ=255)                                                   aCGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATCGACGCCGTAACCCTCGGTTATATCTTTTCCGCCtt                                                                                           <  1:290532/149‑1 (MQ=255)                                                        tgtgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGAt                                                                                      <  8:51795/149‑1 (MQ=255)                                                             gcTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGCACCATTGCTATCGCCTATGTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATggg                                                                                   <  1:100785/149‑1 (MQ=255)                                                              cTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCATCGGTTATATCTTTTCCGCCATCGGATGGGCt                                                                                 >  8:73979/1‑148 (MQ=255)                                                               tACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTa                                                                               >  4:187112/1‑149 (MQ=255)                                                                  cTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGt                                                                            >  5:167708/1‑149 (MQ=255)                                                                        cTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctg                                                                       <  6:285098/148‑1 (MQ=255)                                                                           ttCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGcc                                                                    <  1:34704/148‑1 (MQ=255)                                                                             cATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                     >  6:79185/1‑145 (MQ=255)                                                                             cATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                     <  5:79185/145‑1 (MQ=255)                                                                              aTGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGAt                                                                >  4:145946/1‑149 (MQ=255)                                                                              aTGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGAt                                                                >  8:63926/1‑149 (MQ=255)                                                                                gCTGTTTCTGGTGACAACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGAACAAACCCAACGCCAATGCCGCCGAAACCCGCGGTTATATCTTTTCCGCATACGGATGGGCATACGAGGCTCGCCAGAtac                                                              >  2:85110/1‑147 (MQ=255)                                                                                  tGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctgg                                                                      >  3:80522/1‑139 (MQ=255)                                                                                  tGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctgg                                                                      <  4:80522/139‑1 (MQ=255)                                                                                          gTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctg                                                    >  8:119785/1‑149 (MQ=255)                                                                                          gTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCACATTCCtggtgg                                                       <  2:277068/146‑1 (MQ=255)                                                                                                  gATCAACTATGCCGACCGCGCG‑ACCAGTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTc                                             >  1:171695/1‑148 (MQ=255)                                                                                                      aaCTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATc                                         <  3:187112/148‑1 (MQ=255)                                                                                                      aaCTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATT‑CTGGTGGCTGGCTACTCGATCGc                                       >  1:309471/1‑149 (MQ=255)                                                                                                        cTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCt                                      <  7:63926/149‑1 (MQ=255)                                                                                                            gCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCATCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCgg                                  >  5:316182/1‑149 (MQ=255)                                                                                                             ccGACGGCTCG‑ACCATTGCTTTTG‑TTGGGGCCTGCAGGAAGAAATACCTCGGCATTGCCGGCTTAACCCTCGGTTATATCTTTTCCGCCTTCGGTTGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTGGGt                                 <  1:85110/149‑1 (MQ=255)                                                                                                                 ccGCGCG‑ACCATTGCTTTCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGTGCTTACGTGGCTGGCCAGATTCCCGGTGGCTGGCTACTCGATCGCTTCGGTTCCa                             <  7:73979/149‑1 (MQ=255)                                                                                                                  cgcgcg‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCaa                            <  6:167708/149‑1 (MQ=255)                                                                                                                   gcgcg‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTGCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAg                           <  3:145946/149‑1 (MQ=255)                                                                                                                   gcgcg‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTcc                              >  5:185161/1‑146 (MQ=255)                                                                                                                   gcgcg‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTcc                              >  1:205368/1‑146 (MQ=255)                                                                                                                       g‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACg                       >  4:98993/1‑149 (MQ=255)                                                                                                                                cTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGc                <  3:98993/148‑1 (MQ=255)                                                                                                                                    tg‑CAGGCTGGAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCtt             >  3:40363/2‑147 (MQ=255)                                                                                                                                     g‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCa           <  7:119785/148‑1 (MQ=255)                                                                                                                                              cAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCa  >  5:217007/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTTCA  >  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: ‑