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,7870AC95.8% 139.4 / ‑5.9 49I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (25/21);  total (27/22)
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.98e-01

GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132648‑2132924
                                                                                                                                           |                                                                                                                                         
gacttGGTAGAGGGTCAACAGCTGTTATTTTTTTCTGGCGCGGTCGAATTTGTCTTGGTGCTTGCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                        <  2:154423/139‑1 (MQ=255)
     ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCGCTGCGTc                                                                                                                                   >  1:710554/1‑143 (MQ=255)
      gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCg                                                                                                                                  >  1:1027280/1‑143 (MQ=255)
       tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                 <  2:109435/143‑1 (MQ=255)
         gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTGGCct                                                                                                                               >  2:542826/1‑143 (MQ=255)
            ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGTTCTGCGTCGCCtctt                                                                                                                            >  1:823508/1‑143 (MQ=255)
            ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                            <  1:789264/143‑1 (MQ=255)
               gaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                         <  2:882354/142‑1 (MQ=255)
               aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCAGGAGAATTTTTTCGGTGACGAGGTGATCTGCATATCACCGTGATTGATCTGCAGGCTCTGCGTCGCCtcgtctg                                                                                                                         >  2:911742/1‑142 (MQ=255)
                  cAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGTGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGCGATTCATATGCAGGCGCTGCGTCGCCTCTTCTCCtt                                                                                                                      >  2:620778/1‑143 (MQ=255)
                    gCTGATCTGTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                    <  1:168633/143‑1 (MQ=255)
                         tattcgttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTACGAGTCGAACAACCAGAAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                               <  1:952850/137‑1 (MQ=255)
                           ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGAtt                                                                                                             >  2:710552/1‑143 (MQ=255)
                               ttGCTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                         <  1:523197/143‑1 (MQ=255)
                                   tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGttc                                                                                                     >  1:949141/1‑141 (MQ=255)
                                        cgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                <  1:542678/143‑1 (MQ=255)
                                           tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGGCCGCCGGAC‑ccacccac                                                                                            >  2:705704/1‑140 (MQ=255)
                                             gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGTGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTTCGTCGCCTCTTCAACCTTACTCGGTTCGACGTAACcccccaccc                                                                                           >  2:932425/1‑142 (MQ=255)
                                                        tgcgCTTGCGGTAACCGTGGAAGTACCCGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                                <  1:877601/140‑1 (MQ=255)
                                                          ttcTTGCGTTACACGCGGAAGCACCCGGAGAATTTTTTCGGTGCCGAGGTGACCTGCATTTCACCGTGATCGATATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                              <  2:242694/141‑1 (MQ=255)
                                                            cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGCCCGCCGTACTCCCCCACCATCCGCTggccgccgt                                                                            >  1:724687/1‑136 (MQ=255)
                                                            cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGt                                                                            >  1:421630/1‑143 (MQ=255)
                                                                 gggaaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCCACTTTCCACGGCTCGCCGTACCCCACCACCATCCGCAGGAACCTGACGTaa                                                                       >  2:1056218/4‑143 (MQ=255)
                                                                 ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGACATGCAGGCTGGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTaa                                                                       >  1:149629/1‑143 (MQ=255)
                                                                  gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGAGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCGGCAGGAAGATGTCGTAAt                                                                      >  1:641497/1‑143 (MQ=255)
                                                                   tacgaGTGGAAGTACCCGGAGATGTTGTGCGGTGACGGGGGGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                     <  2:192509/139‑1 (MQ=255)
                                                                         tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGACGTAATCCTgcgc                                                               >  2:1010090/1‑143 (MQ=255)
                                                                           gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGCCATTGCGGACTCCACCACCATCCTCCGGAAGATGGCGACATCCTGTGCCg                                                             >  2:422410/1‑143 (MQ=255)
                                                                                ttgcGGAGTGTTTTTTCGGTGACGACGTTATTTTCATTTCACTGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                        <  2:510325/140‑1 (MQ=255)
                                                                                 cccGGAGAATTTTTTCGGTGACGAGGTGCTCGTCATTTCCCCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCGTTCa                                                       <  1:325240/143‑1 (MQ=255)
                                                                                 cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTCGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       <  1:668118/143‑1 (MQ=255)
                                                                                     gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                   >  1:679182/1‑143 (MQ=255)
                                                                                     gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCAGCTTTccggtc                                                   >  1:162222/1‑138 (MQ=255)
                                                                                              ttCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCTCGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                          <  1:491043/143‑1 (MQ=255)
                                                                                               tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCAGTACTACACAACCAACAGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                         >  1:257344/1‑143 (MQ=255)
                                                                                                 ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                       <  2:589426/143‑1 (MQ=255)
                                                                                                 ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGGAACCCTGTGCCGCTATCACGTCGGGATGGaacaa                                       >  2:893165/1‑143 (MQ=255)
                                                                                                  gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                      <  1:929114/143‑1 (MQ=255)
                                                                                                   tGACGGGGTGATCTGCCGTTCACCGTGACTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGTAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaac                                     <  2:690032/143‑1 (MQ=255)
                                                                                                   tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCATCGaacaaac                                     >  1:21933/1‑143 (MQ=255)
                                                                                                     aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAacacacac                                   >  1:644158/1‑143 (MQ=255)
                                                                                                        aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACActct                                <  2:256380/143‑1 (MQ=255)
                                                                                                         ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                               <  1:985386/143‑1 (MQ=255)
                                                                                                           tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTCTCAGCACGGTCTCGAACAACCACCCTTTg                             >  2:41718/1‑143 (MQ=255)
                                                                                                             aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                           >  1:1038174/1‑143 (MQ=255)
                                                                                                                 gCACTTCACCGTGACTGATATGCAGGATCTACCTCCCCTCTCCCACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                       <  2:512430/143‑1 (MQ=255)
                                                                                                                  cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAGCACTCTTTGAAACGcc                      >  1:384456/1‑143 (MQ=255)
                                                                                                                            ggattgctATGCAGGCTCTGCGGCGCCTCTTGTGCCGTCCACGGTTCGCCGCACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGt            <  2:1034716/142‑1 (MQ=255)
                                                                                                                               ctgagacgCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTACACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa         <  1:118429/136‑1 (MQ=255)
                                                                                                                                      gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg  >  1:1070936/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                         
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132648‑2132924

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