Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,517,761 G→A 61.4% T438T (ACC→ACT 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,7610GA61.4% 18.1 / 26.5 44T438T (ACC→ACTspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base A (18/9);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGACC  >  minE/2517701‑2517823
                                                            |                                                              
cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                                                        >  1:76995/1‑69 (MQ=255)
cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                                                        >  1:1994462/1‑69 (MQ=255)
cAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCgg                                                        >  1:3330162/1‑69 (MQ=255)
  gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCTGc                                                       <  1:2089309/68‑1 (MQ=255)
  gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                                                       <  1:1058932/68‑1 (MQ=255)
  gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                                                       <  1:2316279/68‑1 (MQ=255)
  gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                                                       <  1:206809/68‑1 (MQ=255)
  gCCACGCGGCATTCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGc                                                       <  1:2063352/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)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:2741410/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:2584629/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:2457177/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:2444661/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:3393175/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:3419762/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:596683/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:967569/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:1931039/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:1792374/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:1735348/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:1623864/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:157464/1‑67 (MQ=255)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCgg                                            >  1:2233868/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)
              tCGGGCGAGCACCCGGTGCGGTAATGATTTCAACGGTTTGACCGCTAGTAAGCGGCTGCGACAGCg                                             >  1:863642/1‑66 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc                                                    >  1:1033123/1‑58 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc                                                    >  1:1430251/1‑58 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGc                                                    >  1:1434650/1‑58 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTa                                         >  1:1316419/1‑69 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTa                                         >  1:3093535/1‑69 (MQ=255)
               cGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTAACCGCTGGTAAGCGGCTGc                                                    >  1:522708/1‑58 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:3238019/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:394693/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:506812/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:1336057/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:3191717/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:746503/68‑1 (MQ=255)
                                                       ccGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:2922696/68‑1 (MQ=255)
                                                        cGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGAcc  <  1:2931233/67‑1 (MQ=255)
                                                            |                                                              
CAGCCAAGCGGCATTCGGGCGAGCGCCCGGAGCGGTAATGATTTCAACGGTTTGACCGCTGGTAAGCGGCTGCGACAGCGGGTAAGGCTGGCGGTCAACGCGTGCGCCCACGCAGGCATGACC  >  minE/2517701‑2517823

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: