Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,090,056	C→G	11.3%	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	11.3%	147.4 / 10.3	53	S36C (TCC→TGC) ‡	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 = 5.00e-01

CGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCA  >  NC_002947/4089909‑4090196                                                                                                                                                      |                                                                                                                                             cGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATTCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTc                                                                                                                                              <  2:58540/148‑1 (MQ=255)          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                                                                                                                                      <  1:63650/120‑1 (MQ=255)                                     ccAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGc                                                                                                                                      >  2:63650/1‑120 (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)                                                                                                                                                      |                                                                                                                                             CGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCG‑ACCATTGCTATCG‑CAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCA  >  NC_002947/4089909‑4090196

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