Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	85,088	C→A	13.3%	intergenic (+66/+64)	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,088	0	C	A	13.3%	136.3 / 10.2	51	intergenic (+66/+64)	aroE/PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (18/26);  new base A (3/4);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGG  >  NC_002947/84945‑85230                                                                                                                                                 |                                                                                                                                               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)                                          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                                <  5:23297/147‑1 (MQ=255)                                                                                                              gACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                                >  6:23297/1‑147 (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                  <  8:125579/149‑1 (MQ=255)                                                                                                                          ttCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc                  <  3:88072/149‑1 (MQ=255)                                                                                                                              cacaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc              >  5:6744/1‑149 (MQ=255)                                                                                                                                   aGTA‑CGTCGTACGCAGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCt         >  4:240483/1‑149 (MQ=255)                                                                                                                                          tcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCg   >  5:168237/1‑149 (MQ=255)                                                                                                                                           cgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCACATTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCgg  >  5:2324/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               TTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGG  >  NC_002947/84945‑85230

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