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% 156.2 / NA 57I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (26/28);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132649‑2132925
                                                                                                                                          |                                                                                                                                          
aTCTGGTAGAGGGTCAACAGCTGAGATTTTTTGCTGGCTCGGTGGAATTTGTCTTTGTGTTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                        <  2:317649/143‑1 (MQ=255)
aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACAGTGATTGATATGCAGGCTCTg                                                                                                                                        >  2:195362/1‑143 (MQ=255)
aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATCTCACCGTGATTCATACGCAGCCTCTg                                                                                                                                        >  2:258024/1‑143 (MQ=255)
   tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATGTGTCTTTGTGCTTTCGGTAAAAGGGGAAGTACCCGGAGAATTTTTTAAACGAAGCGAACATCAGTTTTTCACCGTGATTGCTACTCAGGGGCTGCGt                                                                                                                                     >  2:472764/1‑143 (MQ=255)
     gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                   >  2:266300/1‑143 (MQ=255)
       agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAACGTGGAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTTCATTTCACCGTTATTGATATGCCGTCTCTGCTGCGcc                                                                                                                                 >  1:15242/1‑143 (MQ=255)
        gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                >  1:132105/1‑143 (MQ=255)
          gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                              >  1:485989/1‑143 (MQ=255)
              aaaCAGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  1:921891/143‑1 (MQ=255)
                aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        >  1:527310/1‑143 (MQ=255)
                aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTACCGAGGTGGTCTGCATTTCACCGTGCTTGATATGCAGGCGCTGCGTCGCCTCCTCTGCt                                                                                                                        >  2:643960/1‑143 (MQ=255)
                aCAGCTGATATTTCTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        <  1:695546/143‑1 (MQ=255)
                aCAGCTGAGATTGTTTGCTGGCGCGGTGGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        <  2:653475/143‑1 (MQ=255)
                 cAGCTGATATTTTTTGCTGGCGCTGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                       <  2:612796/143‑1 (MQ=255)
                  cgctgctATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      <  1:693776/142‑1 (MQ=255)
                  aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      <  1:211205/143‑1 (MQ=255)
                       atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                 <  1:56046/143‑1 (MQ=255)
                          ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                              <  2:812783/143‑1 (MQ=255)
                          tttgttGCTGGTGGGTTCGTATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGGTACGAGGTGATTTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                              <  2:147191/143‑1 (MQ=255)
                               tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         <  2:527310/143‑1 (MQ=255)
                               tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         <  1:536818/143‑1 (MQ=255)
                                  tGGCGCAGTCGGAGTGGTCTTTGTGCTTGCGGTAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCCACTTTCCCCGGTTCGCAGtcc                                                                                                      >  2:945519/1‑141 (MQ=255)
                                    gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTACGCCGTaatg                                                                                                    >  2:296040/1‑140 (MQ=255)
                                                 tgccTTGGTGCTGGCGTTAAAGGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                       <  2:952866/140‑1 (MQ=255)
                                                 tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGCATGTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                       <  1:90590/143‑1 (MQ=255)
                                                  gTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                      <  2:731402/143‑1 (MQ=255)
                                                           cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                             <  1:350331/143‑1 (MQ=255)
                                                            ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCTCCACCATCCGCAGGAAGATgac                                                                            >  1:324431/1‑141 (MQ=255)
                                                                ggTAAAAGTGGAAGTATCCAGAGAAGTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGACTCTGCCTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCATGAAGATGTCGTaa                                                                        <  1:800463/143‑1 (MQ=255)
                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGACGTaa                                                                        >  2:709656/1‑143 (MQ=255)
                                                                  tagaaGTGGAAGTACCCGGAGAATTTTTTTTTTTTTTTTGTGGTCTGAATCTGACCGTGATTGTTCTGCAGGTTGTGCGTGGCCTCGTCTACTTTCCTCGGTTCCCCGTACTCCGCCAGCACCCGCGGGACGTTGTTGTTAtt                                                                      >  2:1029741/4‑142 (MQ=255)
                                                                      aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGCAAGATGTCGTAATCCtgt                                                                  >  1:403466/1‑143 (MQ=255)
                                                                       gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAATATTTCGGCATCCtgtg                                                                 >  2:419576/1‑143 (MQ=255)
                                                                          gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTCCCACGGTTCCCCCTACTCCACCCCCATCCGCCGGCAGATGTCGTAGTCCTGTGCCg                                                              >  1:243197/1‑143 (MQ=255)
                                                                            ggTACGCGGGGATGTTTTTCGGTGTCGGGGTGATCTGCATTTCACCGTGATTGATTTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  2:80862/142‑1 (MQ=255)
                                                                              tgcccGGAGAGTTTTGTCGGACACGGTGCGGTCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:511704/141‑1 (MQ=255)
                                                                              tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:132105/143‑1 (MQ=255)
                                                                              tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:485117/143‑1 (MQ=255)
                                                                                  cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  1:12408/75‑1 (MQ=255)
                                                                                  cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          >  2:12408/1‑75 (MQ=255)
                                                                                  cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCGCGACCATCCGCAGGCAGATGTCGGAATGCTGTGCCGCttccggt                                                      >  1:943653/1‑138 (MQ=255)
                                                                                     agaATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:329939/143‑1 (MQ=255)
                                                                                         ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCCGGAAGATGTCGTAATTCTGTGCCGCTTTCAGTTCGGTAt                                               >  1:1025515/1‑143 (MQ=255)
                                                                                         ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                               >  1:1009907/1‑143 (MQ=255)
                                                                                          tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                              <  1:616155/143‑1 (MQ=255)
                                                                                                 gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGCTCGGTATCGaccaca                                       >  2:386027/1‑141 (MQ=255)
                                                                                                       aGGTGATCTGCATTTCAACGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTACACGGTTCGACCTAATACACAACCATCCGCCCGAAGATGCAGTCATCCTATGACGATTTAAACTCGACAGCAAACAAacacaag                                 >  2:11200/1‑140 (MQ=255)
                                                                                                        ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                >  2:314875/1‑143 (MQ=255)
                                                                                                          tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTCTCAGTTCGGTATGGAACCAACACTCttat                              >  1:530548/1‑141 (MQ=255)
                                                                                                              cTGCGTTTCACCGTGATTGCTGTGCAGGCTCTGCGCCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                          <  2:681476/143‑1 (MQ=255)
                                                                                                              cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTGCACCCCCATCCGCAGGAAGCTGTCGTAATACCGTGCCGCTGTCAGTTAGGTAGCGCACACACGCTCTTTGaagc                          >  2:536080/1‑141 (MQ=255)
                                                                                                                 cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGGCGTAATCCTGTGCCGCTTGCAGTTGGGTAGCGAACAGACACTCGTTGaagggcc                       >  1:263500/1‑138 (MQ=255)
                                                                                                                     tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACACTCTTTGCAACGCCATGc                   >  1:1068457/1‑143 (MQ=255)
                                                                                                                         cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  1:266300/143‑1 (MQ=255)
                                                                                                                            gcttgatCTGGAGGCTCTGCCGCGCCTCTTCGACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTa            <  2:967291/141‑1 (MQ=255)
                                                                                                                                   atGCAGTCTCTGGCTCGCATGTTCTTGTTTCCTCGGTTCGCCGTACGCCACCAGCGGCCGCAGGAAGAGGACGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt     <  2:521130/143‑1 (MQ=255)
                                                                                                                                      cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTCCAGTTCGGTCTCGAACAAACACTCTTTGAAACCCCATGCCCAGGTAAAGCCCTgct  >  2:364097/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                          
ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132649‑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: ‑