Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,429,321 T→C 39.2% M88V (ATG→GTG)  wzzE ← Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,429,3210TC39.2% 56.5 / 90.2 79M88V (ATG→GTG) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base C (19/12);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

GCCAGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCATATTCACCGTT  >  minE/2429257‑2429391
                                                                |                                                                      
gCCAGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTg                                                                   <  1:547316/70‑1 (MQ=255)
  cAGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAAc                                                                <  1:2622734/71‑1 (MQ=255)
  cAGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAAc                                                                <  1:2447942/71‑1 (MQ=255)
   aGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACg                                                               <  1:924303/71‑1 (MQ=255)
   aGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACg                                                               <  1:185071/71‑1 (MQ=255)
     ttGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAAc                                                                <  1:2895678/68‑1 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACgg                                                              >  1:2504740/1‑68 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACgg                                                              >  1:3315960/1‑68 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACgg                                                              >  1:3500917/1‑68 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACgg                                                              >  1:3525196/1‑68 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACgg                                                              >  1:540473/1‑68 (MQ=255)
       gCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACa                                                           >  1:1263212/1‑71 (MQ=255)
         aTAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATc                                                         <  1:714988/71‑1 (MQ=255)
          tAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATcc                                                        <  1:3246396/71‑1 (MQ=255)
          tAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATc                                                         >  1:31935/1‑70 (MQ=255)
          tAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATcc                                                        <  1:3119291/71‑1 (MQ=255)
          tAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATcc                                                        <  1:2455715/71‑1 (MQ=255)
            aCAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACg                                                               >  1:1346540/1‑62 (MQ=255)
             cAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCAt                                                                        >  1:1156060/1‑52 (MQ=255)
             cAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAgg                                                     <  1:1849367/71‑1 (MQ=255)
              aaaCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGt                                                    >  1:1696006/1‑71 (MQ=255)
              aaaCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACa                                                           >  1:3593754/1‑64 (MQ=255)
               aaCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGt                                                    >  1:3556421/1‑70 (MQ=255)
               aaCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGt                                                    >  1:644498/1‑70 (MQ=255)
                aCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTa                                                  <  1:3196334/71‑1 (MQ=255)
                aCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTa                                                  >  1:1944787/1‑71 (MQ=255)
                aCTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTa                                                  >  1:1759875/1‑71 (MQ=255)
                 ctctTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTAc                                                 <  1:3303564/71‑1 (MQ=255)
                 ctctTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTAc                                                 >  1:834540/1‑71 (MQ=255)
                    tttATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGc                                               >  1:2010807/1‑70 (MQ=255)
                     ttATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCa                                              >  1:1551769/1‑70 (MQ=255)
                      tataGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             >  1:3559689/1‑70 (MQ=255)
                       ataGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGtt                                                                     >  1:1630024/1‑45 (MQ=255)
                       ataGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGtt                                                                     >  1:1584301/1‑45 (MQ=255)
                       ataGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCaaa                                            >  1:704229/1‑70 (MQ=255)
                        taGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaaaaa                                          >  1:1812497/1‑71 (MQ=255)
                        taGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCaaaaa                                          >  1:708006/1‑71 (MQ=255)
                        taGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCaaaa                                           <  1:2379172/70‑1 (MQ=255)
                          ggCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAAt                                         <  1:3306147/70‑1 (MQ=255)
                           gCTTCGTCCATGACCATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgc                                      <  1:3628642/71‑1 (MQ=255)
                            cTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgc                                      <  1:2915930/71‑1 (MQ=255)
                            cTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgc                                      <  1:2384419/71‑1 (MQ=255)
                             ttCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgct                                     >  1:3555059/1‑71 (MQ=255)
                                ggccATGACCGATGGGTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATtgctgct                                  <  1:3533881/69‑1 (MQ=255)
                                gTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgctgct                                  <  1:644628/71‑1 (MQ=255)
                                gTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgctgct                                  <  1:1325029/71‑1 (MQ=255)
                                  ccATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTa                                                  >  1:2673022/1‑53 (MQ=255)
                                  ccATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTa                                                  >  1:2687893/1‑53 (MQ=255)
                                   cATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATtgctgctgc                                >  1:3116106/1‑70 (MQ=255)
                                   cATGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCg                               <  1:2207106/71‑1 (MQ=255)
                                    aTGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGc                                               >  1:2757347/1‑54 (MQ=255)
                                    aTGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGc                                               >  1:3713097/1‑54 (MQ=255)
                                    aTGACCGATGGTTGGTCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGc                                               >  1:2131458/1‑54 (MQ=255)
                                     tGACCGATGGTTGATCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             >  1:2770613/1‑55 (MQ=255)
                                      gACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAg                             >  1:1499991/1‑70 (MQ=255)
                                      gACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAg                             >  1:1880210/1‑70 (MQ=255)
                                      gACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAg                             >  1:1747700/1‑70 (MQ=255)
                                       aCCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGt                            <  1:567342/70‑1 (MQ=255)
                                       aCCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGt                            <  1:3116165/70‑1 (MQ=255)
                                       aCCGATGGTTGGACGGCAGAAGCCATGTTTGAACgg                                                              <  1:583762/36‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:354504/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:2365807/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:339099/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:3649841/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:2757902/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCaa                                             <  1:2105912/42‑1 (MQ=255)
                                                  gtCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTa                           <  1:562441/60‑1 (MQ=255)
                                                  gtCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATccccc                  >  1:3438588/1‑69 (MQ=255)
                                                  gtCGGCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATccccc                  >  1:2553210/1‑69 (MQ=255)
                                                     ggCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCa             >  1:3049481/1‑71 (MQ=255)
                                                      gCAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCa             >  1:14983/1‑70 (MQ=255)
                                                       cAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCata           <  1:1961737/71‑1 (MQ=255)
                                                       cAGAAGCCACGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCata           >  1:972938/1‑71 (MQ=255)
                                                         gAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGa                              <  1:3232260/50‑1 (MQ=255)
                                                          aaGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCAGCGAGTAATATCCCCCCAGCATATTc        <  1:647217/71‑1 (MQ=255)
                                                            gCCACGTTTGAACGGACATCCAGGTTACGCaaaaa                                          >  1:2122894/1‑35 (MQ=255)
                                                            gCCACGTTTGAACGGACATCCAGGTTACGCaaaaa                                          >  1:2377163/1‑35 (MQ=255)
                                                            gCCACGTTTGAACGGACATCCAGGTTACGCaaaaa                                          >  1:3249031/1‑35 (MQ=255)
                                                               aTGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCATATTCACCGt   >  1:2420187/1‑71 (MQ=255)
                                                               aTGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCATATTCACCGt   >  1:2457498/1‑71 (MQ=255)
                                                                tGTTTGAAGGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCATATTCACCGtt  >  1:3131841/1‑71 (MQ=255)
                                                                |                                                                      
GCCAGTTGCATAACAAACTCTTTATAGGCTTCGTCCATGACCGATGGTTGGTCGGCAGAAGCCATGTTTGAACGGACATCCAGGTTACGCAAAAATTGCTGCTGCGAGTAATATCCCCCCAGCATATTCACCGTT  >  minE/2429257‑2429391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: