Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	417,375	A→G	11.6%	Q362R (CAG→CGG)	waaG →	lipopolysaccharide glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	417,375	0	A	G	11.6%	105.3 / 6.0	43	Q362R (CAG→CGG)	waaG	lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (3/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

ACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417228‑417519                                                                                                                                                     |                                                                                                                                                  aCAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgcag                                                                                                                                                  <  7:192257/149‑1 (MQ=255)          tGGTGCTGGACGAGCCGTTCGAGCAAGAGCGGCTCAGCGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGcc                                                                                                                                         >  2:188042/1‑149 (MQ=255)             tGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGAc                                                                                                                                      <  8:248187/149‑1 (MQ=255)             tGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGAc                                                                                                                                      <  5:75605/149‑1 (MQ=255)                   aCGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATc                                                                                                                                <  3:96617/149‑1 (MQ=255)                      aGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTg                                                                                                                             <  8:40606/149‑1 (MQ=255)                             cGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCTATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa                                                                                                                        >  1:44243/1‑148 (MQ=255)                                aCAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacaga                                                                                                                                     >  6:109531/1‑125 (MQ=255)                                aCAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacaga                                                                                                                                     <  5:109531/131‑7 (MQ=255)                                    gagCAGCTCAACGGCTATCTGCAGCGC‑ATGCTAGAAGAAACGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGAt                                                                                                                                 >  7:50674/1‑131 (MQ=255)                                         gcTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGaa                                                                                                          <  1:188042/149‑1 (MQ=255)                                           tCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGc                                                                                                        >  8:306378/1‑149 (MQ=255)                                             aaCGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCt                                                                                                       >  5:238127/1‑148 (MQ=255)                                             aaCGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTg                                                                                                      <  3:139783/149‑1 (MQ=255)                                                     tCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑gg                                                                                              <  2:44243/148‑1 (MQ=255)                                                                aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGTAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCtt                                                                                    >  5:32940/1‑148 (MQ=255)                                                                    tCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAAc                                                                               >  4:160933/1‑149 (MQ=255)                                                                      gAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGc                                                                             >  8:26420/1‑149 (MQ=255)                                                                        aGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCt                                                                           >  7:261910/1‑149 (MQ=255)                                                                            ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGGCGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTgg                                                                       >  2:204511/1‑149 (MQ=255)                                                                            ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGcccggcc                                                                       >  7:89218/1‑144 (MQ=255)                                                                            ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTgg                                                                       >  6:269597/1‑149 (MQ=255)                                                                                  gcgcGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGGCGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCggg                                                                 <  1:204511/149‑1 (MQ=255)                                                                                         ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGAt                                                           <  1:222748/148‑1 (MQ=255)                                                                                         ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATg                                                          <  5:269597/149‑1 (MQ=255)                                                                                            gCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCt                                                       <  7:306378/149‑1 (MQ=255)                                                                                                gTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGa                                                   >  5:1586/1‑149 (MQ=255)                                                                                                   gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCAt                                                                                                             <  8:133518/89‑1 (MQ=255)                                                                                                   gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCAt                                                                                                             >  7:133518/1‑89 (MQ=255)                                                                                                            ccgggCATTTGCTTAAACCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTCGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        <  1:255942/146‑1 (MQ=255)                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGc                             >  7:231196/1‑147 (MQ=255)                                                                                                                          aTAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGcaa                           <  7:26420/147‑3 (MQ=255)                                                                                                                          aTAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATACTGGCGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAAACCTTCAAACGCCTGTGGGCCGGTAGTGATGCCTTCGAAGCAGTGGAGGCGCTGAAAGGCGAg                          >  8:215189/1‑148 (MQ=255)                                                                                                                                  ctgatGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgc                 <  8:11947/145‑1 (MQ=255)                                                                                                                                   aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgtg                >  7:213909/1‑147 (MQ=255)                                                                                                                                      tGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGc             >  3:81572/1‑149 (MQ=255)                                                                                                                                       cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAg                                         <  6:120349/117‑1 (MQ=255)                                                                                                                                       cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAg                                         >  5:120349/4‑120 (MQ=255)                                                                                                                                        tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGc             <  6:238127/147‑1 (MQ=255)                                                                                                                                         gcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTg           >  2:24856/2‑148 (MQ=255)                                                                                                                                         aCAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCt            <  4:509/147‑1 (MQ=255)                                                                                                                                         aCAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCt            >  3:509/1‑147 (MQ=255)                                                                                                                                              aTGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAggg     <  8:231196/149‑1 (MQ=255)                                                                                                                                                 ccGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGgcgc  >  7:89415/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGC‑ATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATT‑GGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417228‑417519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑