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,7870AC96.0% 67.1 / ‑5.4 26I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (10/14);  minor base T (0/1);  total (11/15)
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.70e-01

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132647‑2132920
                                                                                                                                            |                                                                                                                                     
ggATCTGGGAGGGGGTGAACAGCTGATAGTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                     <  1:156723/143‑1 (MQ=255)
  aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGGGAGCTGCACTTCACCGAGGTTGATGTGCAGGCTCTg                                                                                                                                   >  2:283301/1‑143 (MQ=255)
             ggcgaccAGCTGATATGTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                        <  1:48406/137‑1 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTCATATGCAGGATATGAATCACCtctt                                                                                                                        >  2:106008/1‑143 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                       <  1:278848/143‑1 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                       >  2:278848/1‑143 (MQ=255)
              ctaacgaGCTGATAATTTGTGCTGGCGCGGTCGCATTTGTATTTGTGCTTGCGGTAATAGTGGAAGTACGCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                       <  1:151013/137‑1 (MQ=255)
                     tcTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                <  2:294190/142‑1 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCg                                                                                                            >  2:706175/1‑143 (MQ=255)
                                ttacgtgcgcgcTCGAAGTGTTCTTTGTGCTTGCTGTAAATGGGGAAGTAACCTTAGAATTTTTTTGGTGACGAGGTGATCTGCATTTTACCGTGATTGATATGCAGGCTCTGTGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                     <  2:308691/137‑1 (MQ=255)
                                           gtccattttGTCTTTGTTCTTGCCGTAAAACTGCAAGTACCCGGAGAATTTTTTCGGTGACGCGGTGATCTGTATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                          <  2:391349/137‑1 (MQ=255)
                                             cGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                        <  2:469587/143‑1 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTGTACTTTCAACGGTTCGCCGTACTccccccccg                                                                                       >  1:336525/1‑142 (MQ=255)
                                                   tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                  <  1:618389/143‑1 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTTCCCGGAGTATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                       <  1:376907/143‑1 (MQ=255)
                                                                        aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGGTCTGCATTTCACCGTGATTGATATGTAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTCCGCCGTACGCCACCACCCTCCCCAGGAAGACGTCGTAATCAtgt                                                             >  2:77958/1‑143 (MQ=255)
                                                                          tgggtgTACGCGGAGATTTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGTTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCGACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                           <  2:632622/138‑1 (MQ=255)
                                                                                  cccGGAGACTTTTTTCCGTGACGACGTGATCTGCATTTCACCGTGATTGATATCCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                   <  1:408842/143‑1 (MQ=255)
                                                                                           ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAGTTTCCACGGTTCGCCGTACTCCACCACCATCCGCTGGAAGATGTCGTACTGCTGTGACGCGTTCAGTAAGGTaa                                          >  1:7557/1‑142 (MQ=255)
                                                                                              tttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACTTTCTCCGGTTCGCCGTACTCCACGCCCATCCGCAGACCGATGTCGGAATCCTGTGCCGCTTCCAGTCCGGTATCGa                                       >  2:110039/1‑143 (MQ=255)
                                                                                                tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATCCACCACCATCCGCAGAAAGACCGATGAAACCTGTGCCGCTTTCAGTTTGGTATCGcac                                     >  2:342315/1‑143 (MQ=255)
                                                                                                             caTCTGCATTTCTCCGTTATTGATATGCATGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCTCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGa                        <  1:204284/142‑1 (MQ=255)
                                                                                                               tCTGCTTTTCACCGTGATTGATATGCAGGTTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                      <  2:18616/143‑1 (MQ=255)
                                                                                                                      ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGGATCGAACACACACACGTTGAAAAGCCATg               >  2:149729/1‑143 (MQ=255)
                                                                                                                             tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTCCGGTATCGAACAAACACACTTTGAAACGCCATGCCCAGGt        >  2:740864/1‑143 (MQ=255)
                                                                                                                                   cggtGCGGGCTCTGGGGCGGCTCTTCTACTTTCCACGGTGCGCCGTGCGCCCGCGCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc  <  2:593798/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                     
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132647‑2132920

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