Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 416,124 G→C 9.8% 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,1240GC 9.8% 151.3 / 7.7 51D297H (GAC→CAC) waaCADP‑heptose‑‑LPS heptosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base C (2/3);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TTCTGGCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGTTGC  >  NC_002947/415975‑416271
                                                                                                                                                       |                                                                                                                                                   
ttCTGGCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTg                                                                                                                                                     >  1:82318/1‑150 (MQ=255)
ttCTGGCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTg                                                                                                                                                     >  1:58758/1‑150 (MQ=255)
   tGGCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTCACt                                                                                                                                                  <  2:353842/150‑1 (MQ=255)
           gcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtggt                                                                                                                                          <  1:189783/150‑4 (MQ=255)
           gcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtggt                                                                                                                                          >  2:189783/1‑147 (MQ=255)
            cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTgcg                                                                                                                                         >  1:164942/1‑150 (MQ=255)
                  ccTGCGTGGCGCTCGATACCGGCCTTGGCGACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGCCCCTACCAACCCGGGCCTGACTGGAGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcct                                                                                                                                   <  2:91887/150‑1 (MQ=255)
                          gCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCag                                                                                                                           >  2:362317/1‑150 (MQ=255)
                            ggTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtgg                                                                                                                                           >  2:349368/1‑130 (MQ=255)
                            ggTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtgg                                                                                                                                           <  1:349368/132‑3 (MQ=255)
                             gTCGATACCGGCCTTGGCCACCTTGCGGCGGCACTCGATGTGCCCATCCTTGCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaag                                                                                                                        >  1:580410/1‑150 (MQ=255)
                             gTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaag                                                                                                                        >  1:83350/1‑150 (MQ=255)
                                        ccTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTAGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTa                                                                                                             >  1:416824/1‑150 (MQ=255)
                                            ggCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACaaa                                                                                                         >  2:375579/1‑150 (MQ=255)
                                             gCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcctg                                                                                                                              >  2:472340/1‑128 (MQ=255)
                                             gCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcctg                                                                                                                              <  1:472340/128‑1 (MQ=255)
                                                aCCTGGCGGCGGCACTCGATGTGCCCACCCTTTAGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGa                                                                                                     <  2:416824/150‑1 (MQ=255)
                                                     gcggcggcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGcc                                                                                                <  2:164942/150‑1 (MQ=255)
                                                     gcggcggcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGcc                                                                                                >  1:769490/1‑150 (MQ=255)
                                                          ggcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGa                                                                                           <  2:88126/150‑1 (MQ=255)
                                                                      gCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGtt                                                                               <  2:58758/150‑1 (MQ=255)
                                                                           cccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATc                                                                          >  2:727077/1‑150 (MQ=255)
                                                                            ccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCt                                                                         <  2:295628/150‑1 (MQ=255)
                                                                            ccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACCGGACCCAAAACCTAC‑GA‑TAAGATTCCCCGGGCCGTGGATGGGCCTTGCGCTCCCTGCCTGGAGGACAAGGGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCt                                                                         >  1:613656/1‑150 (MQ=255)
                                                                              tttCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGa                                                                       >  1:787573/1‑150 (MQ=255)
                                                                                       ttGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGt                                                              <  2:221558/150‑1 (MQ=255)
                                                                                        tGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTg                                                             >  2:181468/1‑150 (MQ=255)
                                                                                         ggCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGACACGCGAGTgg                                                            <  2:554474/150‑1 (MQ=255)
                                                                                              tACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcct                                                                                                                               <  2:462289/78‑1 (MQ=255)
                                                                                              tACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcct                                                                                                                               >  1:462289/1‑78 (MQ=255)
                                                                                               aCCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGa                                                                                                     <  1:324519/103‑1 (MQ=255)
                                                                                               aCCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGa                                                                                                     >  2:324519/1‑103 (MQ=255)
                                                                                                  aaCCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGc                                                   >  1:582220/1‑150 (MQ=255)
                                                                                                  aaCCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGc                                                   >  2:776524/1‑150 (MQ=255)
                                                                                                  aaCCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGc                                                   <  1:464193/150‑1 (MQ=255)
                                                                                                        ggCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGTAACGCGAGTGGCCGCTGTGCTTCACt                                             <  1:727077/150‑1 (MQ=255)
                                                                                                           cTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGc                                          <  2:83350/150‑1 (MQ=255)
                                                                                                                tGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGaa                                     >  1:754229/1‑150 (MQ=255)
                                                                                                                    gcC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTCCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATccc                                 >  1:286992/1‑150 (MQ=255)
                                                                                                                                    ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCggg                   >  1:502038/1‑148 (MQ=255)
                                                                                                                                    ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                   >  1:488990/1‑150 (MQ=255)
                                                                                                                                    ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                   >  2:540218/1‑150 (MQ=255)
                                                                                                                                    ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                   >  1:653628/1‑150 (MQ=255)
                                                                                                                                     cAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCt                                        >  1:284364/1‑128 (MQ=255)
                                                                                                                                     cAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCt                                        <  2:284364/128‑1 (MQ=255)
                                                                                                                                     cAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGc                  <  2:82318/150‑1 (MQ=255)
                                                                                                                                         ttCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAGCGCGAGTGGCCGCTGTGCTTCACTGGCCTGAATCCCGAGCATGGGGCGAGCCGcc              >  2:765381/1‑149 (MQ=255)
                                                                                                                                            aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAg                                                               <  2:436978/98‑1 (MQ=255)
                                                                                                                                            aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAg                                                               >  1:436978/1‑98 (MQ=255)
                                                                                                                                            aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGcc                                         <  2:501341/120‑1 (MQ=255)
                                                                                                                                            aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGcc                                         >  1:501341/1‑120 (MQ=255)
                                                                                                                                                   cAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGtt    <  1:362317/150‑1 (MQ=255)
                                                                                                                                                     gTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGTtgc  <  2:614513/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TTCTGGCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCC‑‑TACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGTTGC  >  NC_002947/415975‑416271

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