Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,340 C→A 7.7% W72L (TGG→TTG)  wcaJ ← colanic biosynthesis UDP‑glucose lipid carrier transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,121,3400CA 7.7% 174.2 / 3.0 52W72L (TGG→TTG) wcaJcolanic biosynthesis UDP‑glucose lipid carrier transferase
Reads supporting (aligned to +/- strand):  ref base C (20/28);  new base A (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CGCATACCACGCCAGCCAGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGT  >  NC_000913/2121204‑2121476
                                                                                                                                        |                                                                                                                                        
cgcATACCACGCCAGCCAGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATg                                                                                                                                       >  1:222559/1‑140 (MQ=255)
     aCCACGCCAGCCAGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCACATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGa                                                                                                                                  >  2:371234/1‑140 (MQ=255)
          gccagccaGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCTGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCACATTCTGTCGCTGCCCGAACACCGCGCCAGGAGCGATaaaa                                                                                                                             >  2:1454175/1‑140 (MQ=255)
          gccagccaGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGagacaa                                                                                                                             >  2:148529/1‑135 (MQ=255)
          gccagccaGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATaaaa                                                                                                                             >  1:632860/1‑140 (MQ=255)
             agccagATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATc                                                                                                                          >  1:1153509/1‑140 (MQ=255)
             agccagATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATc                                                                                                                          >  1:1153529/1‑140 (MQ=255)
             agccagATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATc                                                                                                                          >  1:1153656/1‑140 (MQ=255)
             agccagATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCAATGAGCGATAAAAATc                                                                                                                          >  1:1540950/1‑140 (MQ=255)
             agccagATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCAATGAGCGATAAAAATc                                                                                                                          >  1:1541099/1‑140 (MQ=255)
                    ttttCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATg                                                                                                                   <  2:544966/140‑1 (MQ=255)
                         aGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATgccgcc                                                                                                              <  1:1159188/140‑1 (MQ=255)
                             gCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTCGCAACAGGGCAAATTCTGTCGCTGCACGAACACCGCGCCCTGAGCGATAAAAATAGGTGATGCCGCCCAGc                                                                                                          >  2:1378156/1‑140 (MQ=255)
                             gCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTATGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGc                                                                                                          >  1:1422723/1‑140 (MQ=255)
                                   cGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTgg                                                                                                    <  2:1200306/140‑1 (MQ=255)
                                           ttgttgAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAAcaccac                                                                                            <  1:1554990/140‑1 (MQ=255)
                                                gAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCg                                                                                       <  1:262337/140‑1 (MQ=255)
                                                gAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCg                                                                                       <  1:264433/140‑1 (MQ=255)
                                                      caccaGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCa                                                                                 >  2:1209837/1‑140 (MQ=255)
                                                       accaGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAg                                                                                >  2:985314/1‑140 (MQ=255)
                                                       accaGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAg                                                                                <  1:586379/140‑1 (MQ=255)
                                                       accaGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAg                                                                                <  1:1666122/140‑1 (MQ=255)
                                                        ccaGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAgc                                                                               >  1:1033874/1‑140 (MQ=255)
                                                              cGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGccacc                                                                         <  1:1146567/140‑1 (MQ=255)
                                                                  gcTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACa                                                                     <  2:230673/140‑1 (MQ=255)
                                                                  gcTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACa                                                                     <  2:1466424/140‑1 (MQ=255)
                                                                  gcTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACa                                                                     <  1:1646710/140‑1 (MQ=255)
                                                                     gAAAATCACGCTTAAGGTCCAGTTATGTATCAACAGGTCAAATTCTGTCGCTGCCCGACCACCGCGGCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGGTTAATCAGCGCCTCCAACAGGt                                                                  <  2:1036313/140‑1 (MQ=255)
                                                                      aaacTCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGAAATCAGCGCCACCAACAGGTg                                                                 <  2:1502983/140‑1 (MQ=255)
                                                                      aaaaTCACGCTTAAGGTCCAGTTTTGTCGCACCAGGGCAAATTCTGTCGCTGCCCGATCACCGCGCCATGAGCGATAATATTCGGTGTTGCCGCCCAGCTTCTGGAACACCACCAGCGTAATCAGCGCCTCCAACGGGTg                                                                 >  2:687143/1‑140 (MQ=255)
                                                                       aaaTCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGc                                                                <  1:1636890/140‑1 (MQ=255)
                                                                          tCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATg                                                             <  2:1297972/140‑1 (MQ=255)
                                                                           cACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCTCTGCCCGTACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGt                                                            <  2:858904/140‑1 (MQ=255)
                                                                           cACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGt                                                            <  2:812740/140‑1 (MQ=255)
                                                                             cGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTag                                                          <  1:654905/140‑1 (MQ=255)
                                                                                ttAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAgg                                                       >  1:1205542/1‑140 (MQ=255)
                                                                                    gctCCAGTCTTGGCGCCCCACGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAACAATCGGTGATGCCGCCCAGCATCTGGACCACCACCAGCGTAATCAGCGCCACCAACAGGTTCATGTAGAGGAATg                                                   <  2:352406/138‑1 (MQ=255)
                                                                                      tCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGAc                                                 <  1:1541262/140‑1 (MQ=255)
                                                                                      tCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGAc                                                 <  1:1541442/140‑1 (MQ=255)
                                                                                                tAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGa                                       >  1:704077/1‑140 (MQ=255)
                                                                                                          gCATATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAAc                             <  1:297762/140‑1 (MQ=255)
                                                                                                          gCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAAc                             <  1:297081/140‑1 (MQ=255)
                                                                                                          gCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAAc                             <  1:440758/140‑1 (MQ=255)
                                                                                                                      gCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTcc                 <  1:442683/140‑1 (MQ=255)
                                                                                                                         gCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGc              >  1:437859/1‑140 (MQ=255)
                                                                                                                            cGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCaaa           <  2:222559/140‑1 (MQ=255)
                                                                                                                              aaCACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACa         <  2:428046/140‑1 (MQ=255)
                                                                                                                                cacCGCGCAATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACctg       <  2:1540950/140‑1 (MQ=255)
                                                                                                                                cacCGCGCAATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACatg       <  2:1541099/140‑1 (MQ=255)
                                                                                                                                    gcgcCATGAGCGTGATAAATCGGTGATGCCGCCCAGCATCTGGACCACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACGTCGCTAACCAGCCATAGTCCGGCTAACTTGAAgg   >  2:1666611/1‑140 (MQ=255)
                                                                                                                                    gcgcCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATgg   >  1:339242/1‑140 (MQ=255)
                                                                                                                                     cgcCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGt  >  2:864757/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CGCATACCACGCCAGCCAGATTTTCAGTTGCGTGTCGAAATCATTGTTGAACGCCACCAGTCCGGCGCTGAAAATCACGCTTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGT  >  NC_000913/2121204‑2121476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑