Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,787 A→C 100% 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% 168.9 / ‑5.0 57I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (20/35);  minor base G (1/1);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑2132926
                                                                                                                                         |                                                                                                                                           
tggggTGGAGGGGAAACAGCGGAGATTGTGTGCGGGCGCGGTCGAAGTGGGCGGGGGGCTGGCGGGAAAAGTGGAAGTACCTGGAGAATTTTAACGCTGACTAAGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  2:496089/140‑1 (MQ=255)
ccgggTAGAGGCTACACACCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTCGTGCTTGCGGTAACAGTGGAAGTAGCCCGAGAATTTTTTCGTTGACGAGGTGATCGGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  1:269409/140‑1 (MQ=255)
  tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                      <  2:216864/143‑1 (MQ=255)
     tAGAGGGTAAACAGCTGGTGTTGTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:526153/143‑1 (MQ=255)
      agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTTATTGGTATGCGGGCTCTggcccgcc                                                                                                                                  >  2:587102/1‑136 (MQ=255)
        agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                <  2:980242/143‑1 (MQ=255)
            tAAACGGCTGATATTTTTTGCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttct                                                                                                                            <  2:57626/143‑1 (MQ=255)
              aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCccttctcc                                                                                                                          >  1:677594/1‑141 (MQ=255)
               aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACGGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTCCt                                                                                                                         >  1:856660/1‑143 (MQ=255)
                   cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc                                                                                                                     >  2:1047924/1‑143 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGATTTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGCTTGGTATGCAGGGTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  2:659963/1‑143 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTTATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcct                                                                                                                    >  2:979074/1‑142 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    <  1:669853/143‑1 (MQ=255)
                         ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCGCTACGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               >  1:461775/1‑143 (MQ=255)
                           tgatGCTGGCGCGGTCGAATTTGCCGTTGTGCGTGCGGTGGAGGTGGAGGTAGCCGGAGGATGTTTTCGGTGGCGAGGTGATCTGCATGTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCg                                                                                                             <  2:181485/140‑1 (MQ=255)
                            tttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                            <  2:181674/143‑1 (MQ=255)
                              tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                          <  2:426142/143‑1 (MQ=255)
                                ggcgcgcGGTCGCGGTTGTCGTTGGGCGCGGGGTGAAAGGGGAAGTGCCCGGAGAGTTGTTTCGGTGACGGGGTGGTCTGCATTTCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                        <  2:37615/140‑1 (MQ=255)
                                cTGGCGCGGTCCAAGTTGTCTTTGTGCTTGCGGTAGAAGTGGAAGTAGCCGGAGAGTTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                        <  2:956372/143‑1 (MQ=255)
                                   gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCAGCGTCGCCTCTTTTACTTTACACGGTTGACCGTACtg                                                                                                     >  2:271468/1‑142 (MQ=255)
                                            aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGCCGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGCTCCGCCGTCCCCCCCCCCCAt                                                                                            >  2:412951/1‑143 (MQ=255)
                                               ttGTCTTCGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                         <  1:341930/143‑1 (MQ=255)
                                                 gtgtttGTGCTTGCGGTAATGGTGGGGGTGCCCGGGGGATTTTGTCGGTGTGGATGTGATCTGCTTTTCTCCGTGATTGATGTGCAGGGTCTGCGTCGGCTCTTCTACTTTCCTCGGTTCGCCGTACTCCACCACCATCCGCa                                                                                       <  2:49029/140‑1 (MQ=255)
                                                  tgtttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  2:562920/141‑1 (MQ=255)
                                                  tCTTTGTGCTTGCGGTAAGAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  2:631246/143‑1 (MQ=255)
                                                    tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                    >  2:1103881/1‑143 (MQ=255)
                                                    tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACAATCCGCAGGa                                                                                    >  1:741301/1‑143 (MQ=255)
                                                     ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  1:260062/143‑1 (MQ=255)
                                                     ttGTGCTTGCGGTAAAAGTGGAAGAACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  1:172741/143‑1 (MQ=255)
                                                       ttgCTTCCAGTTAGAGTGGCAGTACCCGGATAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                 <  2:946639/142‑1 (MQ=255)
                                                        tgCTTGCGGTCACCGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGAGTTATATGCAGGCACTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                                <  1:567904/143‑1 (MQ=255)
                                                        tgCTTGCGGTACAAGTCGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                                <  1:771993/143‑1 (MQ=255)
                                                          cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACGGTTGGCAGTACGCCACCACCATCCGCAGGAAGATGt                                                                              >  2:452987/1‑143 (MQ=255)
                                                            tGCGGTAAAGGTGGCAGCACCCGGCGACTGTTTCCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  1:241123/143‑1 (MQ=255)
                                                                 tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGACTCTGCGTCGCCTCAGCTCCTTTCCCCGGGCCGCCGCACTCACCCACCCCCCCCCGTAAGAcgtcgtcctc                                                                       >  2:64117/1‑139 (MQ=255)
                                                                    gaGTGGACGTCCCCGGAGAATTTTTTCTGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                    <  1:771639/142‑1 (MQ=255)
                                                                      gttgAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  <  1:771983/140‑1 (MQ=255)
                                                                      gTGGAAGTACCCGGAGTTTTTTTTCGGGGACGAGGTGGTCTGCATTTCACCGTGATTGCTAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  <  2:1088693/143‑1 (MQ=255)
                                                                         gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                               <  1:406820/143‑1 (MQ=255)
                                                                           aGTCCCCGGAGAATTTTTGCTGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                             <  1:169476/143‑1 (MQ=255)
                                                                              aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGCACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                          <  1:868275/143‑1 (MQ=255)
                                                                               cccGGAGTATTTTTTCGGTGACGCGGTGATCTGCATTGCCCCGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                         <  2:877042/143‑1 (MQ=255)
                                                                                  ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  2:595412/143‑1 (MQ=255)
                                                                                    aattttttttCGGTGTCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                    <  1:554721/139‑1 (MQ=255)
                                                                                           tttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCGTTCAGTTCGGTATCgg                                             >  1:857780/1‑142 (MQ=255)
                                                                                                gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaac                                        >  1:860258/1‑142 (MQ=255)
                                                                                                gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                        <  2:256462/143‑1 (MQ=255)
                                                                                                   aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                                     >  1:139004/1‑143 (MQ=255)
                                                                                                       ggTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                 <  2:277019/143‑1 (MQ=255)
                                                                                                         tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTg                               >  1:1078612/1‑143 (MQ=255)
                                                                                                           aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCCGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCTCTTTGaa                             >  1:603026/1‑143 (MQ=255)
                                                                                                               gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTATCCACGGTTCGACGTCCTCAACCACCAACCGCCGGAAGATGTCGTAATCCTGTGACGCTTTCAGTTCGGACTAGAACAAACACTCTGTGgaaaac                         >  2:52294/1‑140 (MQ=255)
                                                                                                                cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGACTCTTCTACTTTACACGGTGCGCCGTACTCCACAACCATCGGCAGGAAGATGTCGTAATCCTGTGCCGCTTTACGTTCGACATCGAACAAACACTCTGTGAGAAGcc                        >  2:873715/1‑143 (MQ=255)
                                                                                                                  tttCACCGTGATTGATGTGCAGGCTCTGCGTCGTCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                      <  2:677594/143‑1 (MQ=255)
                                                                                                                   ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCATg                     >  2:516734/1‑143 (MQ=255)
                                                                                                                      aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc                  <  1:709807/143‑1 (MQ=255)
                                                                                                                                      aGGCTCTGCGTCGCTTCTTCTACTTTCCACGGTTCGCCGTACTCCACCGCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt  <  2:218446/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑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: