Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 90,917 T→C 57.0% T291T (ACT→ACC ftsZ → GTP‑binding tubulin‑like cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE90,9170TC57.0% 19.0 / 46.0 65T291T (ACT→ACCftsZGTP‑binding tubulin‑like cell division protein
Reads supporting (aligned to +/- strand):  ref base T (12/16);  new base C (13/24);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TCGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGC  >  minE/90868‑90979
                                                 |                                                              
tCGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTtctc                                             <  1:86790/69‑1 (MQ=255)
tCGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTtctc                                             <  1:3274632/69‑1 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTAt                                                        >  1:340770/1‑57 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:2669205/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:776773/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:71182/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:444187/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:3386032/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:3348100/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:3200143/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:2413364/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:2185895/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:2114612/1‑58 (MQ=255)
 cGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:1628342/1‑58 (MQ=255)
 cGAAACGATAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACCGTGGTTATc                                                       >  1:2621423/1‑58 (MQ=255)
    aaCGGTAGGTAACACCAACCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGGACTTCTCTTg                                          <  1:1396625/68‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:1403157/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:1286550/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:724297/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:545552/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:1479979/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:2414168/55‑1 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:2550597/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:3242637/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:358581/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:423879/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:2691077/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:2579533/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:1545829/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:776286/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:838735/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:192348/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:84618/1‑69 (MQ=255)
                 acCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATatgaat                            >  1:1808354/1‑69 (MQ=255)
                               tGCTTCCGACAACTCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2101981/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1076960/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1152365/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:669667/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:780519/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:3397358/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:3394004/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1326012/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2639197/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1835230/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1873676/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:1966901/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2084267/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2112626/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2179130/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:227299/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2640403/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2880922/69‑1 (MQ=255)
                               tGCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:2918848/69‑1 (MQ=255)
                                gCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:525616/68‑1 (MQ=255)
                                gCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:673415/68‑1 (MQ=255)
                                gCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:113035/68‑1 (MQ=255)
                                gCTTCCGACAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCAGATATGAATGACGAGCTGCGCGt              <  1:1041962/68‑1 (MQ=255)
                                    ccGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTa             <  1:1422463/65‑1 (MQ=255)
                                    ccGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTa             <  1:789798/65‑1 (MQ=255)
                                    ccGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTCGCGTa             <  1:3428690/64‑1 (MQ=255)
                                     cGAAAACGCGACTGTGGTTATCGGTACTTCTCTTg                                          <  1:718352/35‑1 (MQ=255)
                                      gACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCgt        <  1:1320614/68‑1 (MQ=255)
                                       aCAACGCGACCGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGt              <  1:647780/61‑1 (MQ=255)
                                           cgcgACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGc  <  1:2951441/69‑1 (MQ=255)
                                            gcgACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGc  <  1:2968402/68‑1 (MQ=255)
                                                 |                                                              
TCGAAACGGTAGGTAACACCATCCGTGCATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGC  >  minE/90868‑90979

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: