Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 591,494 A→G 23.1% P524P (CCA→CCG ftsK → DNA‑binding membrane protein required for chromosome resolution and partitioning

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE591,4940AG23.1% 109.0 / 27.6 86P524P (CCA→CCGftsKDNA‑binding membrane protein required for chromosome resolution and partitioning
Reads supporting (aligned to +/- strand):  ref base A (27/39);  new base G (10/10);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGAACGTGAACAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGT  >  minE/591437‑591538
                                                         |                                            
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:992537/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:979117/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:93249/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:3388159/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:3333613/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:3205927/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:3109781/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:2921323/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:2821666/1‑69 (MQ=255)
tgaacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCaa                                   >  1:1331677/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:1008427/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:878847/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:3160256/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:1891916/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:1799185/1‑69 (MQ=255)
  aacgtgaacAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                 >  1:1274909/1‑69 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:1636853/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:993666/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:840089/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:581053/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:1143927/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:3251482/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:3191618/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:3150189/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:2676211/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:2437615/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:1437090/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:1785000/1‑68 (MQ=255)
        aacaacTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                            >  1:1484125/1‑68 (MQ=255)
              ttGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAg                     <  1:1140757/69‑1 (MQ=255)
              ttGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAg                     <  1:92276/69‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAagaaccagaac                                          >  1:253529/1‑36 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGAAAGCAGCTTCTGtt          <  1:1924750/68‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGtt          <  1:3140996/68‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCg                                        >  1:1120354/1‑38 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           >  1:1626873/1‑67 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:1902401/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:912975/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:530406/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           >  1:1758481/1‑67 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           >  1:1763028/1‑67 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:1278271/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:2061715/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:2062451/67‑1 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           >  1:2608435/1‑67 (MQ=255)
                          aTCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:2440580/67‑1 (MQ=255)
                             aaCCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgc        >  1:2599332/1‑67 (MQ=255)
                             aaCCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgc        >  1:2952874/1‑67 (MQ=255)
                             aaCCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgc        >  1:2994567/1‑67 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:2601153/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:513010/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:1786686/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:1229149/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:554182/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:681733/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:1890828/48‑1 (MQ=255)
                                cGATTCCAGAACCGGTTAAAGAACCGGAACCGATCAAATCTTCGCTGa                        <  1:1987217/48‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:527110/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:1005421/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2762533/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:3169108/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2831100/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2834891/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:1217871/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCTAATCTTCGCTGAAAGCACATTCTGTTGCAGCAGt  <  1:2132138/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGTACCTTCTGTTGCAGCAGt  <  1:3373899/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGTACCTTCTGTTGCAGCAGt  <  1:1707451/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:3102832/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:1143801/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:854648/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2273379/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2934683/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:675821/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:618675/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:3101848/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:541268/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:515680/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2823543/69‑1 (MQ=255)
                                 gATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:121076/69‑1 (MQ=255)
                                  aTTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGTACCTTCTGTTGCAGCAGt  <  1:2442144/68‑1 (MQ=255)
                                  aTTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt  <  1:2601626/68‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:3223745/45‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:1681830/45‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:644275/45‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:374778/45‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:1186634/45‑1 (MQ=255)
                                                tAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt           <  1:3294034/45‑1 (MQ=255)
                                                      accataaccGATCTAATCTTCGCTGAAAGCACATTCTGTTGCAGCAGt  <  1:27510/48‑1 (MQ=255)
                                                         |                                            
TGAACGTGAACAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGT  >  minE/591437‑591538

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: