Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 594,043 C→T 54.1% intergenic (+131/‑25) ftsK → / → lolA DNA‑binding membrane protein required for chromosome resolution and partitioning/chaperone for lipoproteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE594,0430CT54.1% 0.0 / 28.4 37intergenic (+131/‑25)ftsK/lolADNA‑binding membrane protein required for chromosome resolution and partitioning/chaperone for lipoproteins
Reads supporting (aligned to +/- strand):  ref base C (15/2);  new base T (7/13);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 2.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTAC  >  minE/593983‑594101
                                                            |                                                          
ttttCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGa                                                    <  1:178643/69‑1 (MQ=255)
  ttcttcTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGAc                                                   >  1:2636493/1‑68 (MQ=255)
  ttcttcTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGAc                                                   >  1:1754437/1‑68 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGaataat                                  >  1:508285/1‑68 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:3077595/1‑69 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:465866/1‑69 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:1963053/1‑69 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:1955656/1‑69 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:195000/1‑69 (MQ=255)
                   tGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATTGGGAGTGACGTAATTTGAGGAATAatg                                 >  1:159524/1‑69 (MQ=255)
                                   aTAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2533854/68‑1 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:3181455/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:269940/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:3107655/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:3062062/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:2581561/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:2528029/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:1697624/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:1661439/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:1135693/1‑69 (MQ=255)
                                             gagGGAACTCCCGATCGGGAGTGACGTAAGTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGc       >  1:2410834/1‑69 (MQ=255)
                                                  aaCTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTAc  >  1:2341231/1‑69 (MQ=255)
                                                  aaCTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTAc  >  1:146547/1‑69 (MQ=255)
                                                  aaCTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTAc  >  1:825381/1‑69 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGTAATAATGATGAAAAAAATTGcc                  <  1:2290217/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:384296/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:938569/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:3054378/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2946273/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2857865/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2731407/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2581893/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2260633/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:2092957/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:1967170/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:1325758/46‑1 (MQ=255)
                                                         gATTGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGcc                  <  1:1292978/46‑1 (MQ=255)
                                                            |                                                          
TTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTAC  >  minE/593983‑594101

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