Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative colanic acid biosynthesis glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC98.0% 142.9 / ‑5.7 49I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (22/26);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925
                                                                                                                                            |                                                                                                                                           
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                           <  1:102946/143‑1 (MQ=255)
      ggtggtgggTAACCAGCTGTGATTTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     <  2:806109/137‑1 (MQ=255)
        tAGAGGGTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  1:39823/143‑1 (MQ=255)
        tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:962350/143‑1 (MQ=255)
         agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATATTATAAACAACTTGTTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                                  >  1:259788/1‑143 (MQ=255)
           agGGTAAACAGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                <  1:638550/143‑1 (MQ=255)
             ggTTAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                              <  1:848160/143‑1 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATCTGCAGGCTCTGCGTCGCCtctt                                                                                                                              >  1:64153/1‑143 (MQ=255)
                 aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCATTGTGGGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGGGCCTCTTGTAc                                                                                                                          >  1:903909/1‑143 (MQ=255)
                    aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCGACttt                                                                                                                      >  1:814149/1‑143 (MQ=255)
                            tttttgGCTGGCGCTGTCGATTTTGTCTTTGTGCTTGCGGTAAAAGTGGTAGTACCCGGAGAATTGTTTCGGTGATGATGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                              <  2:127472/143‑1 (MQ=255)
                              ttttGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTCAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                            <  1:450134/143‑1 (MQ=255)
                              tgttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                            <  2:865587/141‑1 (MQ=255)
                                   cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTACCCACGGTTCGTCGTa                                                                                                       >  1:302188/1‑143 (MQ=255)
                                    tGGCGCGGTCGAATTTGTCTTTGTGCTTGCAGTAAAAGTGGAATTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCTCGGTTCGCCGTAc                                                                                                      >  1:726602/1‑143 (MQ=255)
                                          ggTCGAATTTGTCTTTGTGCTTCCGGCATTAGTGGACGTACCCGGATAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                <  1:257338/143‑1 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTAATCCGCATTTCACCGTGATTGATATGCCGGCCCTGCAGCCCCTCCTCTCCTTTCCCCCGTTCGCCGCACCccaccaccc                                                                                            >  2:731757/1‑142 (MQ=255)
                                               aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAt                                                                                           >  2:229727/1‑143 (MQ=255)
                                                aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTTGTTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGCTTCGCCGTACTCCCCCACCCTc                                                                                          >  2:546079/1‑143 (MQ=255)
                                                    gttttttTGCTTGCGGTAAGAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGTCGTACTCCACCACCATCCGCa                                                                                      <  2:215200/140‑1 (MQ=255)
                                                     tCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACTACTTTCCACGGTTCACCGTACCCCACCAACATcaccat                                                                                     >  1:857747/1‑138 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                   <  1:1107104/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACCTTTCGCCCTCCTCCACCACCACCCGccaga                                                                                   >  1:109186/1‑139 (MQ=255)
                                                         tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                                 >  2:618722/1‑143 (MQ=255)
                                                          gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                >  1:1000009/1‑143 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCCCCATCCGCCgggacacgc                                                                             >  2:841981/1‑136 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCACGCTCTGCCTCGCCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGAAGATGt                                                                             >  1:101657/1‑143 (MQ=255)
                                                               tGCGGTAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                           <  1:194834/143‑1 (MQ=255)
                                                                  ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACAGTGATTGATATGCAGGCTCTGCGTCGCCTGTTCTACTTTCCATGGTTCAGCATACTCCACTAACATCCTAAGGGAGATGTCGGaa                                                                        >  1:227417/1‑143 (MQ=255)
                                                                    tAAGAGTCGAGGTACGCGGCGAATTTTGTCGGTGACGCGGTGTTCTGCATTTCACCGTCATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                      <  2:335243/143‑1 (MQ=255)
                                                                                  cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                        <  1:1069897/143‑1 (MQ=255)
                                                                                  cccGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                        <  1:753753/143‑1 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTGCAGGGTTCGGCTTGCTCCACCACCATCCGCAGGAAGCTGTGGTAGGACTGTGCCGCTTTCAGtgc                                                    >  2:688483/1‑141 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACTTTCCACGCTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCGTTCAGTTCg                                                   >  2:574345/1‑143 (MQ=255)
                                                                                             ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAttg                                             >  1:191276/1‑141 (MQ=255)
                                                                                               ttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGCACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                           <  1:1047722/143‑1 (MQ=255)
                                                                                                tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                          <  1:618722/143‑1 (MQ=255)
                                                                                                   gggACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:197466/141‑1 (MQ=255)
                                                                                                    tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATATCCGCGATCTGATCTACTTAAACAGCAATCGTCACTTTTACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCATCGaacaaac                                      >  1:844367/1‑143 (MQ=255)
                                                                                                          ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTAGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  2:814149/143‑1 (MQ=255)
                                                                                                                cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                          >  2:924888/1‑143 (MQ=255)
                                                                                                                   cactgcacaGTGCTTGCTATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                       <  2:935149/138‑1 (MQ=255)
                                                                                                                          ccGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa                <  2:64153/143‑1 (MQ=255)
                                                                                                                            ttgattgatATTCGGGCTCTGTCTCTCCTCCTCTACTTTCCTCGGTTCGCCGTTCTCCACCCCCGTCCGCAGGAAGGTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACTCTCTTTGAAACGCCATGCCCAgg              <  2:776650/142‑1 (MQ=255)
                                                                                                                            gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATGCGCAGGAAAATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAACAAACACTCTTTCAAAAGTCATGcgccgc              >  2:693466/1‑138 (MQ=255)
                                                                                                                            gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGACACGCCATGCCCAgg              <  1:663772/143‑1 (MQ=255)
                                                                                                                             tgcttgctATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCTGGCGCAGGAAGATGTCGTCATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGt             <  2:931036/140‑1 (MQ=255)
                                                                                                                                 tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGACATGCCCAGGTCAAg         >  2:424627/1‑143 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  <  1:1032718/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925

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