Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 636,849 C→G 38.0% L1275L (CTC→CTG mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE636,8490CG38.0% 34.2 / 42.8 50L1275L (CTC→CTGmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base C (25/6);  new base G (17/2);  total (42/8)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAC  >  minE/636798‑636901
                                                   |                                                    
cAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:65708/69‑1 (MQ=255)
cAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:178423/69‑1 (MQ=255)
      gtggTTGCAGAACGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAacgca                               <  1:2324041/67‑1 (MQ=255)
      ggggTTGCAGAACGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAacgca                               <  1:209567/69‑1 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGcc                           >  1:2376107/1‑68 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGcc                           >  1:714900/1‑68 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1000140/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:2602740/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:323119/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1995542/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1828048/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1706891/1‑69 (MQ=255)
           tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1604302/1‑69 (MQ=255)
            gCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCAt                         <  1:3433839/69‑1 (MQ=255)
            gCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCAt                         <  1:3373915/69‑1 (MQ=255)
             cAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:491180/56‑1 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGTCATGCtg                     >  1:665822/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3327198/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:30678/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:444230/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3344189/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3386944/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:374859/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:28173/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:921127/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:971496/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2840384/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2753598/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:1429371/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:152387/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2536804/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2280454/1‑68 (MQ=255)
                aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCt                      >  1:2598300/1‑68 (MQ=255)
                   gTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTg                  <  1:2710583/69‑1 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACTTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2034366/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2912774/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:771969/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:1731215/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:507706/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2178640/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:400830/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2532771/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:3427327/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2701611/1‑69 (MQ=255)
                      tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:3287140/1‑69 (MQ=255)
                                    gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:2566288/1‑68 (MQ=255)
                                    gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:1696169/1‑68 (MQ=255)
                                    gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:325110/1‑68 (MQ=255)
                                    gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:3101440/1‑68 (MQ=255)
                                    gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:1013241/1‑68 (MQ=255)
                                                   |                                                    
CAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAC  >  minE/636798‑636901

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