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.2% 159.4 / NA 57I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/30);  minor base T (0/1);  total (26/31)
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.90e-01

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926
                                                                                                                                            |                                                                                                                                            
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGGCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCAGGCTc                                                                                                                                            >  2:255127/1‑143 (MQ=255)
       gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAGTACGGCACTTGAGAGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                     >  1:874317/1‑143 (MQ=255)
        tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                    >  1:916745/1‑143 (MQ=255)
             ggTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTATTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                               >  2:488692/1‑143 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCctcct                                                                                                                               >  1:529572/1‑141 (MQ=255)
                  aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTATTTGTGATTGCGGTAAAAGTGGAAGTACACGGAGAATTTTTTCGGTGACGAGGTGATATGCATTTCAACGTGATTGATATGCAGGCTCTGCGTAGCCTGTTCTaat                                                                                                                          >  2:439031/1‑141 (MQ=255)
                       ggATATTTTTTGCTGGCGCGGTCGCATTTGTGTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                     <  2:412778/142‑1 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGGTTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                   >  1:829867/1‑143 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTCGCGGTAAAAGTGGAAGTACCCGGAGCATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCTCTGCGTCGCCGCTTCTACTCTCCACg                                                                                                                   >  1:473610/1‑143 (MQ=255)
                         atatCTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtttgcccg                                                                                                                   >  2:327014/1‑138 (MQ=255)
                            tttgttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  1:733861/143‑1 (MQ=255)
                              ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                             <  1:852193/143‑1 (MQ=255)
                                ttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                           <  2:179005/143‑1 (MQ=255)
                                 ggCTGGGGCGGAGGGAGTTGGCTATGGGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                          <  1:13253/142‑1 (MQ=255)
                                  gCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGt                                                                                                         <  1:790906/143‑1 (MQ=255)
                                    tGGCGCGGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                       <  2:803385/143‑1 (MQ=255)
                                       agctgTCGTATTTGTCTGTGTCGTTGCGGTAAAGGGGGCAGTACCTGGAGAATTTTTTCGTTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACGGTTCGCCGTACTcc                                                                                                    <  2:30798/139‑1 (MQ=255)
                                         cgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGtcctcccc                                                                                                  >  1:963571/1‑141 (MQ=255)
                                           gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACACGGTTCGCCGTACTccacca                                                                                                >  2:602701/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacccccc                                                                                             >  1:554298/1‑142 (MQ=255)
                                                aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAAGTTTTTCGGGGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGGGTGGCCTCTTCTACTTTCCACGGTTCGGCGCACTCCGCCCCCATc                                                                                           >  2:620790/1‑143 (MQ=255)
                                                     tCTTTGTGATTCCTGTTAAAGTGGAAGTACCCGGAGAATTTTTTCGTTGACGATGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  2:107336/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGCTCGCCGTACTCCACCACCATCCGCAGGa                                                                                    >  1:118145/1‑143 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCTACTTTCCACGGTTCGCCGTACTCCACCACCGTCCGCAGGa                                                                                    >  1:60618/1‑143 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGAATGTTCTACTTGCCACGGTTCGCCTCACTCCACCACCATCCGCAGGa                                                                                    >  1:1188795/1‑143 (MQ=255)
                                                            gcatGCGTTTAAAGTTGAATCACCCCGAAAATGTTTTCGCTTACGATGTGATCGTCATCTCACCGTGATTGATATGCAGGCTCTTCGTCGCCTCTGCCACTTTCCACACTTCCCCTTACTCCACCACCATCCTCAGTAAGATg                                                                               <  1:829329/140‑1 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                             >  1:985495/1‑143 (MQ=255)
                                                              ttGCGGGAAGAGTGGAATTACCCGGAGAGTTTTCTCGATTACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                             <  2:526022/143‑1 (MQ=255)
                                                                 cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTa                                                                          <  1:240034/143‑1 (MQ=255)
                                                                  ggtacgaGTGGATCTGCCCGCAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         <  2:1108299/137‑1 (MQ=255)
                                                                  ggTAAAAGGGGAAGTACCCGGAGAATGGTTTCGGTGACGAGGTGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         <  1:655313/143‑1 (MQ=255)
                                                                   gTAGAAGGGGAAGTACCCGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                        <  1:494456/143‑1 (MQ=255)
                                                                      aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAGCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACACGGTTCGCCGTACTCCACCACCATCCGCCGGCAGATGTCGTAATCCt                                                                     >  1:149008/1‑143 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                                 >  1:1107268/1‑143 (MQ=255)
                                                                                tgcccGGAGAACTTTCTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACCCCTTTCCGTGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:327014/141‑1 (MQ=255)
                                                                                tACCCGGCGTCTTTTTTCGGTGACGCGGTTATCTGGATTTTACCGTGATTGATATGTAGGTTCTGCGTTTCCTGTTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCGTAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  2:79173/143‑1 (MQ=255)
                                                                                tACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:592409/143‑1 (MQ=255)
                                                                                 gcccGGAGAATGTTCTCGGTGACGGGCGGATCTGCATTTCGCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                          <  1:1128532/142‑1 (MQ=255)
                                                                                  ccccGAGAATTGTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGCGCCGCTTTCa                                                         <  2:918296/143‑1 (MQ=255)
                                                                                  cccGTAGAATTTTGTCGGTGACGAGGTGATCTGCATGTCACCGTGAGTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                         <  2:1135059/143‑1 (MQ=255)
                                                                                     gcagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTTCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:248334/141‑1 (MQ=255)
                                                                                       agagTTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGAAGTAAGATATCGTAATCCTGTGCCGCTTTCAGTTCg                                                    <  2:662652/143‑1 (MQ=255)
                                                                                           ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                                >  2:40894/1‑143 (MQ=255)
                                                                                               ttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                            <  1:426460/143‑1 (MQ=255)
                                                                                               ttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                            <  1:810881/143‑1 (MQ=255)
                                                                                                          tgTGATCTGCATTTCACCGTGTTTGACATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                 <  2:165000/142‑1 (MQ=255)
                                                                                                          ggtggtCTGCATTCCACCGTGTATGACACGCAGGCTCTCCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGCCGTAATCCTGTGCCGCTTTCAGTTCGGTAGCGAACAAACACTCtt                                 <  2:1013325/143‑1 (MQ=255)
                                                                                                          ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                 >  1:947067/1‑143 (MQ=255)
                                                                                                          ggTGATCTGCATTTCACCGTGAATGAAAAGCAGGCTCTGAGTCCACTCACGAATGTCTGTATTGTTGCCATACACTATCACCAACCGCAGGAAGAAGTCGTAATCCGGTGCCGCTTTCAGTTCGGTATCGaataaaagcgctt                                 >  2:252945/1‑136 (MQ=255)
                                                                                                           gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGCACAAACACTCttt                                >  1:280391/1‑143 (MQ=255)
                                                                                                                cTGCATTCCACCCTGGTTGTTTTGCAGGCTCTTCCTCGCCTCTTCTACTTTCTACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                           <  2:329146/143‑1 (MQ=255)
                                                                                                                cTGCAGTCCCCCGTGATTGAGATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGGTCGCCGTACTCCACCACCATCCGCATGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                           <  2:1067829/143‑1 (MQ=255)
                                                                                                                         aCCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGCTCGGCACCGAACCCACCCTCTTTGAACCGCCCTGccc                  >  2:168225/1‑143 (MQ=255)
                                                                                                                          ccGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACGAACACTCTTTGAAACGCCGTGcaca                 >  2:633020/1‑143 (MQ=255)
                                                                                                                                  gatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa         <  1:1045966/143‑1 (MQ=255)
                                                                                                                                    tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTTCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCGTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc       <  1:389855/143‑1 (MQ=255)
                                                                                                                                         aGGCTCTGCGTCGCCTCCTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGGAAAGACCAGGtc  >  1:690354/1‑142 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926

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: