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,7870AC97.9% 137.5 / NA 48I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (17/30);  minor base G (1/0);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925
                                                                                                                                            |                                                                                                                                          
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATCGATATGCAGGagc                                                                                                                                          >  2:218385/1‑141 (MQ=255)
     tggtggAGGGTAGAGAGCTGATATTTTTTGCTGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     <  2:849192/143‑1 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGGGGTGGTTTGTGTGTTTCCTTGATTGATATGCAGGCTCTGCGTc                                                                                                                                    <  2:551541/143‑1 (MQ=255)
       gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATATGTCTTTGTTCTTGCGGTAAAAGTGGAAGTACCCGGAAAATTTTTTCGGTGGCGAGGTGATATGCTGTTGATCGTGATTGTTATGCAGGGTccgtggcg                                                                                                                                   >  2:507725/1‑136 (MQ=255)
        gagagGGTAAGCAGCTGATATTTTTTGCTGGCGCGGTCTAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  2:536796/142‑1 (MQ=255)
          gagGGTCGACAGCTGGTGTGTTTTGCTGGCGCGGTCGAATTTGTTTTTGTGCTTGCGGTAGAAGTGGAAGTACCCGTGGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGGGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                <  2:268696/143‑1 (MQ=255)
             ggTAAACAGCTGATATTCTTTGCTCGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTCGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                             <  1:503980/143‑1 (MQ=255)
                  aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTCGTCTTTTTTCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        >  1:686581/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTCCGTCGCCGCTTCTACTTTccg                                                                                                                   >  2:949234/1‑142 (MQ=255)
                                   cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGGCTCGCCGTa                                                                                                       >  2:620867/1‑143 (MQ=255)
                                    tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                      >  1:324085/1‑143 (MQ=255)
                                    gtgcgcGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                      <  1:338776/141‑1 (MQ=255)
                                      gcgcGGTCGCGTTTGTCTTTGTGTTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                    <  2:710559/143‑1 (MQ=255)
                                          gggcGAATTTGTCGTTGTGCTGGCGGTGAAAGTGGAAGTGCGCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                <  2:21165/140‑1 (MQ=255)
                                            tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                              >  1:700112/1‑143 (MQ=255)
                                                       tgtgtgGTGGCGGGAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGa                                                                                   <  1:653595/141‑1 (MQ=255)
                                                        gtgtgCTGGCGCTAAAAGTGGAAGTACCCGGAGAGTTTTTTCGGTGACGCGGTGCTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  1:125564/142‑1 (MQ=255)
                                                            gCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGAAAGATg                                                                              >  2:487750/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            >  1:113224/1‑143 (MQ=255)
                                                              tggCGGTAAAAGTGGAAGTACCCGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  1:346498/141‑1 (MQ=255)
                                                                    tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGGTGTCGTAGTc                                                                      >  1:490795/1‑143 (MQ=255)
                                                                        cgtgtgaGTACCCGGAGAATTTTTGCGGCCACGAGGTCATCTGCACTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  1:60380/137‑1 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCGTGTGc                                                                >  1:94892/1‑143 (MQ=255)
                                                                            gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCg                                                              >  1:630469/1‑143 (MQ=255)
                                                                            gAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                              <  2:889668/143‑1 (MQ=255)
                                                                              agcaCCCGGATTAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  2:909492/140‑1 (MQ=255)
                                                                                tgcccGGAGAGTGTTTTCGGTGACGGGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:898530/141‑1 (MQ=255)
                                                                                tACCCGGAGCATTCTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:734257/143‑1 (MQ=255)
                                                                                   ccGGAGAATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCCCCACCACCCTCCGCATGACGATGTCGTAGTCCGGTGCCGCCGTCAg                                                       >  2:829734/1‑143 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGAGATTGATAGGCAGGCTCTGCGTCTCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCGCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:370561/143‑1 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  2:785012/143‑1 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:697699/143‑1 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  2:942479/143‑1 (MQ=255)
                                                                                        gaATTGCTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                  <  1:967974/143‑1 (MQ=255)
                                                                                         aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTAATCGTGTGCCGCTTTCAGTTCGGt                                                 >  1:104221/1‑143 (MQ=255)
                                                                                                 cggcgccGAGGTGCTGTGCATTTCATCGTGACTGGTATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCTCCGTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                         <  2:459982/137‑1 (MQ=255)
                                                                                                     gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacacaca                                     >  1:781206/1‑143 (MQ=255)
                                                                                                     gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGATACCGaactgaaa                                     <  1:925601/143‑6 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGTAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACAatc                                  >  1:43814/1‑143 (MQ=255)
                                                                                                          ggTGATCTGCATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  2:792734/143‑1 (MQ=255)
                                                                                                              aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAGACACTCTTTGaa                            >  1:617868/1‑143 (MQ=255)
                                                                                                                     tctccgcGTGAGTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCATGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCTGTATCGAACAAACACTCTTTGAAACGCCAt                     <  2:904797/137‑1 (MQ=255)
                                                                                                                      ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  2:676119/143‑1 (MQ=255)
                                                                                                                      ttCACCGTAATTGATATGCTGTCTCTGCGTCGCTTCTTCTACTTTGCACGGTTCGCCGTACTTCACCACCATCCGCAGGAAGATGTCGTAATTCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTTTTTGAAACGCCATg                    <  2:555969/143‑1 (MQ=255)
                                                                                                                           cGTGATTGGTCGGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  2:687009/143‑1 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               >  1:646974/1‑143 (MQ=255)
                                                                                                                                   atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc       <  1:316409/143‑1 (MQ=255)
                                                                                                                                      tGCAGGCGCTGCGTCGCTTCGTCTCCTTTCCACGGTTCGACGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTGTAAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg    <  2:710516/143‑1 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCGGCGCCACCTGCGATGCGGAGTCGAACAAACACGCTTTGAAACGCCACGCCCAGGTACAgccgaggt  >  2:865998/1‑136 (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: