Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 417,375 C→T 18.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,3750CT18.8% 59.0 / 10.0 32P361L (CCG→CTG) waaGlipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base T (3/3);  total (17/15)
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.93e-01

CGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCC  >  NC_002947/417229‑417523
                                                                                                                                                  |                                                                                                                                                     
cGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAAAAAGTGCCCATCATCAGCCATCTGCAGCGCATGCTCTTAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCAtgctgcg                                                                                                                                                    >  1:575037/1‑149 (MQ=255)
 gACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCTTGCTATGGTCTGGCATGTGCTGATAGCGCCGATCTGTACAGCCTGCCgcag                                                                                                                                                   >  2:312346/1‑150 (MQ=255)
        gCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTAAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtgctg                                                                                                                                            <  1:290260/150‑1 (MQ=255)
           tGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtac                                                                                                                                         <  2:575037/150‑4 (MQ=255)
               gCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGt                                                                                                                                     >  1:115115/1‑150 (MQ=255)
                            ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta                                                                                                                                          >  2:750282/1‑130 (MQ=255)
                            ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta                                                                                                                                          <  1:750282/132‑3 (MQ=255)
                            ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCaa                                                                                                                        >  2:332571/1‑150 (MQ=255)
                                aaCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                            <  1:564929/142‑1 (MQ=255)
                                aaCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                            >  2:564929/1‑142 (MQ=255)
                                   aaGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCg                                                                                                                 <  2:621237/150‑1 (MQ=255)
                                       gcagcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATg                                                                                                             >  1:209930/1‑150 (MQ=255)
                                                            gcgcATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAccc                                                                                        >  1:605847/1‑150 (MQ=255)
                                                            gcgcATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGACCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAccc                                                                                        >  1:809173/1‑150 (MQ=255)
                                                                aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCa                                                                                    >  2:628191/1‑150 (MQ=255)
                                                                   cTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGACCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAAc                                                                                 <  2:809173/150‑1 (MQ=255)
                                                                    tCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACg                                                                                >  1:205826/1‑150 (MQ=255)
                                                                      gAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCCATCTGTACAGCATGCCGAAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGcc                                                                             <  1:700069/150‑1 (MQ=255)
                                                                       aaGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCt                                                                            >  2:744196/1‑150 (MQ=255)
                                                                             cGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta                                                                                                                                          >  2:638964/1‑81 (MQ=39)
                                                                             cGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta                                                                                                                                          <  1:638964/83‑3 (MQ=39)
                                                                                         ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGc                                                          <  1:628191/150‑1 (MQ=255)
                                                                                          cAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGcc                                                         >  2:191208/1‑150 (MQ=255)
                                                                                                               ggCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTg                                    <  1:332571/150‑1 (MQ=255)
                                                                                                                       cTGATAGCGCCGATCTGTAAAGTCTGCCGCTGCATGCAGCCGCTGTGATCCCAGCCCACCCCACCGCTCGGGGCTGATTGTTGCCGATGTTTTTAAATGCTGGTTGGGCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGa                            <  2:736500/150‑1 (MQ=255)
                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                           <  1:91453/150‑1 (MQ=255)
                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                           >  1:736500/1‑150 (MQ=255)
                                                                                                                              cgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGCTCTCGACCCATCCCTTAAACGCTGTTTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTAt                     <  1:744196/150‑1 (MQ=255)
                                                                                                                                   aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGa                >  1:494590/1‑150 (MQ=255)
                                                                                                                                       cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTg            >  1:485748/4‑150 (MQ=255)
                                                                                                                                        tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGa                                                    >  2:4541/1‑109 (MQ=255)
                                                                                                                                        tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGa                                                    <  1:4541/109‑1 (MQ=255)
                                                                                                                                            ccaagCCGAAGCATGCTGCCGACGTGATCCTTGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGAGACGGGCGGGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAgg       <  2:209930/146‑1 (MQ=255)
                                                                                                                                                 ccGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGcgcc  >  2:532571/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCC  >  NC_002947/417229‑417523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: