Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	85,040	G→A	8.0%	intergenic (+18/+112)	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,040	0	G	A	8.0%	140.6 / 2.4	50	intergenic (+18/+112)	aroE/PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

GCCAAGGTACTGGATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCA  >  NC_002947/84894‑85175                                                                                                                                                    |                                                                                                                                         gCCAAGGTACTGGATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTggc                                                                                                                                        >  5:288322/1‑149 (MQ=255)   cAAGGTACTGGATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTggcgg                                                                                                                                      >  2:304733/1‑149 (MQ=255)     aGGTACTGGATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGcc                                                                                                                                    <  6:62129/149‑1 (MQ=255)             gATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATcgc                                                                                                                             >  4:272332/1‑148 (MQ=255)             gATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATcgc                                                                                                                             >  4:285304/1‑148 (MQ=255)                 ggTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGc                                                                                                                        <  7:172491/149‑1 (MQ=255)                        gATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGGCGCCAGGCCg                                                                                                                 >  8:169859/1‑149 (MQ=255)                                  gAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACaaa                                                                                                       >  3:14120/1‑149 (MQ=255)                                     caggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGt                                                                                                     >  8:10002/1‑148 (MQ=255)                                         cagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑cgtc                                                                                               <  3:285304/149‑1 (MQ=255)                                             tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTa                                                                                            >  1:224245/1‑149 (MQ=255)                                                aGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGc                                                                                                                   >  3:224504/1‑123 (MQ=255)                                                aGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGc                                                                                                                   <  4:224504/123‑1 (MQ=255)                                                aGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGc                                                                                         <  4:14120/148‑1 (MQ=255)                                                 gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCCAAGCGGccc                                                                                                                                   >  6:151564/1‑106 (MQ=255)                                                 gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCCAAGCGGccc                                                                                                                                   <  5:151564/106‑1 (MQ=255)                                                    ttcttTATCTGGAGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGAt                                                                                     >  3:62564/1‑148 (MQ=255)                                                        ttATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCCAAGCGGcc                                                                                                                                    >  8:291641/1‑98 (MQ=255)                                                        ttATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCCAAGCGGcc                                                                                                                                    <  7:291641/98‑1 (MQ=255)                                                          aTCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                              >  7:113917/1‑149 (MQ=255)                                                          aTCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                              >  5:2885/1‑149 (MQ=255)                                                           tCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTc                                                                                                              >  5:271838/1‑117 (MQ=255)                                                           tCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTc                                                                                                              <  6:271838/117‑1 (MQ=255)                                                           tCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                              <  3:272332/148‑1 (MQ=255)                                                               gcgtggcgtgAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg                                                                             <  4:62564/145‑1 (MQ=255)                                                               gcgtggcgtgAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg                                                                             <  6:288322/145‑1 (MQ=255)                                                               gcgtggcgtgAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGa                                                                          <  7:10002/148‑1 (MQ=255)                                                               gcgtggcgtgAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCC‑ATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGaac                                                                        <  3:81529/149‑3 (MQ=21)                                                                          ggCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGCTCCTCTGTGGGAGCGGCCTTGCGt                                                              >  1:234621/1‑149 (MQ=255)                                                                          ggCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGt                                                              <  1:304733/149‑1 (MQ=255)                                                                              ggATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                              <  2:228450/129‑1 (MQ=255)                                                                                aTACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                        <  2:224245/149‑1 (MQ=255)                                                                                  aCGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                     <  7:68053/134‑1 (MQ=255)                                                                                  aCGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                     >  8:68053/1‑134 (MQ=255)                                                                                    ggCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATc                                                                                    <  6:79167/117‑1 (MQ=255)                                                                                    ggCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATc                                                                                    >  5:79167/1‑117 (MQ=255)                                                                                      cg‑cCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                  <  2:234621/149‑1 (MQ=255)                                                                                              gtgggCGGAACGGCGCCGGCAACGCGCGCGGGGCTGAGCGGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTCCCCACAACGTA‑CGGCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                           <  7:169859/146‑1 (MQ=37)                                                                                               ttGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg                                                                             >  6:256606/1‑114 (MQ=255)                                                                                               ttGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg                                                                             <  5:256606/114‑1 (MQ=255)                                                                                                  gCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACCCAAGGCCGGTACCCCAAAGTA‑CGGGGTACGCCGCACCGCGGGGGGGGCGGCCGTGCGGGGCGAGGGGGGGCAAAGccgccc                                       >  7:198213/1‑144 (MQ=255)                                                                                                        cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTc                                                                                                              >  7:192694/1‑73 (MQ=255)                                                                                                        cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTc                                                                                                              <  8:192694/73‑1 (MQ=255)                                                                                                        cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGCCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAt                                 <  6:2885/149‑1 (MQ=255)                                                                                                         tGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCCACGCAAGGCCGTTCCCACAAACTA‑CGTCGTACGCCGATCCTCTGTGGGCGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                >  7:196153/1‑149 (MQ=255)                                                                                                              cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCCCAAAGTA‑CGTCGTACGCCGCTCCTCTGTGGGCGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGaagacccac                           >  3:178519/1‑142 (MQ=38)                                                                                                               ggCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATGCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCCAcc                          >  5:94390/1‑149 (MQ=255)                                                                                                                gCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCt                         >  7:19168/1‑149 (MQ=255)                                                                                                                                  aGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATc       >  5:59786/1‑149 (MQ=255)                                                                                                                                     tGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg    <  3:155889/149‑1 (MQ=255)                                                                                                                                      gggcGGGGCGGTTTGGCGGCCCTTTCGGTCCTCATGGCCTTACCCACAATTTA‑GTTGGTAGGCGGTTCTTTTTTGGTTGCGGCCTTGCTTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGc   <  8:196153/146‑1 (MQ=255)                                                                                                                                       gCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca  >  5:71230/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                         GCCAAGGTACTGGATGGGTTGGGGATGCTGGCTGAGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCG‑CCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCA  >  NC_002947/84894‑85175

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