Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 417,372 C→T 10.8% P361L (CCG→CTG)  waaG → lipopolysaccharide glucosyltransferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947417,3720CT10.8% 121.2 / 5.2 46P361L (CCG→CTG) waaGlipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417226‑417519
                                                                                                                                                   |                                                                                                                                                    
cGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgc                                                                                                                                                    <  1:225889/149‑1 (MQ=255)
cGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCg                                                                                                                                                     >  2:222748/1‑148 (MQ=255)
 gACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCg                                                                                                                                                     >  8:192257/1‑147 (MQ=255)
  aCAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgcag                                                                                                                                                  <  7:192257/149‑1 (MQ=255)
           tGGTGCTGGACGAGCCGTTCGAGCAAGAGCGGCTCAGCGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGcc                                                                                                                                         >  2:188042/1‑149 (MQ=255)
              tGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGAc                                                                                                                                      <  8:248187/149‑1 (MQ=255)
              tGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGAc                                                                                                                                      <  5:75605/149‑1 (MQ=255)
                    aCGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATc                                                                                                                                <  3:96617/149‑1 (MQ=255)
                       aGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTg                                                                                                                             <  8:40606/149‑1 (MQ=255)
                              cGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCTATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa                                                                                                                        >  1:44243/1‑148 (MQ=255)
                                 aCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacaga                                                                                                                                     >  6:109531/1‑125 (MQ=255)
                                 aCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacaga                                                                                                                                     <  5:109531/131‑7 (MQ=255)
                                     gagCAGCTCAACGGCTATCTGCAGCGCATGCTAGAAGAAACGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGAt                                                                                                                                 >  7:50674/1‑131 (MQ=255)
                                          gcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGaa                                                                                                          <  1:188042/149‑1 (MQ=255)
                                            tCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGc                                                                                                        >  8:306378/1‑149 (MQ=255)
                                              aaCGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCt                                                                                                       >  5:238127/1‑148 (MQ=255)
                                              aaCGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTg                                                                                                      <  3:139783/149‑1 (MQ=255)
                                                      tCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTgg                                                                                              <  2:44243/148‑1 (MQ=255)
                                                                 aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGTAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCtt                                                                                    >  5:32940/1‑148 (MQ=255)
                                                                     tCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAAc                                                                               >  4:160933/1‑149 (MQ=255)
                                                                       gAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGc                                                                             >  8:26420/1‑149 (MQ=255)
                                                                         aGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCt                                                                           >  7:261910/1‑149 (MQ=255)
                                                                             ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGGCGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTgg                                                                       >  2:204511/1‑149 (MQ=255)
                                                                             ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGcccggcc                                                                       >  7:89218/1‑144 (MQ=255)
                                                                             ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTgg                                                                       >  6:269597/1‑149 (MQ=255)
                                                                                   gcgcGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGGCGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCggg                                                                 <  1:204511/149‑1 (MQ=255)
                                                                                          ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGAt                                                           <  1:222748/148‑1 (MQ=255)
                                                                                          ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATg                                                          <  5:269597/149‑1 (MQ=255)
                                                                                             gCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCt                                                       <  7:306378/149‑1 (MQ=255)
                                                                                                 gTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGa                                                   >  5:1586/1‑149 (MQ=255)
                                                                                                    gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCAt                                                                                                             >  7:133518/1‑89 (MQ=255)
                                                                                                    gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCAt                                                                                                             <  8:133518/89‑1 (MQ=255)
                                                                                                             ccgggCATTTGCTTAAACCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTCGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        <  1:255942/146‑1 (MQ=255)
                                                                                                                         tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGc                             >  7:231196/1‑147 (MQ=255)
                                                                                                                           aTAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGcaa                           <  7:26420/147‑3 (MQ=255)
                                                                                                                           aTAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATACTGGCGCAAAAAACCGCATGAAGCTGATATTGGCCGAAACCTTCAAACGCCTGTGGGCCGGTAGTGATGCCTTCGAAGCAGTGGAGGCGCTGAAAGGCGAg                          >  8:215189/1‑148 (MQ=255)
                                                                                                                                   ctgatGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgc                 <  8:11947/145‑1 (MQ=255)
                                                                                                                                    aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgtg                >  7:213909/1‑147 (MQ=255)
                                                                                                                                       tGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGc             >  3:81572/1‑149 (MQ=255)
                                                                                                                                        cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAg                                         <  6:120349/117‑1 (MQ=255)
                                                                                                                                        cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAg                                         >  5:120349/4‑120 (MQ=255)
                                                                                                                                         tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGc             <  6:238127/147‑1 (MQ=255)
                                                                                                                                          gcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTg           >  2:24856/2‑148 (MQ=255)
                                                                                                                                          aCAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCt            <  4:509/147‑1 (MQ=255)
                                                                                                                                          aCAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCt            >  3:509/1‑147 (MQ=255)
                                                                                                                                               aTGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAggg     <  8:231196/149‑1 (MQ=255)
                                                                                                                                                  ccGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGgcgc  >  7:89415/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417226‑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: