Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,489,400 A→G 53.7% A133A (GCA→GCG dnaN → DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,489,4000AG53.7% 5.3 / 43.0 54A133A (GCA→GCGdnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (14/11);  new base G (9/20);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 9.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGCGGCGGATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCTATTACTTA  >  minE/2489334‑2489466
                                                                  |                                                                  
gCCAGCGGCGGATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAAc                                                                  <  1:3195787/69‑1 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAAc                                                                  >  1:3138873/1‑63 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:1084966/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:2581062/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:30521/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:511931/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:2085746/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:3154334/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:3410053/1‑64 (MQ=255)
      ggcggATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACg                                                                 >  1:1989944/1‑64 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCTAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1867562/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACTCTGCCGCAGGCGACg                                                                 <  1:2886170/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1111838/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:3183954/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:354115/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:357813/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:478584/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:593913/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:3028224/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2956978/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:626178/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2547695/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1732976/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1706029/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1767769/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2033373/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2428376/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:1383751/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2314188/55‑1 (MQ=255)
               cccGAACCTCGATGACTGGCAGAGCGAAGTCGAATTTACCCTGCCGCAGGCGACg                                                                 <  1:2101370/55‑1 (MQ=255)
                               tGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGttt                                   >  1:434489/1‑69 (MQ=255)
                               tGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGttt                                   >  1:998194/1‑69 (MQ=255)
                                      gTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATgg                            >  1:800683/1‑69 (MQ=255)
                                      gTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATgg                            >  1:3243445/1‑69 (MQ=255)
                                              gAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTg                                                     >  1:200591/1‑36 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:2495846/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:1320146/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:358806/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:1411766/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:98546/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:3080734/69‑1 (MQ=255)
                                                  ttACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACg                <  1:2855350/69‑1 (MQ=255)
                                                       cTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGc           <  1:2880675/69‑1 (MQ=255)
                                                       cTGCCGCAGGCAACGATGAAGCGTCTGAGTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGc           <  1:2975771/69‑1 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:3350707/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:1770468/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:369288/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:493261/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:1364224/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:2436019/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:1253881/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:766978/1‑69 (MQ=255)
                                                        tGCCGCAGGCGACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCt          >  1:250340/1‑69 (MQ=255)
                                                                gCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGTGCATCAGGACGTTCGCTATTACTTa  <  1:3252754/69‑1 (MQ=255)
                                                                  |                                                                  
GCCAGCGGCGGATTTCCCGAACCTCGATGACTGGCAGAGTGAAGTCGAATTTACCCTGCCGCAGGCAACGATGAAGCGTCTGATTGAAGCGACCCAGTTTTCTATGGCGCATCAGGACGTTCGCTATTACTTA  >  minE/2489334‑2489466

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: