Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,461,741 T→C 100% S327S (TCT→TCC gidA → glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,461,7410TC78.3% 89.1 / 16.1 64S327S (TCT→TCCgidAglucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base T (3/11);  new base C (50/0);  total (53/11)
Fisher's exact test for biased strand distribution p-value = 4.90e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGA  >  minE/2461697‑2461806
                                            |                                                                 
tcCACCGGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:4310/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:275725/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:2611104/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:2586610/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:2462982/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:314098/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:3169241/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:1859496/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:399715/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCCATGCAGGGGATGGAAAAcg                                                >  1:547827/1‑64 (MQ=255)
tcCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCAGCTCCATGCAGGGGATGGAAAACGCGaa                                            >  1:1513834/1‑68 (MQ=255)
               ttCGATGTGCATATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:3335919/69‑1 (MQ=255)
               ttCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:1921697/69‑1 (MQ=255)
               ttCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:3243418/69‑1 (MQ=255)
               ttCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:1462978/69‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:3306251/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:1822708/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:2659805/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:1149149/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:1311780/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:3210801/62‑1 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCa                     >  1:2962138/1‑69 (MQ=255)
                      tgCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCa                     >  1:1153640/1‑69 (MQ=255)
                             gCAAATCGTCCGCTCTATGCAGGGGATGTAAAACGCGAAGATCGTGCGTCCGGGt                            <  1:843553/55‑1 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:2782671/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:3098757/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:2711509/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:1072484/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:2590895/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:2350577/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:818882/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:1908925/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:1343304/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGt               >  1:1812106/1‑68 (MQ=255)
                             gCAAATCGTCCGCTCCATGAAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGACATTGAGt               >  1:926996/1‑68 (MQ=255)
                                   cGTCCGCTCCATGGAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:421635/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATg            >  1:2608362/1‑65 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAct         >  1:794865/1‑68 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:737243/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:702522/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:353135/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:614871/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:3221797/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:53898/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:561742/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:375526/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2902990/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1031666/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:890141/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1369618/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1484118/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1580020/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1582363/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:167301/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1755788/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1791402/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:1971344/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2044769/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2077749/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2148001/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2239317/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2724670/1‑67 (MQ=255)
                                   cGTCCGCTCCATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGAc          >  1:2786325/1‑67 (MQ=255)
                                       cGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACttcttc    >  1:477927/1‑69 (MQ=255)
                                         ctctTTGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3122300/69‑1 (MQ=255)
                                         ctctATGGAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:2927695/69‑1 (MQ=255)
                                         ctctATGCAGTGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3284576/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACTCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3392683/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:767604/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:990106/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3418089/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3325437/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:3177879/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:300887/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:2932381/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:2021157/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:1971311/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:1305459/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:1190590/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:1042356/69‑1 (MQ=255)
                                         ctctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:1035999/69‑1 (MQ=255)
                                          tctATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGa  <  1:2356508/68‑1 (MQ=255)
                                            |                                                                 
TCCACCAGCCTGCCGTTCGATGTGCAGATGCAAATCGTCCGCTCTATGCAGGGGATGGAAAACGCGAAGATCGTGCGTCCGGGTTATGCCATTGAGTATGACTTCTTCGA  >  minE/2461697‑2461806

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: