Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,517,725 G→A 64.3% G450G (GGC→GGT spoT ← bifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,517,7250GA64.3% 22.9 / 22.3 42G450G (GGC→GGTspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base A (18/9);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GCGCGCGCTTTCGAGCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTA  >  minE/2517672‑2517784
                                                     |                                                           
gcgcgcgcTTTCGAGCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAAt                                               <  1:147905/68‑1 (MQ=255)
 cgcgcgcTTTCGAGCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTGATg                                              <  1:607088/68‑1 (MQ=255)
 cgcgcgcTTTCGAGCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATg                                              <  1:754106/68‑1 (MQ=255)
              gCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAAt                                               <  1:869278/54‑1 (MQ=255)
              gCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTg                                <  1:129476/69‑1 (MQ=255)
                tAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGAc                              <  1:1028307/69‑1 (MQ=255)
                             cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                 >  1:1994462/1‑69 (MQ=255)
                             cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                 >  1:76995/1‑69 (MQ=255)
                             cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                 >  1:3330162/1‑69 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCTGc                <  1:2089309/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:206809/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:2063352/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:2316279/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:2436211/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:1501337/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:609152/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:95149/68‑1 (MQ=255)
                               gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                <  1:1058932/68‑1 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:2584629/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:596683/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:2741410/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:3393175/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:3419762/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:2233868/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:967569/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:2444661/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1196785/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1333099/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1411940/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:157464/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1623864/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1735348/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1792374/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:1931039/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg     >  1:2457177/1‑67 (MQ=255)
                                           tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCg      >  1:863642/1‑66 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc             >  1:1434650/1‑58 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc             >  1:1430251/1‑58 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc             >  1:1033123/1‑58 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTa  >  1:3093535/1‑69 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTa  >  1:1316419/1‑69 (MQ=255)
                                            cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTAACCGCTGGTAAGCGGCTGc             >  1:522708/1‑58 (MQ=255)
                                                     |                                                           
GCGCGCGCTTTCGAGCTAACGACAAAGTTCAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTA  >  minE/2517672‑2517784

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