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,7870AC94.5% 146.9 / NA 56I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (25/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑2132925
                                                                                                                                       |                                                                                                                                          
tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTTTTTGTTTTAGTACCCGGAGAATGTGTTCGGTGACGAGGGGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     <  1:63091/143‑1 (MQ=255)
 ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACAATAATTGATATACAGGCTCTACGTc                                                                                                                                    >  1:545640/1‑143 (MQ=255)
   tAGAGGGTAAACAGCTGATATTTTTTCTTGGAGCGGTTTAATTTTTCTTTGTGCTTGCTGTAAAAGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  1:216212/143‑1 (MQ=255)
   tAGAGGCTAAACAGCTGATATTTTTTGCTGGCGCGGTCTAATTTGTTTTTGTGCTTGCGGTTAAAGTGTAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  1:585622/143‑1 (MQ=255)
     gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAGAGTGGAAGTACCCGGAGGAGTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                <  2:483016/143‑1 (MQ=255)
         gTAAACAGCTTATATTTTTTGCTGTCGTTGTTGTTTTTGTCTTTTTGCTTGTGGTAAAAGTGGAAGTACCCGGAGAATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  1:1066256/143‑1 (MQ=255)
         gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTGAACGAAAAAAAAATGCAGGCTCTGCGACGCCtcttc                                                                                                                            >  1:503196/1‑143 (MQ=255)
         gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            >  1:539990/1‑143 (MQ=255)
             aCAGCTGCTATGTTTTGCTGGCGTGGTCGAATTTGTCTTTGTGCGTGCGGTATTAGGTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        <  2:1001765/143‑1 (MQ=255)
             aCAGCTGATCTTTTTTGCTGGCTCGGTCGAATTTGTCTTTTTGCTTGCGGTATTTTTGGTAGTACCCGGAGAATTTGTTCGGTGACGAGGGGATCGGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        <  1:50487/143‑1 (MQ=255)
               aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCTGTGACGATGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCTCTGCTTCTCCttt                                                                                                                      >  1:628638/1‑143 (MQ=255)
                gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                     <  1:415309/143‑1 (MQ=255)
                         ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAACTTTTTCGATGACGAGGTGATCTACATTTCACCGCGATTGATATGTAGGATCAGCCCCGCCTCTTCTACTTTCCACGGTTCg                                                                                                            >  1:1187596/1‑143 (MQ=255)
                            tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGTGATTGATATGCAGGATCTGCGTCGCCTATTCTACTTTCCACGGTTCGCCg                                                                                                         >  1:56503/1‑143 (MQ=255)
                               tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTCCGCCGtca                                                                                                      >  2:221219/1‑141 (MQ=255)
                                  cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                                   <  1:1001116/143‑1 (MQ=255)
                                   gcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGTAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                  >  1:803915/1‑143 (MQ=255)
                                       tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTGCTTTCCACGGTTCGCCGTACTccaccac                                                                                              >  2:578485/1‑143 (MQ=255)
                                           aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGTCGCGGTGATCTGCATTTCACCGTGATTGATATACAAGATCTGCGTCACCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                          >  1:40584/1‑143 (MQ=255)
                                                  tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTAGACGTACTTCACCATCATCCGccaga                                                                                   >  2:292678/1‑139 (MQ=255)
                                                   ttGTGCTTGCGGTAAAAGTAGAAGTACCCGGAGAATTTTTTCGGTGTCGTGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  1:101872/143‑1 (MQ=255)
                                                     gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                >  2:1038563/1‑143 (MQ=255)
                                                      tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCTTTTTACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                               <  1:463452/143‑1 (MQ=255)
                                                        cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                             <  1:432549/127‑1 (MQ=255)
                                                        cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                             >  2:432549/1‑127 (MQ=255)
                                                        cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGACGt                                                                             >  2:586624/1‑143 (MQ=255)
                                                               tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCTCTCCCCAGGAACATGTCGTAATc                                                                      >  2:233953/1‑143 (MQ=255)
                                                                   aGTGTAAGTACCCTGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  1:539007/143‑1 (MQ=255)
                                                                   aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:1185368/143‑1 (MQ=255)
                                                                     tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCCTTTGc                                                                >  2:1214653/1‑143 (MQ=255)
                                                                        aaGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                             >  2:212210/1‑143 (MQ=255)
                                                                          gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  >  2:370974/1‑120 (MQ=255)
                                                                          gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  1:370974/120‑1 (MQ=255)
                                                                           ttcccGGAGAGTTTTGTCGGTTACGGCGTGACCTGCATTTCACCGTGTTTGATATGCGGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:1247012/141‑1 (MQ=255)
                                                                                ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGGTGTCGTGATCCTGTGCCGCTTTCAGtt                                                     >  1:1209123/1‑143 (MQ=255)
                                                                                ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                     >  1:934785/1‑143 (MQ=255)
                                                                                ggAGAATTTTTGCGGTGCCGCTGTTATCTGCATTTCACCGTGTTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                     <  1:498452/143‑1 (MQ=255)
                                                                                 gtgaATTTTTTCGGTTACGATGTGATTTTCATTTCTCTGTGATTGATATTCTGGTTCTGCGTCGCCACTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:151107/141‑1 (MQ=255)
                                                                                   gagtttgttCGGTGACTAGGTGATCTGTATTTCATCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                  <  1:766267/140‑1 (MQ=255)
                                                                                      ttttttCGGTGACGAGGTGATTTGCATTTTACCTTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                               <  1:540009/143‑1 (MQ=255)
                                                                                       gttttCGGTGCCGAGGTGATCGGCATTTCCCCGTGATTGAGATGCAGGCTCTGCGTCGCGTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTTTGCCGCTTTCAGTTCGGTATc                                              <  2:933915/142‑1 (MQ=255)
                                                                                          ttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCAACGGATCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATGCTGTGCCGCTTTCAGTTCGGTATCGaa                                           >  1:1044739/1‑143 (MQ=255)
                                                                                             ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGTGTTGCCTCTTCTACTTTTCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                        <  1:933980/143‑1 (MQ=255)
                                                                                              gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:539990/143‑1 (MQ=255)
                                                                                                gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCTCCGTACTCCACCACCATCTGAAGGAAGATGTAGTAATCCTGTAAAGCAATCAGTTAGGTAGCGaacaaaca                                     >  1:170963/1‑143 (MQ=255)
                                                                                                gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGACGTCGTAATCCTGTGCCGCTTTCAGTTCGATATCGaacacacc                                     >  1:1243118/1‑142 (MQ=255)
                                                                                                   gatgTGATCTGCATTTCACCGTTATTGATATGCAGGCTCTGTGTCGTCTTTTCTTCTTTCTACGGTTCGCCGTGCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACActc                                  <  1:923788/140‑1 (MQ=255)
                                                                                                      gTGCTCTGCAGTTCACCGTCATTACTATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                               <  2:1152715/143‑1 (MQ=255)
                                                                                                      gTGATCTGCATTTCTACATGATTGAGATGCAGGCGCGGCGGGGCGTGTTGTAGTGGCAGCGGGTCGCCGGACTCGACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                               <  1:185871/143‑1 (MQ=255)
                                                                                                      gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCATTCAGTTCGGAATCGAACAAACACTAttt                               >  1:129890/1‑143 (MQ=255)
                                                                                                      ctggtCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCTGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                               <  2:722021/139‑1 (MQ=255)
                                                                                                       tGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACAGTGCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGTACTCCTGTGCCGCTTTCAGTTCGGTATCGCACAAACCCTCTTTg                              >  2:555094/1‑143 (MQ=255)
                                                                                                            tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTGTTTGaaaag                         >  1:547825/1‑141 (MQ=255)
                                                                                                             gCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTCCCGCTTTCAGTTCGGTATCGAACAAACACCATTTGAAAcac                        >  2:44045/1‑141 (MQ=255)
                                                                                                                     ccGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTTCTTTCCACGGTTCGCCGTACTCCATCATCTTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa                <  1:571000/143‑1 (MQ=255)
                                                                                                                                   cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACGCTTTAAAACGCCATGCCCAGGCAAAGACCAGGt  >  2:1173228/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                          
TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑2132925

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