Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 416,121 G→C 7.7% D297H (GAC→CAC)  waaC → ADP‑heptose‑‑LPS heptosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947416,1210GC 7.7% 162.4 / 3.5 52D297H (GAC→CAC) waaCADP‑heptose‑‑LPS heptosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base C (2/2);  total (22/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 = 7.71e-01

GCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGT  >  NC_002947/415977‑416265
                                                                                                                                                 |                                                                                                                                                 
gcggcggcAAAGGCTTGCTTGGCGGTCGATACGGGCCTTGGCCACCAGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGAc                                                                                                                                                 <  5:173615/147‑1 (MQ=255)
gcggcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTg                                                                                                                                               <  5:213529/149‑1 (MQ=255)
  ggcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGc                                                                                                                                             <  5:11371/149‑1 (MQ=255)
  ggcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGc                                                                                                                                             >  8:67739/1‑149 (MQ=255)
     ggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCtt                                                                                                                                          <  7:41028/149‑1 (MQ=255)
      gcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTg                                                                                                                                         >  4:1022/1‑149 (MQ=255)
       cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGCGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtgg                                                                                                                                         <  1:60818/149‑3 (MQ=255)
          aGGCCTGCGGGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCt                                                                                                                                     <  7:87302/149‑1 (MQ=255)
                  gTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcct                                                                                                                             <  8:26859/149‑1 (MQ=255)
                     gCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCa                                                                                                                          <  5:90296/149‑1 (MQ=255)
                        gTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaa                                                                                                                       <  7:67739/149‑1 (MQ=255)
                                cGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCa                                                                                                               <  6:111010/149‑1 (MQ=255)
                                             cTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAg                                                                                                  >  7:26763/1‑149 (MQ=255)
                                                 cggcggcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCTcc                                                                                              <  6:189149/149‑1 (MQ=255)
                                                         cTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTcc                                                                                                                                   <  4:327292/104‑1 (MQ=255)
                                                         cTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTcc                                                                                                                                   >  3:327292/1‑104 (MQ=255)
                                                          tCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCTGGGCCTGACTGGCGCCTACG‑GACTGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCAACTACAAACCGAGCGCCGATGACCTg                                                                                     >  5:295473/1‑149 (MQ=255)
                                                                  cccACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGtt                                                                             >  7:71846/1‑149 (MQ=255)
                                                                  cccACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGtt                                                                             <  8:26763/149‑1 (MQ=255)
                                                                       ccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGa                                                                                            >  4:146263/1‑129 (MQ=255)
                                                                       ccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGa                                                                                            <  3:146263/129‑1 (MQ=255)
                                                                       ccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATc                                                                        >  3:186479/1‑149 (MQ=255)
                                                                           tCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTg                                                                      <  3:1022/147‑1 (MQ=255)
                                                                             gCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGaaa                                                                   <  7:3642/148‑1 (MQ=255)
                                                                             gCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAc                                                                  <  8:302783/149‑1 (MQ=255)
                                                                                 tttGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAcgc                                                                >  8:239157/1‑147 (MQ=255)
                                                                                 tttGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAcgc                                                                <  4:186479/147‑1 (MQ=255)
                                                                                    ggCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGt                                                            <  7:311058/148‑1 (MQ=255)
                                                                                     gCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTgg                                                          >  7:53988/1‑149 (MQ=255)
                                                                                      cccTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGc                                                         <  2:132998/149‑1 (MQ=255)
                                                                                       ccTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGAcc                                                                                       >  2:151977/1‑118 (MQ=255)
                                                                                       ccTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGAcc                                                                                       <  1:151977/118‑1 (MQ=255)
                                                                                           ccAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCtg                                                    >  2:100818/1‑149 (MQ=255)
                                                                                            cAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGc                                                         >  5:236661/1‑143 (MQ=255)
                                                                                            cAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGc                                                         <  6:236661/143‑1 (MQ=255)
                                                                                                         aCTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGcc                                                                                  <  2:225220/105‑1 (MQ=255)
                                                                                                         aCTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCt                                     >  2:185496/1‑149 (MQ=255)
                                                                                                           tGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCATGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCCCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTg                                    >  8:189538/1‑148 (MQ=255)
                                                                                                             gcgcCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAAt                                 >  8:85184/1‑149 (MQ=255)
                                                                                                                ccTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATccc                              <  7:85184/149‑1 (MQ=255)
                                                                                                                    cg‑gacggaCCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGa                            <  8:71846/147‑1 (MQ=255)
                                                                                                                     g‑gacggacCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCa                         >  4:120373/1‑149 (MQ=255)
                                                                                                                                 gATTCCCCCGGCCAGTGTCTTGCCGTGCGCTCCCTGCCGGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGcc              <  4:234822/149‑1 (MQ=255)
                                                                                                                                      aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCtt          <  3:120373/148‑1 (MQ=255)
                                                                                                                                        aaGGCCAGTCACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTg         <  5:54444/146‑1 (MQ=255)
                                                                                                                                        aaGGCCAGTCACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTg         >  6:54444/2‑147 (MQ=255)
                                                                                                                                        aaGGCCAGTCACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTGAg       >  8:55238/2‑149 (MQ=255)
                                                                                                                                          ggCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAgc      <  1:218094/149‑1 (MQ=255)
                                                                                                                                          ggCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCg                  >  4:54967/1‑136 (MQ=255)
                                                                                                                                            ccAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTc                                                                            >  8:112459/1‑77 (MQ=255)
                                                                                                                                            ccAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTc                                                                            <  7:112459/77‑1 (MQ=255)
                                                                                                                                             cAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGt  >  3:177298/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACG‑GACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCG‑CCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGT  >  NC_002947/415977‑416265

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