Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,180 (T)6→7 coding (150/1047 nt) wzzE → enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,1741.T98.0% 158.7 / ‑3.6 49coding (144/1047 nt)wzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (21/27);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTC  >  NC_000913/3969035‑3969299
                                                                                                                                            |                                                                                                                              
cacaACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACg                                                                                                                                 >  1:13369/1‑140 (MQ=255)
     ccAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGttttt                                                                                                                          >  2:155999/1‑140 (MQ=255)
         tGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGTTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGATAATCGAGCTGGCATATACGtttttttgc                                                                                                                      >  1:290648/1‑137 (MQ=255)
         tGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGtttttttgc                                                                                                                      >  1:690065/1‑137 (MQ=255)
             tGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGCTTGAGGTCATCGCGCTGGCGTGAACGTTTTGTTGCTCGt                                                                                                                   >  2:418210/1‑140 (MQ=255)
               ggAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCa                                                                                                                 >  2:473661/1‑140 (MQ=255)
               ggAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGCCTGGCGTTCGCGGTAAACGCGCTGGAGTATACGTTTTTTTGCTCGTCa                                                                                                                 >  2:347045/1‑140 (MQ=255)
                   aCCGGCCGAAGACGCAGAAAATGAAATGGATATTCGTAGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAg                                                                                                             <  1:347045/140‑1 (MQ=255)
                           aaGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTc                                                                                                     >  1:549469/1‑140 (MQ=255)
                            aGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCg                                                                                                    <  1:155999/140‑1 (MQ=255)
                              aCGCTGAAAATGAACTGCATATGCGTGGGTTGTTTCGTATCTTGTGGGGTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGAc                                                                                                  <  2:290648/140‑1 (MQ=255)
                               cGCTGAAAATGAACTGGTTATTAGTGGGTTGTTTCTTACCTTGTGGGCAGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGGGCTCGACg                                                                                                 <  2:263636/140‑1 (MQ=255)
                                gCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACgg                                                                                                <  2:549469/140‑1 (MQ=255)
                                    aaaaTGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGAt                                                                                            <  1:570122/140‑1 (MQ=255)
                                          aaCTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGa                                                                                      <  1:825529/140‑1 (MQ=255)
                                             tGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCg                                                                                   <  1:473661/140‑1 (MQ=255)
                                               gATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGa                                                                                      >  2:705670/1‑135 (MQ=255)
                                               gATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGa                                                                                      <  1:705670/135‑1 (MQ=255)
                                                 tatTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCa                                                                               <  1:485413/140‑1 (MQ=255)
                                                 tatTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCa                                                                               <  2:798252/140‑1 (MQ=255)
                                                  atTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCaa                                                                              <  2:475176/140‑1 (MQ=255)
                                                  atTCGTGGGTTGTTTCGAACCTTGTGGGCTGCGAAGCTATGGATTATTAGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCaa                                                                              <  2:257791/140‑1 (MQ=255)
                                                        gggTTGTTTCGTACCGTGTGGGGTGCGAATCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGa                                                                        <  1:418210/140‑1 (MQ=255)
                                                         ggTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGGGGAGCTCAACGGCGATTACCGCTCGTCCAACGGTGaa                                                                       >  1:257791/1‑140 (MQ=255)
                                                             gttTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATg                                                                   >  1:629614/1‑140 (MQ=255)
                                                             gttTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATg                                                                   <  2:596792/140‑1 (MQ=255)
                                                             gttTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAACATg                                                                   <  2:717858/140‑1 (MQ=255)
                                                             gttTCGTACCTTGTGCGCTGGGAAGCTATTGATTATTTGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATg                                                                   <  1:254612/140‑1 (MQ=255)
                                                              tttCGTACCTTGTGGGGTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGc                                                                  >  1:64711/1‑140 (MQ=255)
                                                              tttCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGc                                                                  >  1:720520/1‑140 (MQ=255)
                                                                   tACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTgggg                                                             >  1:650054/1‑140 (MQ=255)
                                                                   tACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCCACGGCGATTACCGATCCTCCAGCGGTGAATATGCTgggg                                                             >  2:305970/1‑140 (MQ=255)
                                                                   tACCTTGTGGGCTGGGAAGCTATGGATTATTGACATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACg                                                                            >  2:520339/1‑125 (MQ=255)
                                                                   tACCTTGTGGGCTGGGAAGCTATGGATTATTGACATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACg                                                                            <  1:520339/125‑1 (MQ=255)
                                                                    aCCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATGCGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAAc                                                                             <  1:217666/123‑1 (MQ=255)
                                                                    aCCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAAc                                                                             >  2:217666/1‑123 (MQ=255)
                                                                        tgtgGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTggg                                                              >  2:583003/1‑134 (MQ=255)
                                                                        tgtgGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTggg                                                              <  1:583003/134‑1 (MQ=255)
                                                                               tGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCg                                                 >  1:136872/1‑140 (MQ=255)
                                                                                      cTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCgcagcagc                                          >  1:209844/1‑140 (MQ=255)
                                                                                                  ggCATGGGGCTGGCGTTTGCGTGAATCGCGCTGGCTTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTa                              <  2:355464/140‑1 (MQ=255)
                                                                                                      tGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCt                          >  1:809278/1‑140 (MQ=255)
                                                                                                            tGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGt                    <  2:720520/140‑1 (MQ=255)
                                                                                                            tGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCCGGATGt                    <  2:209844/140‑1 (MQ=255)
                                                                                                            tGGCGTTTGCGTTAAGCGCGCTTGCTTATCCGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGt                    <  1:305970/140‑1 (MQ=255)
                                                                                                             ggCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTc                   >  2:122040/1‑140 (MQ=255)
                                                                                                             ggCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTc                   <  2:136872/140‑1 (MQ=255)
                                                                                                              gCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTcc                  <  2:13369/140‑1 (MQ=255)
                                                                                                                   tGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCa             >  1:627643/1‑140 (MQ=255)
                                                                                                                   tGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCa             >  1:25425/1‑140 (MQ=255)
                                                                                                                   tGCGTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCa             >  1:201528/1‑140 (MQ=255)
                                                                                                                      gTTAATCGCGCTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAAc          <  2:629614/140‑1 (MQ=255)
                                                                                                                              cgcTGGCGTATACGTTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTc  >  1:42616/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                              
CACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTC  >  NC_000913/3969035‑3969299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: