Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 593,687 T→C 41.8% G1255G (GGT→GGC 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
*minE593,6870TC41.8% 27.6 / 63.1 86G1255G (GGT→GGCftsKDNA‑binding membrane protein required for chromosome resolution and partitioning
Reads supporting (aligned to +/- strand):  ref base T (36/14);  new base C (18/18);  total (54/32)
Fisher's exact test for biased strand distribution p-value = 4.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTC  >  minE/593621‑593750
                                                                  |                                                               
cTGGAAATCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:2139329/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGgtggtg                                                                >  1:1367604/1‑68 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:2900940/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:2697625/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:482203/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:765233/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:2494718/1‑69 (MQ=255)
cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                                               >  1:1395902/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGATGGTGGTGCggg                                                            >  1:3092876/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:2490657/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:3334696/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:1753946/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:1595759/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:1457678/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:2393318/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:3009888/1‑69 (MQ=255)
   gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                                            >  1:739999/1‑69 (MQ=255)
     aaGCGCGTGGTCGCCCACTGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:3416103/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1118903/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:360298/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1068362/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:467809/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:2262812/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:3145067/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:16639/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:728158/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:2604721/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1487927/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1473721/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1367202/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                                          >  1:1406333/1‑67 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:2640226/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:2753364/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:2539448/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:2812044/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:400548/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1966494/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:971188/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:903404/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1263174/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1841722/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:780889/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1273479/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1389078/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:140743/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:661280/67‑1 (MQ=255)
     aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:1523884/67‑1 (MQ=255)
     aaGCGCGTGGGCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                                            <  1:224141/67‑1 (MQ=255)
             ggTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGAt                                                  <  1:1865064/69‑1 (MQ=255)
              gTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGAt                                                  >  1:2245536/1‑68 (MQ=255)
              gTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATg                                                 >  1:2696398/1‑69 (MQ=255)
              gTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATg                                                 >  1:2018888/1‑69 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:1704452/69‑1 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:1784796/69‑1 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:2203765/69‑1 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:2039137/69‑1 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:2782489/69‑1 (MQ=255)
                       aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                                        <  1:1913924/69‑1 (MQ=255)
                        gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                                       <  1:633093/69‑1 (MQ=255)
                        gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                                       <  1:2420283/69‑1 (MQ=255)
                        gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                                       <  1:1184226/69‑1 (MQ=255)
                        gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                                       <  1:1112506/69‑1 (MQ=255)
                             ttGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:3270948/1‑69 (MQ=255)
                             ttGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:980917/1‑69 (MQ=255)
                             ttGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTg                                   >  1:1048500/1‑68 (MQ=255)
                                  ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATcc                              >  1:1418665/1‑68 (MQ=255)
                                  ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg                             >  1:826381/1‑69 (MQ=255)
                                  ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg                             >  1:3082089/1‑69 (MQ=255)
                                  ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg                             >  1:787323/1‑69 (MQ=255)
                                  ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg                             >  1:291914/1‑69 (MQ=255)
                                        aCCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTc                       >  1:2171024/1‑69 (MQ=255)
                                        aCCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTc                       >  1:1999302/1‑69 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1804752/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:110786/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  <  1:1282902/46‑1 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:137620/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1395200/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1509894/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1816107/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1928556/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:3344883/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:2304648/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:2305110/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:2376527/1‑46 (MQ=255)
                                                    gaaagcgaaGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTgg                                  >  1:1021099/1‑46 (MQ=255)
                                                            aGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGt   <  1:520262/69‑1 (MQ=255)
                                                            aGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGt   <  1:1268282/69‑1 (MQ=255)
                                                             ggtggtGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGt   <  1:967047/68‑1 (MQ=255)
                                                             ggtggtGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTc  <  1:2434851/69‑1 (MQ=255)
                                                              gtggtgCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGt   <  1:1414537/67‑1 (MQ=255)
                                                                  |                                                               
CTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTC  >  minE/593621‑593750

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: