Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	85,077	T→G	12.9%	intergenic (+55/+75)	aroE → / ← PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	85,077	0	T	G	12.9%	112.7 / 11.4	54	intergenic (+55/+75)	aroE/PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base G (3/4);  total (23/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 = 5.40e-01

CAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCT  >  NC_002947/84930‑85223                                                                                                                                                    |                                                                                                                                                    caggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCCGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTa                                                                                                                                                    >  6:283171/1‑149 (MQ=255)  aggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTTAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑c                                                                                                                                                  >  4:213553/1‑149 (MQ=255)        cTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTa                                                                                                                                            <  2:76359/149‑1 (MQ=255)         tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGCCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTAc                                                                                                                                           >  6:87380/1‑149 (MQ=255)                ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATc                                                                                                                                    >  1:56113/1‑149 (MQ=255)                ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCC‑GGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATcc                                                                                                                                   >  1:229484/1‑149 (MQ=255)                  cttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCACTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCct                                                                                                                                  >  1:167162/1‑149 (MQ=255)                       tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                                                                              <  5:166151/148‑1 (MQ=255)                       tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                                                                              <  5:124593/148‑1 (MQ=255)                             gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCg                                                                                                                       >  5:104323/1‑149 (MQ=255)                              tggcgtgAGGCTGGATACGGCGCCGGTGTTGGATGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACgg                                                                                                                      <  4:220408/149‑1 (MQ=21)                               ggcgtgAGGCCGGATACGGCGCCGGTTTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                                                                     <  3:273437/149‑1 (MQ=255)                                  gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACGGCCtt                                                                                                                  <  2:75527/149‑1 (MQ=21)                                   tgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACGGCCTTg                                                                                                                 >  7:80393/1‑149 (MQ=21)                                     aGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc                                                                                                                <  2:274105/148‑1 (MQ=255)                                      ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACGGCCTTGCGt                                                                                                              <  8:80393/149‑1 (MQ=21)                                       gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGGCGTTCCCGCAAAGTA‑CGTCGTACGCCGCTCCGCTGTGGGAGCGGCCGTGCGtg                                                                                                             >  5:207252/1‑148 (MQ=39)                                       gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc                                                                                                             <  7:122378/149‑1 (MQ=255)                                       gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc                                                                                                             >  2:243775/1‑149 (MQ=255)                                                   gcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCCCAAAGTA‑CGTTGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                                                                  >  1:123106/1‑148 (MQ=255)                                                   gcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGCCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGc                                                                                                 <  7:2544/149‑1 (MQ=255)                                                     cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTACCACAAAGTA‑cgt                                                                                                                                                >  4:98358/1‑100 (MQ=255)                                                         gtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTAGTAAGCCCAACCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                           >  4:112842/1‑149 (MQ=255)                                                          ttGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                           <  1:243775/148‑1 (MQ=255)                                                            ggCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCGGTGGGAGCGGCCTTGCGTCGCGATGGGc                                                                                                    >  3:171519/1‑137 (MQ=255)                                                              cTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAccc                                                                                      <  2:56113/149‑1 (MQ=255)                                                                gAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAcc                                                                                       >  7:125579/1‑146 (MQ=255)                                                                  aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                  <  7:35256/149‑1 (MQ=255)                                                                      cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CCTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                <  5:283171/147‑1 (MQ=255)                                                                       gcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCa                                                                             <  6:104323/149‑1 (MQ=255)                                                                           gCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCt                                                                         >  4:88072/1‑149 (MQ=255)                                                                              acgcgctggGGGCTGAGCGGGCCGGGGCCGCTTGGCGGCCCTATCGCGCCGCACGGCCGTCCCCACAATGTA‑CGGCGTACGCCGATCCGCTGTGGCAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAt                                                                      <  8:285958/146‑1 (MQ=38)                                                                                  gctcgGGGCGGAGGTGGCCGGGGCCGGTTGGCGGCGCTATCGGGACGCAAGGGCGTTCCCACAAAGTA‑CGTCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                  <  5:87380/149‑1 (MQ=38)                                                                                       gggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAAc                                                             <  3:213553/149‑1 (MQ=255)                                                                                       gggCTGAGCTGGCCGGGGCCACTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGGGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAAc                                                             <  2:167162/149‑1 (MQ=255)                                                                                             aGCTGGCC‑GGGCCGCTTTGCGTCCCTATCGCCACGCAAGGCCTTTCCCACAACGTA‑CGTCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATc                                                       <  2:229484/148‑1 (MQ=38)                                                                                                           gCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                         >  4:73637/1‑149 (MQ=255)                                                                                                                  ggcCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTcc                                  <  1:208345/149‑1 (MQ=255)                                                                                                                         tCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt                           <  7:115845/149‑1 (MQ=255)                                                                                                                           gcgACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                         <  3:73637/149‑1 (MQ=255)                                                                                                                             gACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                         >  6:23297/1‑147 (MQ=255)                                                                                                                             gACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                         <  5:23297/147‑1 (MQ=255)                                                                                                                              aCGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                      <  3:112842/149‑1 (MQ=255)                                                                                                                                 cacggCCGCCCCCACAACGAC‑CGCCGTCCGCCCAACCTCTGTGGCAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                   <  3:174175/146‑1 (MQ=37)                                                                                                                                    ggCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                 >  8:72917/1‑148 (MQ=25)                                                                                                                                    gcccGTTCCCACAACGAC‑CCTCGTACGCCCACCCACTGTGGGACCGGCGTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                <  3:24058/147‑1 (MQ=255)                                                                                                                                      ccGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                    >  3:164887/1‑111 (MQ=14)                                                                                                                                      ccGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                    <  4:164887/111‑1 (MQ=14)                                                                                                                                        gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                            >  5:93321/1‑133 (MQ=255)                                                                                                                                        gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                            <  6:93321/133‑1 (MQ=255)                                                                                                                                         ttCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc           <  3:88072/149‑1 (MQ=255)                                                                                                                                         ttCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc           <  8:125579/149‑1 (MQ=255)                                                                                                                                             cacaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc       >  5:6744/1‑149 (MQ=255)                                                                                                                                                  aGTA‑CGTCGTACGCAGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCt  >  4:240483/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCT  >  NC_002947/84930‑85223

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