Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 593,672 C→T 50.0% S1250S (AGC→AGT 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,6720CT50.0% ‑3.0 / 64.7 72S1250S (AGC→AGTftsKDNA‑binding membrane protein required for chromosome resolution and partitioning
Reads supporting (aligned to +/- strand):  ref base C (25/11);  new base T (18/18);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTC  >  minE/593611‑593729
                                                             |                                                         
tggtgCATGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2397344/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgagagc                                                     <  1:2818392/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:675390/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2398910/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2337514/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2316038/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2181628/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:853155/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2165816/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:209921/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:2095735/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:85335/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:3143666/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:171629/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:3220732/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:3307711/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:1472463/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:1093012/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:854800/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:407230/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:523444/68‑1 (MQ=255)
tggtgCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagcgaaagc                                                     <  1:5603/68‑1 (MQ=255)
tggtgCACGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCGGACagtgaaagc                                                     <  1:82012/68‑1 (MQ=255)
         aCTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACagtgaaagc                                                     <  1:3159099/59‑1 (MQ=255)
          cTGGAAATCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:2139329/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGgtggtg                                           >  1:1367604/1‑68 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:2494718/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:765233/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:2900940/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:1395902/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:482203/1‑69 (MQ=255)
          cTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGc                                          >  1:2697625/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGATGGTGGTGCggg                                       >  1:3092876/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:2490657/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:2393318/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:1457678/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:3009888/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:1753946/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:3334696/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:1595759/1‑69 (MQ=255)
             gAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCggg                                       >  1:739999/1‑69 (MQ=255)
               aaGCGCGTGGTCGCCCACTGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:3416103/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:467809/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:360298/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:2604721/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:3145067/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:728158/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1068362/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1367202/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1406333/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1473721/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1487927/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:1118903/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:16639/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTgggcg                                     >  1:2262812/1‑67 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:903404/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:1273479/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:400548/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:1389078/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:140743/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:1523884/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:1966494/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:2812044/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:2753364/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:971188/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:2640226/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:661280/67‑1 (MQ=255)
               aaGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTggg                                       <  1:2539448/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:2018888/1‑69 (MQ=255)
                        gTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATg                            >  1:2696398/1‑69 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:1784796/69‑1 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:1704452/69‑1 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:1913924/69‑1 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:2782489/69‑1 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:2039137/69‑1 (MQ=255)
                                 aGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaag                   <  1:2203765/69‑1 (MQ=255)
                                  gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                  <  1:2420283/69‑1 (MQ=255)
                                  gTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTgaaga                  <  1:633093/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:3082089/1‑69 (MQ=255)
                                            ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg        >  1:291914/1‑69 (MQ=255)
                                            ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg        >  1:787323/1‑69 (MQ=255)
                                            ggCATCACCTCCGACAGTGAAAGCGAAGGTGGCGTGGGCGGTTTCGATGGCGCTGAAGAACTGGATCCg        >  1:826381/1‑69 (MQ=255)
                                                  aCCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTc  >  1:2171024/1‑69 (MQ=255)
                                                  aCCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTc  >  1:1999302/1‑69 (MQ=255)
                                                             |                                                         
TGGTGCAGGACTGGAAAGCGCGTGGTCGCCCACAGTATGTTGATGGCATCACCTCCGACAGCGAAAGCGAAGGTGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTC  >  minE/593611‑593729

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: