Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,106,518 (C)6→7 coding (709/1167 nt) wbbH ← putative O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,106,5121.C95.6% 156.2 / 0.1 45coding (715/1167 nt)wbbHputative O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (16/27);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑CCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACA  >  NC_000913/2106373‑2106645
                                                                                                                                            |                                                                                                                                      
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  1:185318/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  1:86560/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  1:756922/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  2:272498/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  2:393556/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         >  1:604926/1‑140 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  2:897985/140‑1 (MQ=255)
aCTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  1:30895/140‑1 (MQ=255)
  tCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑cc                                                                                                                                     >  1:88699/1‑140 (MQ=255)
   ccACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         <  1:437570/137‑1 (MQ=255)
   ccACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                                                         >  2:437570/1‑137 (MQ=255)
   ccACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑ccc                                                                                                                                    >  2:860001/1‑140 (MQ=255)
   ccACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑ccc                                                                                                                                    >  1:837844/1‑140 (MQ=255)
             aTAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAg                                                                                                                           <  1:296027/140‑1 (MQ=255)
               aGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTc                                                                                                                         <  1:568354/140‑1 (MQ=255)
                 ccccATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAGAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTc                                                                                                                         >  2:754066/1‑138 (MQ=255)
                 ccccATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAGAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTc                                                                                                                         <  1:754066/138‑1 (MQ=255)
                    cATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGt                                                                                                                    >  2:496930/1‑140 (MQ=255)
                    cATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGt                                                                                                                    >  1:807091/1‑140 (MQ=255)
                    cATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGt                                                                                                                    >  1:810679/1‑140 (MQ=255)
                              tCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCa                                                                                                          <  2:88699/140‑1 (MQ=255)
                              tCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCa                                                                                                          <  2:196443/140‑1 (MQ=255)
                                      ccAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAg                                                                                                  <  2:355003/140‑1 (MQ=255)
                                          aaGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGgaga                                                                                              <  2:850858/140‑1 (MQ=255)
                                            gACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAAt                                                                                            <  1:708888/140‑1 (MQ=255)
                                              cATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAg                                                                                          <  1:194577/140‑1 (MQ=255)
                                                      cTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                                                                  <  1:51852/140‑1 (MQ=255)
                                                      cTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                                                                  <  2:919935/140‑1 (MQ=255)
                                                      cTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                                                                  >  2:63125/1‑140 (MQ=255)
                                                      cTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                                                                  <  2:238953/140‑1 (MQ=255)
                                                       tGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                                                                  <  2:20797/140‑1 (MQ=255)
                                                       tGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTAc                                                                                 >  1:697506/1‑140 (MQ=255)
                                                                tAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATAGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaat                                                                        <  2:807091/140‑1 (MQ=255)
                                                                    tACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGa                                                                    >  2:205879/1‑140 (MQ=255)
                                                                       ttGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑CCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATaaa                                                                >  1:510456/1‑140 (MQ=255)
                                                                         gCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAAc                                                               >  1:556697/1‑140 (MQ=255)
                                                                         gCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAAc                                                               >  1:810793/1‑140 (MQ=255)
                                                                          cTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACa                                                              <  2:343836/140‑1 (MQ=255)
                                                                                   aaaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGtttt                                                     <  1:695066/140‑1 (MQ=255)
                                                                                    aaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTa                                                    <  2:708675/140‑1 (MQ=255)
                                                                                      cTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTAct                                                  <  1:63125/140‑1 (MQ=255)
                                                                                         cTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTa                                               >  1:308459/1‑140 (MQ=255)
                                                                                           gAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTAtt                                             <  1:619521/140‑1 (MQ=255)
                                                                                           gAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTAtt                                             <  1:324421/140‑1 (MQ=255)
                                                                                                         gggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGa                               <  2:810679/140‑1 (MQ=255)
                                                                                                          ggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGaa                              <  2:308459/140‑1 (MQ=255)
                                                                                                          ggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGaa                              >  2:392328/1‑140 (MQ=255)
                                                                                                          ggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGCAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGaa                              <  2:837844/140‑1 (MQ=255)
                                                                                                            cTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATg                            <  1:392328/140‑1 (MQ=255)
                                                                                                               aCCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCAt                         >  1:847556/1‑140 (MQ=255)
                                                                                                                  aaaTACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAag                      <  1:496930/140‑1 (MQ=255)
                                                                                                                   aaTACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAaga                     <  2:714927/140‑1 (MQ=255)
                                                                                                                   aaTACATGGATAGATAATATGCCAT‑‑CCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAagag                    <  2:510456/140‑1 (MQ=255)
                                                                                                                          ggATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATa              >  1:790767/1‑140 (MQ=255)
                                                                                                                          ggATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATa              <  2:604926/140‑1 (MQ=255)
                                                                                                                                ataataTGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATa        >  2:632118/1‑140 (MQ=255)
                                                                                                                                     atGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGAc   >  1:739145/1‑140 (MQ=255)
                                                                                                                                      tGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACa  >  2:373084/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                      
ACTCCACCTGTTAATAGCCCCATCAGCCTTTCAAAAAACCAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAT‑‑CCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACA  >  NC_000913/2106373‑2106645

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