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.7% 132.8 / NA 48I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (27/18);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑2132925
                                                                                                                                       |                                                                                                                                          
tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGACAATTTTTTCGGTGACGAGGTGATCTACATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     >  2:572483/1‑143 (MQ=255)
    agagCGTAAACCGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTCGTGCTTGCGGTACAAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                                 <  1:702421/143‑1 (MQ=255)
        ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCgggtctg                                                                                                                             >  1:421625/1‑137 (MQ=255)
        ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                             >  1:1021553/1‑143 (MQ=255)
         gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  1:181422/143‑1 (MQ=255)
              cAGCTGATATGTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                       <  2:676696/143‑1 (MQ=255)
                     tatGTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTCGCGGTAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                                <  1:113430/143‑1 (MQ=255)
                         ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                            >  1:912043/1‑143 (MQ=255)
                           ttGGTGGCGCGGTCGAATTTGGCTTTGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                          <  1:816426/143‑1 (MQ=255)
                              cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       >  1:333526/1‑143 (MQ=255)
                                ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGGTCGCCGTAcg                                                                                                     >  2:418921/1‑142 (MQ=255)
                                  cgcgGTCCAATTTGTCTCTGTGCGTGCGGTAAAAGTGGAAGTGCCCGGAGAGTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCCGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                                   <  1:1110681/143‑1 (MQ=255)
                                      gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCCccacca                                                                                               >  2:219989/1‑143 (MQ=255)
                                      gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccccca                                                                                               >  1:178386/1‑143 (MQ=255)
                                       tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTtcaccac                                                                                              >  1:1068268/1‑143 (MQ=255)
                                       tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                              <  1:876610/143‑1 (MQ=255)
                                       tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGATTCGACGTACGccgccac                                                                                              >  2:881234/1‑143 (MQ=255)
                                         gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTTTCCACGGTTCGCCCTACTccaccacca                                                                                            >  2:184762/1‑143 (MQ=255)
                                              tgtccgggtgGTTGCGGTGAAAGTGTGAGGACCCGGAGTATGTTTTCGGTGACGAGGTGATCTGTATTTCACCGTGTTTGTTATGCCGGCTCCGCGGCGCCGCTTCTACTTTCCGCGGTGCGCCGTACTCCACCACCATCCGc                                                                                       <  2:157697/136‑1 (MQ=255)
                                                 cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCCgg                                                                                    >  1:1050076/1‑143 (MQ=255)
                                                 cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAgg                                                                                    >  1:819590/1‑143 (MQ=255)
                                                  tttGTGCTGGCGGTAAGCGTGGAAGTACTCGGCGTTTTTTTTCGGTGACGCGGTGATCTGCATTTCACCGTGATTGATATGCAGCCTCTGCGCCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCTGCAGGa                                                                                   <  2:1026619/143‑1 (MQ=255)
                                                    tgtgCTTGCGGTAAAAGTGGAAGTAACCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATATCAAAGTGATTGATATGAAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCCCACGAAg                                                                                 >  1:824930/1‑143 (MQ=255)
                                                         atGCGGTAAAAGCGGAAGTCCCCGGACAACTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  1:375497/142‑1 (MQ=255)
                                                           gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                          <  2:626476/143‑1 (MQ=255)
                                                            cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCCTCCGCAGGAAGAAGTCgaa                                                                         >  2:643802/1‑141 (MQ=255)
                                                               tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGACCTGCGTCGCCTATTCTACTTTCCACGGTACGCCGTACTCCACCAACATCCGCAGGAAGATGTCGTAGTc                                                                      >  2:436581/1‑143 (MQ=255)
                                                                aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGCACTCCCCCACCATCCGCAGGCAGCTGGCGCCATcc                                                                     >  2:611073/1‑143 (MQ=255)
                                                                       gtaGTACCCGGAGGAGTGTGTCGGTGACGAGGTGATCGGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                              <  1:886787/141‑1 (MQ=255)
                                                                         aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCCTCCGCCGGCAGATGTCGTAATCCTGTGCCGCt                                                            >  2:125425/1‑143 (MQ=255)
                                                                           tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATACTGTGCCGCttt                                                          >  1:1158967/1‑143 (MQ=255)
                                                                           tACCCGGAGAATCTTTTCGGTGACGAGGTGATCGGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  1:854406/143‑1 (MQ=255)
                                                                                ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCATGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                     >  1:772621/1‑143 (MQ=255)
                                                                                ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCTTAATTCTGTACAGCTTTCACtt                                                     >  1:892908/1‑143 (MQ=255)
                                                                                  agaATGCTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  2:939484/143‑1 (MQ=255)
                                                                                      ttttttCGGTGACGAGGTGATCTGCATTTCACCGAGATTGAAATAAAAAAACTGCGTCGCCTCTTCTACTTTCTAAGGTTCGCCGGACTCCACCACCATCCACAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                               >  1:1032454/1‑143 (MQ=255)
                                                                                       tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCGCCAGCATCCGCAGGAAGATGTCGTAATCCTGTGCGGCTTTCGGTGGGGTATc                                              >  1:1033742/1‑143 (MQ=255)
                                                                                                         aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCTTCTTCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                            >  1:963064/1‑143 (MQ=255)
                                                                                                         aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCTCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                            <  1:756650/143‑1 (MQ=255)
                                                                                                            tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTTGCGCTTTCAATTCGGTATAGAACAAACAATCTATGAAACg                         >  1:496684/1‑143 (MQ=255)
                                                                                                              cATTTCACCGTGATTGATATGCAGACTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAAAAAACACTCTTTGAAACGcc                       >  1:505207/1‑143 (MQ=255)
                                                                                                                tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCATGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGTCac                     <  2:772621/143‑2 (MQ=255)
                                                                                                                tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCACTTCGGTATCGAACAAACACTCTTTGAACCGCCAt                     >  2:846338/1‑143 (MQ=255)
                                                                                                                 ttCATCTTGATTGATGTGCGGGGTCTGCGTCGCCTCTTCTGCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACTCTTTGAAACGCCATg                    <  2:418345/143‑1 (MQ=255)
                                                                                                                       ctcattgctATGCAGGCCCTGCGTCGCCTCTTCGACTTTCCACGGTTGTCGGTACGCCACCAACATCCGCAGGAAGATGTCGTAGTCCTGTGCCGCTGTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAgg              <  2:243003/140‑1 (MQ=255)
                                                                                                                           ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa          <  1:1145368/143‑1 (MQ=255)
                                                                                                                            tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAg         >  2:330003/1‑143 (MQ=255)
                                                                                                                                   cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  >  2:408441/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                          
TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑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: